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INTRODUCTION 

Nature  of  the  Problem.  Contamination  of  the  observed  electroencephalogram  (EEG)  by 
physiological  artifacts  such  as  eye  movements  and  blinks  (electrooculogram,  or  EOG)  is  a  classical 
problem  in  electrophysiological  studies.  EOG  artifacts  are  a  major  impediment  in  the  recording  and 
analysis  of  the  EEG;  the  problem  was  first  reported  as  early  as  1941  and  reports  and  attempted 
solutions  continue  to  be  published  as  recently  as  this  year.  The  EEG  (both  steady  state  and  evoked 
potential)  is  an  important  tool  in  the  diagnosis  of  neurological  dysfunction,  such  as  epilepsy, 
cerebrovascular  trauma  and  brain  tumor,  as  well  as  in  sleep  studies  and  determining  workload, 
mental  state,  and  limitations  of  sensory  information  processing.  Our  approach  uses  a  direct 
interrogation  (signal  injection)  technique  which  is  substantially  different  from  the  a  posteriori 
techniques  published  by  others.  Direct  interrogation  permits  on-line,  virtually  real  time  correction 
of  eye  movement  artifacts  on  the  observed  EEG.  We  conducted  a  research  program  to  add  the 
capability  of  blink  artifact  removal  to  our  current  artifact  rejection  technique  and  implement  this 
in  a  hardware  device.  The  capability  of  completely  correcting  the  EEG  for  EOG  artifacts  offers  a 
number  of  significant  benefits.  It  simplifies  the  acquisition  and  analysis  of  the  steady  state  EEG  and 
evoked  potentials  (EP)  in  clinical  studies.  In  addition,  it  simplifies  the  development  of  algorithms 
for  machine  analysis  of  the  EEG  and  EP;  such  analysis  has  been  thwarted  by  the  presence  of  EOG 
artifacts. 

Background.  A  biological  signal  is  propagated  from  its  site  of  origin  to  a  site  of  measu’ement 
through  a  medium  that  can,  in  principle,  be  described  by  a  transfer  function.  A  transfer  function 
is  merely  a  mathematical  construct  describing  the  relationship  between  the  "input"  and  the  "output" 
of  the  medium.  This  description  can  be  accomplished  equivalently  in  the  time  domain  or,  after 
Fourier  transformation,  in  the  frequency  domain.  Theoretically,  the  transfer  function  can  be 
completely  determined  by  measuring  the  response  of  the  medium  to  a  unit  impulse;  the  unit  impulse 
response  is  the  inverse  transform  of  the  transfer  function.  Convolving  the  unit  impulse  response 
with  any  input  always  yields  the  output  (Van  Valkenburg,  1964). 

It  is  important  to  note  that  unique  transfer  functions  can  only  be  defined  for  linear  systems. 
All  naturally-occurring  systems  are,  in  the  final  theoretical  analysis,  distributed-parameter,  non¬ 
stationary,  quantized,  stochastic,  nonlinear  systems.  However,  it  has  been  repeatedly  proven  that 
modelling  naturally-occurring  systems  as  nearly  lumped-parameter,  stationary,  continuous, 
deterministic,  linear  systems  can  yield  practically  useful  results  (Riggs,  1970). 

Observed  biological  signals,  such  as  EEG,  are  composites  of  the  signal  of  interest  and  other 
unwanted,  but  nevertheless,  real  and  physiological  signals.  These  unwanted  signals  are  traditionally 
termed  "artifacts".  In  theory,  there  are  three  ways  of  eliminating  unwanted  signals  -  prevention, 
traditional  filtering,  and  filtering  after  estimation  of  the  transfer  function. 

Prevention  is  possible  for  many,  but  not  all,  unwanted  signals.  Technical  artifacts,  such  as 
60  Hz  "hum"  and  electrochemical  effects  at  the  electrode  -  tissue  interface,  are  preventable. 
Certainly,  artifact  prevention  is  the  ideal.  However,  physiological  artifacts  are  not  preventable;  it 
is  clearly  "undesirable"  to  prevent  eye  movements  in  a  visual  tracking  task  OR  to  stop  the  maternal 
heart  during  an  attempt  to  obtain  the  fetal  ECG! 

Traditional  filtering  is  useful  for  many,  but  not  all,  unwanted  signals.  When  there  is  good 
frequency  separation  between  the  signal  of  interest  and  the  unwanted  signals,  filters  such  as 
Butterworth,  Tschebychev,  Bessel,  and  Cauer  are  usef ul.  However,  in  the  absence  of  good  f requency 
separation  (or  when  the  bandwidth  of  the  signal  of  interest  overlaps  the  bandwidth  of  the  unwanted 
signals),  traditional  filtering  results  in  loss  of  data  and  signal  distortion.  Steep  roll-off  of  these 
filters  (higher  order  forms)  results  in  significant  phase  distortion,  which  is  often  overlooked  in  visual 
analysis  (Johnson  et  al,  1979).  Often,  EEG  data  which  is  deemed  contaminated  based  on  some 
subjective  or  pre-programmed  criterion  is  simply  rejected  (Gratton  et  al,  1983)  or  masked  (Barlow, 
1985). 

The  problem  of  determining  the  state  of  a  system  from  noisy  measurements  is  called 
estimation  (or  filtering).  With  a  state-space  approach,  the  dynamical  system  is  modeled  by  a  finite- 
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dimensional  Markov  process;  the  conditional  probability  density  function  of  the  state  embodies  all 
the  information,  which  is  available  from  the  measurements  (Jazwinski,  1970).  All  estimates  of  the 
state  can  be  constructed  from  the  density  function,  allowing  formulation  of  linear  and  nonlinear 
filters  and  predictors.  Removal  of  eye  movement  "artifacts"  from  observed  EEG  measurements  is 
an  example  of  such  a  problem.  This  is  not  merely  a  problem  of  theoretical  interest;  eye  movements 
and  eye  blinks  are  a  permanent  source  of  serious  unwanted  signals  in  the  measurement  and  analysis 
of  the  electroencephalogram  (Lyman,  1941;  Case,  1959;  Corby  &  Kopell,  1972;  Girton  &  Kamiya, 
1973;  Matsuo  et  al,  1975;  Gevins  et  al,  1977;  Whitten  et  al,  1978;  Barlow  &  Remond,  1981;  Verleger 
et  al,  1982;  Gratton  et  al,  1983;  Fortgens  &  De  Bruin,  1983;  Woestenburg  et  al  1983;  Elbert  et  al, 
1985)  and  are  especially  troublesome  in  event-related  brain  potential  measurements  (e.g.  CNV  and 
P300)  (Hillyard  &  Galambos,  1970;  Wasman  et  al,  1970;  Girton  &  Kamiya,  1973). 

Before  we  ask  the  question  "What  is  the  transfer  function  which  describes  the  coupling  of  the 
EOG  to  the  observed  EEG?",  let  us  first  review  the  biophysics  underlying  the  process. 

In  electrodynamic  terms,  the  eye  can  be  modeled  as  a  dipole  (Mowrer  et  al,  1935;  Barry  & 
Jones,  1965).  In  this  representation,  the  corneo-retinal  potential  difference  is  the  result  of  charge 
separation,  with  the  corneal  aspect  being  positive  and  the  retinal  aspect  being  negative.  The  corneo- 
retinal  potential  is  known  to  vary  considerably  between  individuals  (Shackel,  1967)  as  well  as  within 
individuals  as  a  function  of  illumination  (Rubin  &  Walls,  1969)  and  time  (Shackel  &  Davis,  1960). 
It  is  further  assumed  (Fortgens  and  De  Bruin,  1983;  Elbert  et  al,  1985)  that  under  normal  conditions 
both  eyes  move  conjugately.  The  motion  of  these  linked  dipoles  creates  potentials,  which  are 
observable  at  distant  EEG  recording  sites.  In  addition  to  the  potentials  created  by  this  dipole 
motion,  potentials  are  also  generated  by  the  eyelid,  acting  like  a  sliding  electrode,  picking  up  positive 
potential  moving  across  the  positively  charged  corneal  surface  (Matsuo  et  al,  1975).  The  resultant 
change  in  charge  distribution  caused  by  closure  of  the  eyelid  can  also  be  described  by  a  change  in 
dipole  moment.  Thus,  the  EOG  can  be  described  as  an  electrical  potential  resulting  from  a  change 
in  ocular  dipole  moment  (Girton  &  Kamiya,  1973;  Elbert  et  al,  1985). 

The  EOG  potential  is  propagated  through  the  medium  to  sites  all  over  the  body;  this  includes 
the  head,  which  may  be  modeled  as  a  four  layer  sphere  of  different  conductivities  (Cuffin  &  Cohen, 
1979).  The  relationship  between  the  EOG  potential  (measured,  say,  near  its  source)  and  the 
propagated  EOG  potential,  at  some  distant  site,  is  completely  described  by  the  transfer  function  of 
the  medium  (assumed  linear).  Numerous  workers  have  attempted  to  exploit  this  relationship  in  order 
to  remove  the  unwanted  EOG  from  the  observed  EEG. 

In  its  simplest  form,  the  transfer  function  might  be  assumed  to  be  constant  and  unity. 
However,  if  we  merely  subtract  the  measured  EOG  from  an  EEG  measured  at  distant  scalp  sites,  it 
is  obvious  that  this  "correction"  will  yield  an  erroneous  estimate  of  the  real  EEG  (Gratton  et  al, 
1983).  The  transfer  function  could  be  assumed  to  be  constant,  but  not  unity  (Barlow  &  Remond, 
1981);  however,  this  would  not  take  into  consideration  the  known  dependence  on  the  distance  to  the 
EEG  electrode  site  from  the  eye  (Girton  &  Kamiya,  1973).  The  transfer  function  could  be  assumed 
to  be  distance  dependent;  however,  this  does  not  take  into  consideration  the  known  frequency  and 
phase  angle  dependence  of  propagated  volume  conductor  potentials  (Gevins  et  al,  1977;  Whitton  et 
al,  1978;  Woestenburg  et  al,  1983;  Elbert  et  al,  1985). 

From  a  theoretical  point  of  view,  the  transfer  function  describing  the  propagation  of  the 
EOG  potential  (v;  [t])  through  the  medium  to  the  distant  recording  site  (v„  [t]),  is  a  function 
h[t]=h[d,A,f,^,t],  such  that: 

VoM  =f  Vi  [fi]  h[t-/i]  d(i 
where:  J 

Vj  =  input  potential 
v0=  output  potential 
d  =  distance 
A  =  amplitude 
f  =  frequency 
<t>  =  phase  angle 
t  =  time 
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From  a  practical  point  of  view,  distance  dependence  of  the  transfer  function  can  be  ignored 
for  specific,  fixed  electrode  sites  (such  as  in  a  single  recording  session).  Frequency  dependence 
(Gevins  et  al,  1977;  Whitton  et  al,  1978;  Elbert  et  al,  1985)  and  phase  angle  dependence  (Whitton  et 
al,  1978;  Woestenburg  et  al,  1983)  cannot  be  ignored.  We  can  find  no  evidence  in  the  literature  for 
amplitude  dependence;  yet,  this  does  not  mean  that  it  can  be  arbitrarily  ignored  without 
investigation.  Finally,  time  dependence  of  the  transfer  function  should  not  be  ignored  (except 
possibly  in  very  short  duration  recording  sessions),  since  it  is  a  fundamental  premise  that  biological 
systems  change  with  time.  An  intuitive  illustration  of  this  might  be  the  temporally-dependent 
impedance  changes  resulting  from  perspiration. 

An  analysis  of  the  published  literature  clearly  indicates  that  thinking  in  this  area  has  been 
slowly  evolving  to  the  aforementioned  full  theoretical  form  of  the  transfer  function  and  the 
parameters  it  depends  upon.  In  fact,  based  on  the  1985  work  of  Elbert  et  al,  it  appears  that  direct 
interrogation  (as  described  later  in  this  proposal)  is  the  next  logical  area  of  investigation  for  the 
removal  of  the  EOG  from  the  EEG. 

Numerous  attempts  have  been  made  to  estimate  the  transfer  functions  of  unwanted  biological 
signals.  Bergveld  &  Meijer  (1981)  have  reported  a  technique  for  removing  the  maternal  ECG  from 
abdominal  electrocardiograms,  in  order  to  obtain  a  fetal  ECG  as  well  as  a  technique  for  determining 
the  ideal  electrode  position  (Meijer  &  Bergveld,  1981).  They  postulate  a  transfer  function  composed 
of  the  linear  combination  of  three  independent  observation  sites  and  attempt  to  estimate  the 
coefficients  of  this  linear  combination.  Johnson  et  al  (1979)  have  reported  a  technique  for  removal 
of  muscle  artifact  from  the  electroencephalogram.  They  formulate  a  nonlinear  estimator  (filter) 
based  on  an  a  priori  model  of  the  EEG  (represented  as  the  superposition  of  four  lightly  damped 
oscillators,  operating  in  the  alpha,  beta,  theta,  and  delta  bands,  driven  by  independent  white 
Gaussian  noises)  and  an  a  priori  model  of  the  muscle  artifact  (represented  by  the  superposition  of 
"action  potentials"  of  three  different  durations  generated  as  impulse  responses  of  three  linear  systems 
driven  by  independent  Poisson  processes).  Techniques  for  removing  the  EOG  "artifact"  from  the 
EEG  have  been  reported  by  Verleger  et  al  (1982),  Gratton  et  al  (1983),  Woestenburg  et  al  (1983),  and 
Elbert  et  al  (1985). 

Verleger  et  al  report  "completely  correcting  for  blink  effects",  but  only  partial  correction  of 
eye  movement  artifact;  this  is  in  contrast  to  Weerts  &  Lang  (1973)  who  "presumably  removed  the 
eye  movement  effect  correctly,  but  overcompensated  for  the  blink  effect"  (Verleger  et  al,  1982). 
They  use  a  regression  approach  consisting  of: 

a.  identifying  maximum  variance  EOG  segments; 

b.  estimating  a  linear  regression  coefficient; 

c.  estimating  a  general  transmission  rate; 

d.  correcting  the  EOG  for  DC  bias;  and 

e.  subtracting  the  weighted  EOG  from  the  observed  EEG. 

Gratton  et  al  (1983)  use  a  somewhat  different  approach.  Their  procedure  consists  of: 

a.  estimating  correction  factors  derived  from  EOG  and 
EEG  data  obtained  during,  rather  than  before,  the 
experiment; 

b.  estimating  separate  correction  factors  for  blinks  and 
eye  movements; 

c.  removing  event-related  EOG  and  EEG  activity  from 
the  data;  and 

d.  subtracting  the  weighted  EOG  from  the  observed  EEG. 

They  state  that  their  approach  has  six  clear  advantages:  it  distinguishes  between  blink  and  eye 
movement  artifact;  it  provides  corrections  that  are  insensitive  to  stimulus-locked  activity;  it  retains 
all  data  for  use  in  subsequent  analyses;  it  does  not  require  special  data  collection;  the  subjects  need 
not  control  or  minimize  eye  movements;  and,  the  estimate  is  based  on  a  large  sample,  rather  than  a 
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few  data  obtained  from  a  few  prescribed  eye  movements.  They  also  properly  point  out  that  "noise" 
in  the  measured  EOG  may  significantly  alter  the  magnitude  of  the  estimated  correction  factor. 

Woestenburg  et  al  (1983)  report  a  technique  for  removing  the  eye  movement  artifact  from 
the  EEG  by  regression  analysis  in  the  frequency  domain.  They  explicitly  recognize  and  demonstrate 
that  the  transfer  of  eye  movement  activity  to  EEG  can  have  frequency  dependent  amplitude  and 
phase  characteristics  and  they  attempt  to  determine  the  transfer  function.  They  assume  that  the 
medium  is  passive  and  constant  and  that  there  is  no  linear  correlation  between  EOG  and  EEG 
activity.  Furthermore,  they  state  that  "a  successful  method  for  removing  the  EOG  artifact  from  the 
EEG  should  be  able  to  handle  the  following  phenomena: 

a.  Transfer  from  EOG  on  EEG  is  frequency  dependent. 

Some  frequencies  may  be  attenuated  more  than  other 
frequencies. 

b.  The  EOG  artifact  as  measured  at  the  scalp  can  be 
distorted  by  phase-shifts. 

c.  Both  vertical  and  horizontal  eye  movements  may 
contribute  to  the  artifact." 

Woestenburg  et  al  (1983)  applied  their  technique  to  simulated  data  as  well  as  to  real  data.  The 
principal  limitation  of  their  technique  is  that  it  is  an  a  posteriori  approach  typically  requiring  two 
blocks  of  36  complex  visual  stimulus  presentations  and  about  one  hour  of  computer  time  for  data 
analysis. 

Elbert  et  al  (1985)  use  a  biophysical  approach  to  the  theoretical  formulation  of  the 
electrodynamic  equations,  which  allow  a  complete  description  of  the  ocular  influence  in  the  EEG. 
They  separate  the  transfer  function,  describing  the  ocular  influence  in  the  EEG,  into  vertical, 
lateral,  and  radial  components  and  attempt  to  identify  (but  do  not  adequately  support)  the  minimum 
necessary  and  sufficient  EOG  electrodes  and  their  anatomical  positions.  Elbert  et  al  explicitly 
recognize  the  frequency  dependence  of  the  transfer  function;  they  report  the  form  of  the  vertical 
component  (g[w,CJ)  as  a  function  of  radial  frequency  (w  =  2zrf )  as  measured  at  the  Cr.  They  both 
report  theoretical  and  empirical  forms.  There  are  two  empirical  forms  reported.  One  form, 
attributed  to  Gasser  et  al,  is  derived  from  naturally  occurring  ocular  artifacts.  The  other  form  was 
derived  following  application  of  an  (unspecified)  artificial  drive  signal  to  the  EOG  electrodes. 

The  application  of  this  artificial  drive  signal,  by  Elbert  et  al,  forms  the  published 
"springboard"  of  our  research  efforts.  The  artificial  drive  signal,  applied  to  the  EOG  electrodes,  is 
an  example  of  direct  interrogation  of  the  biological  system  under  consideration.  In  keeping  with  the 
theoretical  approach  to  determining  the  transfer  function,  it  allows  us  to  apply  a  "unit  impulse",  so 
as  to  completely  describe  the  real  transfer  function.  Judicious  selection  of  an  externally  applied 
drive  signal,  when  properly  utilized,  can  be  a  safe,  effective,  and  noninvasive  means  of  determining 
the  transfer  function  of  the  ocular  influence  on  the  EEG.  An  artificial  drive  signal  has  already  been 
applied  by  Elbert  et  al  (1985)  and  by  us  (unpublished,  1985  and  Falk  et  al,  1987).  Sullivan  (1965) 
reported  use  of  a  40  KHz  drive  signal  for  measuring  impedance  in  order  to  determine  the  direction 
of  the  ocular  dipole. 

All  the  previously  cited  literature  (with  the  exception  of  Elbert  et  al)  attempt  to  determine 
the  transfer  function  (correction  factor,  weighting  factor,  regression  coefficient,  etc.)  through  the 
use  of  naturally  occurring  ocular  motions.  Since  there  are,  potentially,  an  infinite  number  of 
different  ocular  motions,  selection  of  specific  motions  (Weerts  &  Lang,  1973;  Verleger  et  al,  1982; 
Fortgens  &  De  Bruin,  1983)  obviously  lacks  generality  and  completeness.  The  work  of  Woestenburg 
et  al  (1983),  and  then  Gratton  et  al  (1983),  begins  to  circumvent  this  problem  by  basing  the  estimate 
on  a  large  sample,  rather  than  a  few  data  obtained  from  a  few  prescribed  eye  movements.  But  even 
this  approach  does  not  fully  address  the  problem.  Our  approach  is  to  apply  an  external  drive  signal 
which  describes  all  possible  ocular  motions;  these  ocular  motions  are  merely  an  electrical  signature 
composed  of  particular  ap'plitudes  at  particular  frequencies  with  particular  phase  relations.  In  fact, 
because  the  biopotentials  generated  by  ocular  motion  are  not  unbounded,  the  EOG  does  NOT 
contain  all  possible  amplitudes  and  frequencies;  the  EOG  is  constrained  to  frequencies  below,  say 
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for  example,  30  Hz  and  to  amplitudes  below,  say  for  example,  5  mV.  Therefore,  practically,  the 
"unit  impulse"  required  to  theoretically  determine  the  transfer  function  need  not  be  an  impulse 
input,  $(t);  instead,  it  can  be  a  relatively  short  time  duration  pulse  whose  frequency  transform 
includes  those  frequencies  of  interest. 

Purpose  of  the  Present  Work.  The  purposes  of  this  research  study  were  twofold:  first,  to 
refine  our  existing  mathematical  technique,  and  second,  to  implement  it  in  a  portable,  battery- 
operated  twelve  channel  device. 

Methods  of  Approach  -  Mathematical.  The  mathematical  technique  that  allows  us  to  remove 
the  eye  movements  from  the  on-going  EEG  is  called  the  direct  interrogation  technique.  Three  basic 
assumptions  are  made  in  order  to  utilize  this  technique.  We  assume  that  the  eye  movement  signal 
propagates  only  on  the  surface  of  the  head  to  the  distant  EEG  sites.  Since  the  skull  is  approximately 
eighty  times  the  resistivity  of  the  scalp,  the  path  of  least  resistance  is  the  scalp.  Depth  electrode 
studies  have  been  conducted  and  there  was  no  evidence  of  EOG  artifact  in  the  EEG  (Cooper,  1971). 
We  also  assume  that  the  medium  is  linear  (thus  the  theory  of  superposition  holds)  and  the  medium 
is  non-dispersive  (no  shift  in  frequency).  We  have  tested  both  these  assumptions  and  we  find  them 
to  be  true.  With  this  as  our  base,  we  can  model  this  system  as  an  input  (EOG),  an  output  (EOG 
artifact  on  the  EEG),  and  a  medium  (scalp)  and  its  transfer  function. 

Before  we  discuss  the  models  of  eye  movement  and  eye  blink,  we  must  discuss  some 
terminology  regarding  the  ocular  dipole.  An  electric  dipole  is  an  electric  potential  source  arising 
from  the  separation  of  equal  and  opposite  charges  and  resulting  in  an  electric  field  whose  magnitude 
is  nonzero  at  all  points  in  space  except  those  equidistant  from  both  charges.  These  equidistant  points 
define  a  unique  zero- potential  plane  orthogonal  to  the  line  connecting  both  charges.  The  ocular 
dipole  is  an  electric  dipole  with  the  positive  chaige  on  the  cornea  and  the  negative  charge  on  the 
retina.  The  conjugate  eye  dipole  pair  is  comprised  of  the  two  linked  ocular  dipoles  that  move  in 
parallel.  The  surface  image  of  a  dipole  is  that  portion  of  the  electric  field  residing  on  a  surface 
transecting  the  three  dimensional  dipole  electric  field.  The  image  of  the  zero-potential  plane  on  the 
surface  is  a  zero- potential  line.  The  surface  image  of  the  ocular  dipole  is  the  image  on  that  surface 
defined  by  the  skin  on  th“  h»»ad  (including  the  face).  The  zero-potential  line  forms  an  angle  <f>  with 
the  x-axis  of  our  geometrical  coordinate  system.  A  direct  interrogation  stimulus  dipole  or  "surface 
stimulus  dipole"  is  the  electric  source  resulting  from  the  application  of  two  spaced  surface  elec¬ 
trodes  driven  by  a  floating  voltage  source  (a  floating  battery). 

We  made  the  following  explicit  assumptions:  (a)  the  EOG  signal  reaches  the  EEG  recording 
site  via  surface  propagation  (propagation  by  other  means  is  negligible);  (b)  the  medium  is  passive 
and  constant  (over  relatively  short  time  periods);  (c)  the  principle  of  superposition  holds  (the  system 
is  linear  or  nearly  linear)  and  a  unique  transfer  function  does  exist;  (d)  the  medium  is  non- 
dispersive  (frequencies  don’t  change  during  propagation);  and  (e)  our  mathematical  model  properly 
represents  the  electrodynamic  behavior  of  the  conjugate  eye  dipole  pair.  While  further  investigation 
is  required,  we  presently  believe  that  no  other  implicit  assumptions  have  been  made. 

CORRection  of  the  observed  EEG  for  EOG  artifacts  (corr  VEEG  [w])  is  accomplished  in  the 
frequency  domain  and  is  based  on  (a)  measurement  of  the  OBServed  EEG  (obaVEEG  [w])  and 
OBServed  EOG  (ob*VEOG[w])i  (*>)  measurement  of  the  system  response  to  STIMulation  ("""VeegM 
&  *tim  VEOG  M)  for  direct  interrogation,  (c)  a  mathematical  model  that  describes  the  electrody¬ 
namic  behavior  of  the  system  for  Theoretical  Eye  Movements  (tem  VEEg  [w]  &  tem  mV eog  (wJ)  an^ 
Theoretical  Direct  Interrogation  (tdlVEEC{w]  &  td>VEodw]K  and  (d)  measurement  of  the  CALibration 
of  each  recording  channel  (cal  VEEg  [w]  &  cal  VEOg  [w]).  The  mathematical  derivation  is  summa¬ 
rized  here. 
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The  formula  for  implementing  the  EEG  correction,  on  a  frequency  per  frequency  basis,  is: 
c°"  VEEG  [w]  =  °b‘  VEEG  M  -  (SxD/G)  ob8  Veog  M 

where: 

s  _  *tim  vEEG  [w]  +  itim  VEOg  [w]  (using  20  /iA  stimulus  pulse) 

G  =  cal  VEEG  [u>]  +  cal  VEOg  M  (using  1  mV  calibration  pulse) 

D  =  {tem  V EEg  [^]  +  tem  VEqg  lw])  +  (tdl  V EEg  (w)  +  td'  Veog  tw]) 

S  is  a  measure  of  the  system  response  to  direct  interrogation  and  is  the  ratio  of  the  signals  observed 
at  the  EEG  and  EOG  recording  sites;  it  is  the  putative  transfer  function.  G  is  a  measure  of  the 
discrepancy  between  the  recording  channels  and  is  the  ratio  of  the  calibration  signals  observed  at  the 
EEG  and  EOG  recording  sites;  G  would  not  be  necessary,  if  and  only  if  the  recording  channels  were 
absolutely  identical.  D  is  a  geometrical  correction  factor  that  interrelates  the  theoretical  electrody¬ 
namic  behavior  of  the  (non-coilocated)  direct  interrogation  stimulus  dipoles  and  the  ocular  dipoles; 
it  is,  in  fact,  our  mathematical  model.  It  must  contain  both  magnitude  and  phase  information,  so  it 
has  the  form: 


D  =  D’e* 

where  D’  describes  the  magnitude  correction  due  to  geometry  and  £  describes  the  phase  correction 
due  to  geometry.  The  geometrical  correction  factor  D  would  not  be  necessary,  if  and  only  if  the 
direct  interrogation  stimulus  dipole  exactly  and  completely  emulated  the  ocular  dipoles  geometrically 
and  electrodynamically. 

D’  was  derived  by  obtaining  the  general  solution  of  the  general  differential  equation  that 
describes  the  propagation  of  a  potential  generated  by  any  source.  The  general  solution  was 
constrained  to  model  a  dipole  source.  Using  this  equation,  the  conjugate  eye  dipoles  were  resolved 
into  a  single  equivalent  theoretical  source  located  at  the  origin  of  our  selected  coordinate  system. 
Similarly,  by  coordinate  transformation,  the  stimulus  dipoles  were  converted  to  an  equivalent 
theoretical  source  also  located  at  the  origin  of  our  coordinate  system.  With  these  two  source 
equations,  the  magnitude  relationship  of  the  signals  expected  at  the  EEG  and  EOG  recording  sites 
(as  a  result  of  eye  movements  versus  surface  dipole  stimulation)  was  computed.  This  permits 
computation  of  the  magnitude  portion  of  the  geometrical  correction  factor;  it  is  used  to  correct  the 
empirical  transfer  function  (found  by  direct  interrogation  stimulation)  for  the  difference  in 
geometry  between  the  stimulus  dipoles  and  oc**'ar  dipoles. 

Phase  changes  due  to  propagation  through  the  medium  and  this  information  is  contained  in 
the  empirical  transfer  function  obtained  by  direct  interrogation.  Additionally,  there  is  a  relative 
phase  shift  between  the  EEG  and  EOG  recording  sites.  It  is  due  solely  to  the  changing  geometric 
orientation  of  the  isopotential  lines  caused  by  rotation  of  the  surface  image  of  the  ocular  dipole. 
This  information  is  not  contained  in  the  direct  interrogation  data  and  must  be  independently 
corrected.  The  equation  describing  the  single  equivalent  theoretical  source  of  the  conjugate  eye 
dipole  is  a  function  of  the  angle  of  rotation  <f>  of  the  surface  image  of  the  ocular  dipole  pair. 
Differentiation  of  this  equation  with  respect  to  <f>  yields  an  equation  describing  the  change  in 
potential  at  an  EOG  electrode  due  to  a  change  in  0.  When  the  change  in  potential  with  respect  to 
<f>  is  zero,  the  potential  is  at  an  extremum  (maximum  or  minimum)  and  the  corresponding  </>,  at  a 
particular  electrode  site,  can  be  computed.  This  value  of  <t>  is  the  value  of  the  angle  of  rotation  that 
creates  an  extremum  at  the  particular  electrode  site  under  consideration.  It  will  have  different 
values  for  different  electrode  sites.  The  geometrically-dependent  relative  phase  shift  between  an 
arbitrary  pair  of  electrode  sites  is  the  difference  of  their  corresponding  <f>' s.  A  change  in  <f>  can  not 
be  determined  from  one  EOG  electrode;  in  general,  an  orthogonal  pair  is  preferred. 
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Our  mathematical  technique  can  be  summarized  as  follows.  Integral  to  our  technique  are  the 
following  four  (4)  explicit  assumptions: 

a.  the  EOG  artifact  on  the  EEG  is  the  result  of  an  electrodynamic  process,  arising  from  the 
movement  of  the  eye  dipoles  and  from  the  eyelids  across  their  surface  (Elbert  et  al,  1985); 

b.  the  EOG  artifact  reaches  the  EEG  recording  site  primarily  via  surface  propagation 
(propagation  by  other  means  is  negligible)  (Cooper  et  al,  1965,  1971;  Cuffin  and  Cohen, 
1979); 

c.  the  surface  propagation  medium  is  passive,  linear,  and  non-dispersive  (over  relatively 
short  time  periods);  thus,  a  unique  transfer  function  exists  -  this  was  shown  in  our  feasibility 
demonstration;  and 

d.  all  possible  eye  movements  and  blinks  are  completely  described  by  their  Fourier 
components,  and  these  consist  of  a  bounded  set  of  frequencies,  amplitudes,  and  phases. 

Therefore,  EOG  propagation  between  the  site  of  EOG  generation  and  the  EEG  electrodes  can  be 
characterized  by  a  transfer  function;  the  transfer  function  in  turn  can  be  characterized  by  injecting 
a  signal  at  the  EOG  generation  site  and  recording  the  resultant  signal  at  the  EEG  electrodes  (direct 
interrogation). 


This  method  of  rejecting 
mathematically  expressed  as: 

ocular  motion  artifacts  on  the  EEG  recording  can 

EEG^t)  = 

EEG^t)  - 

I FT  (EOGc(s)xS(s)yG(s)xD) 

where:  EEG^t) 

= 

Corrected  EEG  (time  domain) 

EEG'Xt) 

= 

Observed  EEG  (time  domain) 

I  FT 

= 

Inverse  Fourier  Transform 

EOG^s) 

= 

Observed  EOG  (frequency  domain) 

S(s) 

= 

Transfer  function  (frequency  domain) 

G(s) 

= 

Channel  response  correction  factor  (frequency  domain) 

D 

= 

Geometric  correction  factor 

Methods  of  Approach  -  Engineering.  The  mathematical  technique  is  an  algorithm  for 
removing  the  unwanted  influence  of  the  EOG  on  the  observed  EEG.  This  can  be  implemented  in 
hardware  by  constructing  a  microprocessor-based  device  which  can  be  programmed  to  execute  this 
algorithm.  Amplifiers  and  filters  are  used  to  condition  biopotential  signals  which  can  then  be 
digitized  and  processed.  These  processed  data  can  be  once  again  converted  to  analog  signals  for 
display  and  recording.  The  digital  processing  time  will  introduce  a  finite  delay  due  to  the  time 
required  for  the  microprocessor  to  execute  the  necessary  computations.  Standard  engineering 
techniques  permit  implementation  of  the  analog  and  digital  circuitry  in  a  form  that  requires  minimal 
power,  and  thus  can  be  battery  operated. 
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Refinement  of  the  Mathematical  Technique.  In  this  research  study,  we  have  expanded  the 
biophysical  model  to  include  the  blink.  This  yields  a  general  electrodynamic  model  for  both  the 
source  and  the  propagating  electric  field  from  the  eye  for  all  possible  eye  movements  and  blinks. 

The  Biophysical  Model.  The  transection  of  the  face  across  the  three  dimensional  ocular 
dipole  field  (caused  by  the  corneo-retinal  potential  in  the  eye)  yields  a  surface  image  dipole 
propagating  on  the  scalp.  This  surface  image  dipole  can  be  modelled  to  incorporate  both  the  eye 
movement  and  the  eye  blink.  The  eye  movement  produces  a  symmetric  dipole,  while  the  blink 
produces  an  asymmetric  dipole. 

General  Dipole  Representation.  A  dipole  source,  symmetric  or  asymmetric,  is  the 
superposition  of  two  point  sources  separated  by  a  distance.  The  point  source’s  electric  field 
propagates  as  a  function  of  1/r2.  The  voltage  at  any  point  is  described  by  V=kq/r,  where  k  is 
Boltzmann's  constant,  and  q  is  the  amount  of  charge.  The  surface  image  dipole  is  described  here. 


+  -mq 


-nq 


The  figure  on  the  left  shows  two  point  sources  separated  by  a 
distance  (L).  The  voltage  (V)  appearing  at  point  Q  is  derived  as  follows. 

V  =  kq([m/(r-)Lsin0)]+[-n/(r+)Lsin0)]).  (1) 

Rearranging  Equation  1  yields 

V  =  kq{[(m-n)r+(m+n))Lsin0]/(r2-i  L2sin20)}  (2) 

Since  r»L,  we  can  simplify  Equation  2: 


V  =  n[(kq/r2)fLsin0(a+I)  +  (kq/r)(a- 1 )],  (3) 

where  a  =  m/n. 

As  a  note,  if  a=l  (eye  movement)  and  the  dipole  is  symmetric. 
Equation  3  reduces  to, 

V  =  (nkq/r2)Lsin#  =  Ar'linfl.  (4) 

where  A  =  nrqL. 


Furthermore,  it  is  important  to  note  that  the  zero-potential  line  of  the 
dipole  is  the  x-axis  when  a=l  (sin0=O).  When  a*i,  the  zero-potential  line  becomes  a  circle  described 
by, 

x2+  (y+G)2=  G2,  (5) 

where  G=)L[(a+l)/(a-l)j. 
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Selection  of  a  Facial  Coordinate  System.  In  order  to  spatially  represent  the  ocular  dipoles 
or  the  stimulus  dipole  in  planar  geometry,  we  must  select  a  coordinate  system.  Tiiis  is  shown  in 
the  figure  below.  The  point,  Q,  in  this  figure  represents  an  electrode.  The  subscript  L  is  used  to 
show  reference  to  the  left  eye,  the  subscript  R  is  used  to  reference  the  right  eye.  The  electrode  is 
a  distance  /9j,from  the  left  eye  and  /?afrom  the  right  eye.  Reference  to  the  stimulation  dipole  is 
indicated  by  S.  Our  facial  coordinate  system  has  its  origin  at  the  geometric  center  of  the  two  eyes 
(which  are  separated  by  a  distance  z).  The  stimulus  dipole  is  located  a  distance  h  above  the  origin. 
The  electrode  is  a  distance  r  from  the  origin  of  the  coordinate  system.  The  stimulation  dipole  lies 
on  the  y-axis  as  the  eyes  lie  on  the  x-axis. 


Spatial  Resolution  of  the  Bi-ocular  Dipoles.  Using  our  coordinate  system  and  dipole 
representation  described  above,  we  will  model  the  two  ocular  dipoles  as  one  complex  mathematical 
function  based  at  the  origin  of  the  coordinate  system. 

From  the  law  of  cosines: 

/?L=  [r2  +  iZ2  +  rZcos#]*  (6) 

/3r=  (r2  +  iZ2  -  rZcos#]*  (7) 

From  the  law  of  sines: 

7l=  ~<f>  +  sin_1[(r//?|)sin0]  (8) 

7r  =  *-<f>  -  sin^Kr/Z^sinfl]  (9) 

where:  <f>  =  angle  of  the  ocular  dipoles  (zero-potential  line). 
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We  now  substitute  these  identities  (Eqns.  6-9)  into  the  dipole  equation  defined  earlier  (Eqn.  3)  and 
sum  the  two  ocular  dipoles  to  result  in  one  function.  We  obtain  this,  in  a  general  form,  for  any 
electrode: 

V  =  *(<*+ 1  )A{[/9l'2[-(  1  -7L2 )1/2]sin<£+7Lcos^]+[/?R 2[(  1  -7R2)1/2]sin0+7Rcos^]}+(A/L)(a- 1  ){/?L"  1+0r  1 ), 

(10) 


where: 


P\,~  [(r)2+(iZ)2+rZcos 6]\'f 
/?R=  [(r)2+(iZ)2-rZcos0]1/2 
7l=  (r//?i)sin  6 
7r=  (r//?H)sine 


Spatial  Representation  of  the  Stimulus  Dipoles.  We  can  similarly  describe  the  representation 
of  the  stimulus  dipole  in  our  new  coordinate  system  as  we  have  described  the  ocular  dipoles  above. 
Using  Equation  4: 


V  =  B{<7~2[f7sint/>-(l-f32)1/2cosV>]}  (11) 

where: 

a  =  [r2  +  h2  -  2rhsin0]1/2 
r i  =  (-r/a)cos0 

\l>  =  angle  of  the  stimulus  dipole  (an  analog  of  <f>) 

Equations  10  &  1 1  and  the  associated  identities  are  the  basis  of  the  mathematical  model  which  will 
be  used  in  the  calculation  of  the  transfer  function. 

Mathematical  Relationship  Between  Resultant  Ocular  Dipoles  &  Stimulus  Dipoles.  We  can 
now  use  our  basic  equations  and  our  coordinate  system  to  correct  the  putative  transfer  function 
measured  by  surface  dipoles.  The  EOG  artifact  correction  equation  in  the  frequency  domain  is: 

corr  Veeg  =  °b"  ^eeg  -  ob“  VEOG  [Seeg  /Seog  1  d  (12) 

where  D  is  the  geometrical  correction  factor  between  the  stimulus  dipole  and  the  ocular  dipoles.  S 
denotes  the  surface  dipole  stimulation  response  and  obs  V  denotes  the  naturally  occurring  response. 
The  subscript  "EEG"  and  "EOG"  refer  to  the  electrode  recording  the  response;  the  superscripts  "corr" 
and  "obs"  refer  to  the  corrected  and  observed  potential,  respectively. 
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VVe  can  obtain  D  by  manipulation  of  the  equations  described  above. 

{)(<*+ 1  )[[^l  2  H 1  -'YL)1/2]sin^+'rLCOS<^]+[/9R'2  [(1  -7R2)1/2]sin^+7Rcos^]]+((a- 1  1  J)eeg 

[-(i-TfL2)1^5]sin^+TfLCOS^]+[^R5'i(l-'7R)1^5]sin^+'YRCOS^]]+((a-|)7L)[^I;T+^R  ■ j])eog 
C7'?[r?sin0-(i -^WIeeg  ( 1 3) 

a'2  [ijsin0-(T-»i5j1^5cos^]EoG 


where 

0L=  [(r)2+(*Z)2+rZcos0]y2 
Pr=  [(rf+tfZ^-rZcos*]1'*, 
7l  =  (r//9i)sin0, 

7r=  (r/^infl, 
o  =  [r2  +  h2  -  2rhsintf]1'2, 
rj  =  (-r/a)cos0. 


<t>  = 


0°  for  vertical  EOG 
90°  for  horizontal  EOG 


0°  for  vertical  interrogation  pulse 
[  90°  for  horizontal  interrogation  pulse 


a  =  m/n  =  (obtained  in  real  time  as  EOGvertiCa?pe7  EOGverticla|weI). 


We  measure  z  and  h,  as  well  as  ^and  /%for  each  electrode  (EEG  and  EOG).  We  then 
calculate  r  and  6  for  each  electrode,  and  then  calculate  yy  7R  a,  and  tj  for  each  electrode.  Finally, 
we  calculate  D’s  for  each  EEG/EOG  electrode  combination.  This  permits  correction  of  the  observed 
EEG  in  accordance  with  Equation  12.  The  next  figure  shows  an  example  of  a  correction.  The  cross- 
correlation  between  the  observed  EEG  and  the  vertical  EOG  was  0.75;  the  cross-correlation  between 
the  corrected  EEG  and  the  vertical  EOG  was  0.017. 
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VEOG 


Corrected 

EEG 

- 1 - 1 - 1 - 1 - 1 

TIME  (secs) 


Implementation  of  the  Technique  in  Hardware.  The  implementation  of  the  technique 
described  above  required  an  extremely  fast  microprocessor.  The  specification  that  the  portable, 
light-weight  device  must  fit  in  a  flight  suit  pocket  required  a  low-power  CMOS  microcontroller. 
These  specifications  resulted  in  the  selection  of  the  Intel  80C196  microcontroller.  This  chip  contains 
a  very  fast  microprocessor,  an  on-board  analog-to-digital  converter,  extremely  low  power 
consumption,  and  an  already  written  and  tested  Fast  Fourier  algorithm. 

There  were  many  obstacles  encountered  with  the  use  of  this  microcontroller.  There  is  a 
design  flaw  in  the  chip.  Intel  has  since  published  this  flaw  and  has  an  updated  chip.  The  flaw  is 
in  the  unsigned  divide  instruction.  The  result  from  this  instruction  is  either  the  correct  answer  or 
one  least  significant  bit  away  from  the  correct  answer.  This  doesn’t  seem  like  a  major  problem  on 
the  surface,  however  in  a  thirty  two  bit  divide  algorithm,  the  unsigned  divide  is  used.  What 
intermittently  occurs  is  an  incorrect  answer  which  is  off  by  one  least  significant  bit  in  the  HIGH 
word;  the  result  is  that  the  numerical  answer  is  off  by  65,536! 
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Another  flaw  in  the  chip  is  that  the  on-board  eighty  bytes  of  RAM  is  sporadically 
overwritten.  If  variables  located  in  the  onboard  RAM  are  forced  into  the  external  RAM  space,  the 
problem  seems  to  disappear. 

There  is  a  flaw  in  the  C  compiler  written  for  the  80C196.  A  locally  defined  variable  is  being 
overwritten  by  a  subroutine  containing  the  same,  but  locally  defined,  variable.  Using  identically 
defined  but  locally  defined  variables  is  standard  and  "legal"  in  C,  yet  this  compiler  does  not  seem 
to  properly  handle  this  situation. 

There  is  another  flaw  in  Intel’s  system.  The  in-circuit  emulation  system,  used  to  develop 
software  for  the  C196,  defines  the  ROM  as  zero  wait  state  memory.  This  causes  major  timing 
problems,  because  the  ROM  should  be  activated  with  the  user  programmable  wait  states,  which  can 
be  either  one,  two,  or  three.  Yet,  the  emulator  disregards  this  programmable  wait  state  number  and 
accesses  the  ROM  in  zero  wait  states.  Intel  has  been  notified  of  this  timing  fhw  This 
undocumented  discreptancy  makes  the  software  created  on  the  emulator  incompatable  with  the  Intel 
target  hardware  that  would  be  used  in  a  portable  device. 

There  also  have  been  other  general  problems  plaguing  this  effort.  This  microprocessor  is  an 
integer  based  machine.  This  leads  to  two  hurdles.  First,  the  resolution  of  the  mathematics  is 
truncated  to  digital  steps  and  not  continuous  functions.  The  ratio  of  one  to  one  half  is  two  in  the 
continuous  world,  yet  the  ratio  of  one  to  zero  (one  half  is  truncated  to  zero)  is  infinity  in  the  digital 
world.  The  second  hurdle  is  that  the  integer  set  is  bounded  at  -32,768  to  32,767.  This  constraint 
causes  the  programmer  to  scale  numbers  down  as  they  grow  close  to  the  bounds.  This  is  a  double- 
edged  sword,  since  the  function  of  scaling  is  division,  which  leads  to  truncation! 

This  constraint  of  integer  math  caused  us  to  require  the  use  of  a  host  PC  and  to  perform  the 
calculations  of  the  model  parameters  and  the  transfer  function  on  the  host  PC  (since  floating  point 
arithmetic  is  necessary  here).  This  eliminates  the  option  of  continually  interrogating  the  medium 
while  the  subject  is  ambulatory. 

The  blink  component  of  the  model,  although  an  excellent  advance  in  the  biophysical  model, 
added  complexity  to  the  correction  technique.  This  complexity  added  a  significant  amount  of 
computation  time  to  the  microprocessor  based  program.  This  result  was  that  there  was  only  time  for 
one  channel  to  be  corrected  with  the  full  model. 

In  order  to  achieve  this  correction  of  one  channel  in  the  allotted  (real)  time,  there  were 
several  "shortcuts"  that  were  necessary.  The  vertical  and  horizontal  transfer  functions  were  reduced 
from  an  array  of  complex  numbers  (one  for  each  frequency)  to  one  complex  number.  We  showed 
that  the  transfer  function  varied  less  than  ten  percent  over  the  frequency  spectrum.  This  allowed 
us  to  reduce  these  arrays,  yet  it  is  a  practical  variation  from  the  theoretical  ideal.  Another  shortcut 
was  the  elimination  of  the  square  root.  In  calculating  the  absolute  value  of  the  ratio  of  the  upper 
VEOG  to  the  lower  VEOG,  a  square  root  was  necessary.  We  showed  that  the  imaginary  component 
of  the  complex  ratio  was  very  close  to  zero,  so  we  took  the  real  component  of  the  ratio  instead  of 
the  absolute  value. 

There  were  many  technical  obstacles  that  complicated  this  research  study  and  prevented  us 
from  correcting  twelve  channels  of  EEG  for  EOG  artifact.  We  have  successfully  fabricated  a  device 
that  will  accurately  correct  one  channel  of  EEG  for  EOG  artifact.  As  indicated  in  the  previous 
discussion,  it  is  suceptable  to  sporadic  failures  caused  by  the  Intel  80CI96.  Application  specific 
circuits  and  chips  can  be  used  to  implement  this  correction  technique  on  multiple  channels,  however 
the  power  consumption  will  cause  the  battery  size  and  weight  to  increase  significantly.  This  would 
result  in  a  device  too  large  and  heavy  to  place  in  a  pocket  or  wear  on  the  body. 
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CONCLUSIONS 

From  a  scientific  point  of  view,  this  project  was  a  great  success  in  that  the  mathematical 
technique  was  extended  to  handle  blink  artifacts  in  a  non-arbitrary  biophysically  based  manner. 
From  an  engineering  point  of  view,  the  project  was  not  a  great  success  in  that  technological 
limitations  (computing  speed  of  CMOS  processors)  prevented  the  microprocessor  from  correcting 
more  than  one  EEG  channel  in  nearly  real-time. 

There  were  many  technical  obstacles  that  complicated  this  research  study  and  prevented  us 
from  correcting  twelve  channels  of  EEG  for  EOG  artifact.  We  have  successfully  fabricated  a  device 
that  will  correct  one  channel  of  EEG  for  EOG  artifact.  Application  specific  circuits  and  chips  could 
be  used  to  implement  this  correction  technique  on  multiple  channels,  however  the  power 
consumption  will  cause  the  battery  size  and  weight  to  increase  significantly.  This  would  result  in 
a  device  too  large  and  heavy  to  place  in  a  pocket  or  wear  on  the  body.  Full  implementation  of  a 
multichannel  man-borne  device  must  wait  advances  in  computer  hardware  technology. 
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APPENDIX 

The  following  appendix  is  the  operations  manual  for  the  EEG  Artifact  Rejection  System 
(EARS)  device. 
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EEG  Artifact  Rejection  System 


I.  PREFACE 


This  manual  contains  information  needed  for  the  operation  and  maintenance  of  an 
experimental  prototype  device,  which  provides  nearly  real-time  correction  of  artifacts  on  the 
steady  state  electroencephalogram  (EEG).  The  prototype  EEG  Artifact  Rejection  System 
(EARS)  is  a  battery-operated,  portable  device  that  is  designed  to  operate  in  a  variety  of 
experimental  operational  settings  (laboratory,  simulators,  and  aircraft)  and  is  intended  to  fit 
in  the  calf  pocket  of  a  flight  suit.  While  the  experimental  device  has  been  tested  in  a 
laboratory  setting,  it  has  NOT  been  evaluated  in  a  simulator  or  on  an  aircraft.  Furthermore, 
while  it  has  been  subjected  to  limited  human  testing,  it  is  not  an  approved  clinical  device  - 
it  is  an  experimental  prototype.  The  EARS  device  should  only  be  used  on  humans  under  the 
auspices  of  an  experimental  protocol  approved  by  a  duly  constituted  Internal  Review  Board. 


FIGURE  1.1:  ELECTROENCEPHALOGRAM  ARTIFACT  REJECTION  SYSTEM 
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II.  DESCRIPTION 


The  EARS  device  is  a  portable  unit  consisting  of  a  processor  unit  (which  fits  in  a  flight 
suit  calf  pocket)  and  an  electrode  input  selector  unit  (a  "relay  box",  which  is  intended  to  be 
slung  about  a  subject’s  neck,  as  shown  in  Figure  II.  1). 

The  relay  box  accepts  a  standard  multipin  EEG  electrode  input  connector  as  well  as  nine 
(9)  additional  electrode  connections  (one  pair  of  ear  electrodes,  one  pair  vertical  stimulation 
electrodes,  one  pair  horizontal  stimulation  electrodes,  one  pair  vertical  EOG  electrodes,  and 
one  single  horizontal  EOG  electrode).  These  nine  electrodes  and  any  one  of  the  24  EEG 
electrodes  (see  Figure  IV.  1)  are  transmitted  to  the  processing  unit  via  a  cable  using  standard 
25  pin  D  connectors. 

The  processing  unit  communicates  with  the  relay  box  via  its  25  pin  D  connector,  with  a 
host  PC  via  its  9  pin  D  connector,  and  with  an  analog  output  signal  recorder  via  its  9  pin 
round  connector.  The  operating  mode  of  the  processing  unit  can  be  selected  from  a  menu, 
displayed  when  the  device  is  connected  to  a  host  computer.  The  processing  unit  operates 
in  one  of  three  major  modes:  real-time  mode,  interrogation  mode,  and  correction  mode.  In 
the  "real-time"  mode,  the  unit  acts  as  a  conventional  biopotential  amplifier  system.  In  the 
"interrogation"  mode,  the  unit  acquires  EOG  data  and  EEG  data  (only  from  the  one  selected 
channel)  for  use  by  the  host  PC  to  compute  model  parameters  for  the  "correction"  mode.  In 
the  "correction"  mode,  the  unit  acquires  EEG  and  EOG  data,  computes  the  EOG  contribution 
to  the  EEG  data  (using  the  model  and  model  parameters),  subtracts  the  EOG  contribution 
from  the  EEG  signal,  and  outputs  the  corrected  signals  with  a  few  seconds  delay. 
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FIGURE  II.  1:  ELECTRODE  INPUT  SELECTOR  UNIT  (RELAY  BOX) 
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III.INSTALLATION 


The  hardware  consists  of  four  components.  There  is  the  EARS  main  processor  unit,  an 
electrode  input  selector  unit,  a  interconnection  cable  that  connects  these  two  units,  and  an 
analog  output  cable.  A  separate  RS-232  9-pin  "D"  communications  cable  must  be  provided 
for  the  communication  with  the  host  PC.  The  "relay  box"  has  two  screws  on  the  top,  that 
when  unscrewed,  allow  the  lid  to  open  and  the  batteries  to  be  replaced  as  shown  in  Figure 
III.  1 .  Figure  III.2  shows  the  parallel  battery  terminal  connectors  that  permits  replacing  the 
batteries  without  interrupting  operation  of  the  unit. 

Software  for  the  host  PC  is  contained  on  a  3.5"  disk.  One  can  install  this  software  by 
copying  the  disk  onto  a  hard  disk.  The  PC  must  be  an  IBM  PC/ AT/XT  with  a  numeric 
coprocessor.  This  is  done  by  typing  "copy  a:*.*"  when  in  the  desired  directory  on  the  hard 
disk.  The  user  is  now  ready  to  operate  the  EARS  system. 
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FIGURE  III.l:  POSITIONING  OF  THREE  9V  BATTERIES 


CMS  Engineering  Gorp. 


FIGURE  III.2:  BATTERY  CONNECTORS  FOR  UNINTERRUPTED  OPERATION 
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IV.  OPERATION 


The  operation  of  the  EEG  Artifact  Rejection  System  is  relatively  simple,  yet  there  is  a 
specific  protocol  that  must  be  followed  to  ensure  proper  function. 

Once  the  EEG  and  reference  electrodes  (Figure  VII.B.15)  and  the  nine  additional 
electrodes  (shown  in  Figure  II.  1)  are  attached  to  the  subject,  the  unit  can  be  powered  up. 
This  is  achieved  by  connecting  the  "relay  box"  to  the  main  processor  unit  using  the  cable 
provided.  Whe :  this  is  done,  the  power  indicator  on  the  "relay  box"  will  illuminate.  A 
flashing  light  indicates  a  low  battery.  The  RS-232  cable  must  now  be  connected  from  the 
processor  unit  to  a  IBM  PC/XT/AT  personal  computer  with  a  numeric  coprocessor.  A 
numeric  coprocessor  is  required  for  the  Fortran  software  to  operate. 

Type  "EARS"  and  then  a  carriage  return  to  enter  the  first  of  the  two  programs  (see 
Section  VII.C  for  software  description).  When  the  screen  goes  completely  blank 
(approximately  five  seconds  later),  push  the  carriage  return  once  again.  The  EARS  menu  will 
appear  on  the  screen. 

GMS  Engineering  Corporation 
EEG  Artifact  Rejection  System 

N  -  Channel  Nunber  Selection 
L  -  LED  Light  Level 
R  -  Real  Time  Monitoring 
P  -  Calibration  Pulses 
I  -  Interrogation 
C  -  Correction 

Enter  RESPONSE  > 

One  can  select  the  channel  that  the  processor  unit  will  correct  by  typing  "N".  The  system 
will  prompt  the  user  for  the  desired  channel  number  (1-24).  If  a  carriage  return  is  pushed 
without  entering  a  number,  the  current  channel  number  is  selected.  The  default  is  Chann 
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#1.  Figure  IV. 1  delineates  the  correspondence  between  the  channel  numbers  (1-24),  the  37- 
pin  "D"  connector  pins  and  the  normal  EEG  derivations  connected  to  those  pins. 


GMS  EEG 

37  PIN  D 

NORMAL  EEG 

CHANNEL  # 

CONNECTOR  PIN 

DERIVATION 

1 

1 

FP1 

2 

20 

FP2 

3 

2 

F3 

4 

21 

F4 

5 

3 

C3 

6 

22 

C4 

7 

4 

P3 

8 

23 

P4 

9 

5 

01 

10 

24 

02 

11 

6 

F7 

12 

25 

F8 

13 

7 

T3 

14 

26 

T4 

15 

8 

T5 

16 

27 

T6 

17 

28 

Cz 

18 

10 

Fz 

19 

29 

Pz 

20 

11 

Fpz 

21 

30 

Oz 

22 

14 

C3* 

23 

33 

Cz' 

24 

15 

C4' 

FIGURE  IV.  1:  CHANNEL  NUMBERING  SYSTEM 


One  can  change  the  power  indicator  light  level  that  appears  on  the  "relay  box"  by  typing 
"L".  This  will  permit  low  level  light  operations.  The  system  will  prompt  the  user  for  the 
desired  light  level  (1-255).  The  smaller  the  number  the  less  intense  is  the  light.  If  a  carriage 
return  is  pushed  without  entering  a  number,  the  current  light  level  is  selected.  The  default 
is  128. 


The  option  "R"  will  allow  the  user  to  monitor  the  three  EOG  channels  and  the  selected 
EEG  channel  from  the  outputs  on  the  EARS  unit.  To  exit  this  routine,  just  push  any  key 
on  the  PC  keyboard.  This  will  bring  the  user  back  to  the  main  menu. 
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The  option  "P"  will  allow  the  user  to  monitor  the  three  EOG  channels  and  the  selected 
EEG  channel  from  the  outputs  on  the  processor  unit.  A  train  of  calibration  pulses  will  ride 
on  the  outputs  for  approximately  five  minutes.  This  is  caused  by  application  of  a  single 
calibration  pulse  applied  to  all  the  input  channels.  The  amplitude  of  this  calibration  pulse 
is  1  mV.  Each  pulse  is  fifty  milliseconds  in  duration,  and  there  is  approximately  one  second 
between  pulses.  This  aids  the  user  in  adjusting  the  desired  gain  for  each  channel.  To  exit 
this  routine,  just  push  any  key  on  the  PC  keyboard.  This  will  bring  the  user  back  to  the 
main  menu. 

The  option  T  is  used  for  interrogating  the  medium  (subject).  This  routine  requires 
approximately  two  minutes.  The  direct  drive  signals  will  be  output  on  the  interrogation 
quadropole  (the  four  electrodes  on  the  forehead).  After  this  interrogation  process  is 
completed,  the  host  PC  screen  will  prompt  the  user  to  store  the  appropriate  data  for 
processing.  The  prompts  provided  on  the  PC  screen  are:  push  ’PgDn’,  then  type  "7",  and  then 
type  "drive”  and  carriage  return.  Then  push  the  uppercase  "A".  The  data  will  stream  across 
the  screen  and  into  a  file  on  the  disk. 

When  the  screen  prompts  the  user  to  exit  EARS,  push  ’ALT-X’  and  then  "Y".  The  user 
will  now  be  in  DOS.  The  second  program  should  be  run  by  typing  "EEG"  and  a  carriage 
return.  This  program  will  prompt  the  user  to  enter  the  geometrical  distances  (in  mm)  from 
the  eyes  to  the  EOG  and  selected  EEG  electrodes,  as  well  as  the  distance  between  the  eyes 
and  the  distance  from  the  center  of  the  eyes  to  the  quadropole.  These  should  be  carefully 
measured  using  a  soft  cloth  tape  measure.  When  these  parameters  are  entered,  the  program 
then  calculates  and  fine  tunes  the  model  coefficients  and  the  medium  transfer  function. 

When  this  program  is  finished,  the  screen  will  prompt  the  user  to  run  the  EARS  program 
once  again.  One  does  this  by  following  the  same  instructions  as  above.  When  the  main  menu 
appears,  choose  the  "C"  option  to  begin  correcting  the  selected  EEG  channel.  The  user  will 
be  instructed  to  push  ’PgUp’,  ”7",  and  type  "correct"  and  a  carriage  return.  The  appropriate 
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model  parameter  data  will  be  transferred  to  the  main  EARS  processor  unit,  and  the  EEG 
correction  will  begin. 

The  RS-232  cable  can  now  be  disconnected.  The  analog  outputs  are  as  follows. 


Output  #1 

Real  Time  EEG 

Output  #2 

Real  Time  Event  Trigger 

Output  #3 

Delayed,  Corrected  EEG 

Output  #4 

Delayed  Event  Trigger 

Output  #5 

Delayed,  Uncorrected  EEG 

Output  #6 

Delayed  Horizontal  EOG 

Output  #7 

Delayed  Vertical  Upper  EOG 

Output  #8 

Delayed  Vertical  Lower  EOG 

The  blinking  of  the  power  light  on  the  "relay  box"  means  that  the  batteries  are  getting 
low,  and  must  be  changed  within  the  hour.  The  batteries  can  be  changed  WITHOUT 
interrupting  operation  by  putting  three  new  batteries  on  the  reverse  side  of  the  battery  clip, 
and  then  taking  out  the  three  old  batteries.  The  power  light  should  then  be  continuously  on. 
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V.  STORAGE 


Turn  off  the  battery  power  to  the  device  by  disconnecting  the  cable  between  the 
processor  unit  and  the  relay  box.  The  LED  indicator  will  extinguish. 

Disconnect  all  electrodes  and  cables  from  the  device. 

Wipe  off  any  debris  from  the  external  surfaces  of  the  EARS  unit  before  storage.  A  soft 
cloth  dampened  with  water  or  a  mild  soap  and  water  solution  can  be  used.  Do  not  apply 
organic  solvents  to  this  prototype  unit. 

To  conserve  battery  life,  remove  the  three  9V  batteries  from  the  unit.  Do  not  leave  the 
batteries  in  the  unit,  if  long  term  storage  is  intended. 

Return  the  unit  to  its  original  transport  container  or  another  equivalent  storage/protection 
container. 
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VI.  THEORY  OF  OPERATION 


The  EARS  device  is  based  on  the  idea  that  eye  movements  and  blinks  contribute  to  the 
observed  EEG  signals.  If  the  electrical  signal  characteristics  of  these  eye  movements  and 
blinks  are  known  and  the  medium  through  which  these  signals  propagate  to  the  EEG 
observation  sites  (EEG  electrode  sites)  is  characterized,  then  this  unwanted  influence  can  be 
mathematically  removed.  The  removal  of  this  influence  is  the  correction  process.  The 
mathematical  model  describing  this  correction  process  which  is  implemented  in  the  software 
in  the  EARS  device  is  described  in  this  section. 

The  Biophysical  Model.  The  transection  of  the  face  across  the  three  dimensional  ocular 
dipole  field  (caused  by  the  corneo-retinal  potential  in  the  eye)  yields  a  surface  image  dipole 
propagating  on  the  scalp.  This  surface  image  dipole  can  be  modelled  to  incorporate  both  the 
eye  movement  and  the  eye  blink.  The  eye  movement  produces  a  symmetric  dipole,  while  the 
blink  produces  an  asymmetric  dipole. 

General  Dipole  Representation.  A  dipole  source,  symmetric  or  asymmetric,  is  the 
superposition  of  two  point  sources  separated  by  a  distance.  The  point  source’s  electric  field 
propagates  as  a  function  of  1/r2.  The  voltage  at  any  point  is  described  by  V=kq/r,  where  k 
is  Boltzmann’s  constant,  and  q  is  the  amount  of  charge.  The  surface  image  dipole  is 
described  here. 

Figure  VI.  1  shows  two  point  sources  separated  by  a  distance  (L).  The  voltage  (V) 
appearing  at  point  Q  is  derived  as  follows. 

V  =  kq{[m/(r-iLsin0)]+[-n/(r+{Lsin0)]}.  (1) 
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Rearranging  Equation  1  yields 

V  =  kq{[(m-n)r+(m+nHLsin0]/(r2-±L2sin20)}.  (2) 

Since  r»L,  we  can  simplify  Equation  2: 

V  =  n[(kq/r2)iLsin0(a+l )  +  (kq/r)(a-l)],  (3) 

where  a  =  m/n. 

FIGURE  VI.  1:  TWO  As  a  note,  if  a=l  (eye  movement)  and  the  dipole  is  symmetric, 
POINT  SOURCES 

SEPARATED  BY  Equation  3  reduces  to, 

DISTANCE  L. 

V  *  (nkq/r2)Lsinfl  ■  Ar  2sin0.  (4) 

where  A  =  nrqL. 

Furthermore,  it  is  important  to  note  that  the  zero-potential  line  of  the  dipole  is  the  x- 
axis  when  a=l  (sin0=O).  When  a*l,  the  zero-potential  line  becomes  a  circle  described  by, 

x2  +  (y+G)2  =  G2,  (5) 

where  G=iL[(a+l )/(«-!)]. 
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Selection  of  a  F acial  Coordinate  System.  In  order  to  spatially  represent  the  ocular  dipoles 
or  the  stimulus  dipole  in  planar  geometry,  we  must  select  a  coordinate  system.  This  is  shown 
in  Figure  VI.2.  The  point,  Q,  in  this  figure  represents  an  electrode.  The  subscript  L  is  used 
to  show  reference  to  the  left  eye,  the  subscript  R  is  used  to  reference  the  right  eye.  The 
electrode  is  a  distance  from  the  left  eye  and  /?„  from  the  right  eye.  Reference  to  the 
stimulation  dipole  is  indicated  by  S.  Our  facial  coordinate  system  has  its  origin  at  the 
geometric  center  of  the  two  eyes  (which  are  separated  by  a  distance  z).  The  stimulus  dipole 
is  located  a  distance  h  above  the  origin.  The  electrode  is  a  distance  r  from  the  origin  of  the 
coordinate  system.  The  stimulation  dipole  lies  on  the  y-axis  as  the  eyes  lie  on  the  x-axis. 


FIGURE  VI.2:  COORDINATE  TRANSFORMATION  APPLIED  TO  POINT 

Q- 
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Spatial  Resolution  of  the  Bi-ocular  Dipoles.  Using  our  coordinate  system  and  dipole 
representation  described  above,  we  will  model  the  two  ocular  dipoles  as  one  complex 
mathematical  function  based  at  the  origin  of  the  coordinate  system. 


From  the  law  of  cosines: 


From  the  law  of  sines: 


fit  =  [r2  +  iZ2  +  rZcosfl]* 

(6) 

=  [r2  +  iZ2  -  rZcos*]* 

(7) 

7l  =  ~<t>  +  sin  '[(r//9L)sin0] 

(8) 

7r  =  tt-(f>  -  sin  ’[(r/£R)sin0] 

(9) 

where:  <f>  =  angle  of  the  ocular  dipoles  (zero-potential  line). 

We  now  substitute  these  identities  (Eqns.  6-9)  into  the  dipole  equation  defined  earlier  (Eqn. 
3)  and  sum  the  two  ocular  dipoles  to  result  in  one  function.  We  obtain  this,  in  a  general 
form,  for  any  electrode: 

V  =  i(a+l  )A{[/?L2[-(l-7L2),/2]sin$+7Lcos$]+[/9R2[(l  -7R2)!/z]sin0+7Rcos$]}+(A/L)(a- 1  )iPi  '+Pr  '), 

(10) 


where: 

0l  =  [( r  )2+(  i  Z  )2+r  Zcosfl] ,/2 
P*  =  [(r)2+(iZ)2-rZcos0]1/2 
7l  =  (r/£L)sin0 
7r  =  (r//3„)sin0 
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Spatial  Representation  of  the  Stimulus  Dipoles.  We  can  similarly  describe  the 
representation  of  the  stimulus  dipole  in  our  new  coordinate  system  as  we  have  described  the 
ocular  dipoles  above.  Using  Figure  VI.2  and  Equation  4: 

V  =  Btcr^fjsin^-Cl-fj^’^cosV']}  (11) 

where: 

o  =  [r2  +  h2  -  2rhsin0],/2 
T)  =  (-r/a)cosfl 

rj)  =  angle  of  the  stimulus  dipole  (an  analog  of  <f> ) 

Equations  10  &  11  and  the  associated  identities  are  the  basis  of  the  mathematical  model 
which  will  be  used  in  the  calculation  of  the  transfer  function. 

Mathematical  Relationship  Between  Resultant  Ocular  Dipoles  &  Stimulus  Dipoles.  We  can 
now  use  our  basic  equations  and  our  coordinate  system  to  correct  the  putative  transfer 
function  measured  by  surface  dipoles.  The  EOG  artifact  correction  equation  in  the 
frequency  domain  is: 

TOrTV EEG  =  065 V EEG  -  EOC  [Seeg/SeOGJ  D  (12) 

where  D  is  the  geometrical  correction  factor  between  the  stimulus  dipole  and  the  ocular 
dipoles.  S  denotes  the  surface  dipole  stimulation  response  and  ob'V  denotes  the  naturally 
occurring  response.  The  subscript  "EEG"  and  "EOG"  refer  to  the  electrode  recording  the 
response;  the  superscripts  "corr"  and  "obs"  refer  to  the  corrected  and  observed  potential, 
respectively. 
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We  can  obtain  D  by  manipulation  of  the  equations  described  above. 


{*(<*+ 1  -'TL2)v2]sin^+7Lcos^]+I^„2[(l  -'yR2),/2]sin^+'ifRcos^]]+((a- 1  )/L)[j8L \p^])aa 

CT^sin^-(f-jj2)1/2cos^]EEG  (13) 

oY»?sin^-(l  -»?2),/2cosV»]eog 


where. 

Pi  =  [(r)2+(iZ)2+rZcos(?]1/a, 
Pr  =  t(r)2+(  i  Z)2-  rZcosfl] 1/2, 
TTl  =  (r/POsind, 

1*  =  (r/£„)sin0, 
a  =  [r2  +  h2  -  2rhsin0]1/2, 
q  =  (-r/cr)cos0, 


^  = 


0°  for  vertical  EOG 
[  90“  for  horizontal  EOG 


V>  = 


0“  for  vertical  interrogation  pulse 
[  90“  for  horizontal  interrogation  pulse 


a  =  m/n  =  (obtained  in  real  time  as  EOGvrr,icjI,',pp"  /  EOGTertcaitew"  ). 


We  measure  z  and  h,  as  well  as  pL  and  /?R  for  each  electrode  (EEG  and  EOG).  We  then 
calculate  r  and  6  for  each  electrode,  and  then  calculate  7L,  7r,  and  r\  for  each  electrode. 
Finally,  we  calculate  D’s  for  each  EEG/EOG  electrode  combination.  This  permits  correction 
of  the  observed  EEG  in  accordance  with  Equation  12.  Figure  VI.3  shows  an  example  of  a 
correction.  The  cross-correlation  between  the  observed  EEG  and  the  vertical  EOG  was  0.75; 
the  cross-correlation  between  the  corrected  EEG  and  the  vertical  EOG  was  0.017. 
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FIGURE  VI.3:  OBSERVED  AND  CORRECTED  EEG. 


I 
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VII. 


TROUBLESHOOTING  GUIDE 


A.  GENERAL 

The  general  troubleshooting  protocol  of  the  EARS  system  is  extremely  simple. 
If  the  power  light  is  flashing,  the  batteries  need  to  be  changed.  If  the  output  levels 
become  "flat"  (no  signal),  the  batteries  need  to  be  changed.  If  the  batteries  are  new,  then 
the  system  must  be  serviced  by  authorized  personnel. 


B.  HARDWARE  DESCRIPTION 

The  hardware  schematics  are  shown  in  Figures  VII.B.2  through  VII.B.15.  Figure 
VII.B.l  is  a  block  diagram  of  the  complete  system.  Each  analog  channel  contains  four 
user-adjustable  potentiometers.  One  controls  the  gain  of  the  channel;  and  one  controls 
the  offset  of  the  channel.  The  other  two  are  for  fine  tuning  the  60  Hz  notch  filter. 
These  are  factory  calibrated,  and  should  hold  their  calibration  for  several  months  or 
longer. 


The  calibration  circuit  is  a  floating  voltage  source  that  is  switched  into  series  with 
the  inverting  input  of  the  instrumentation  amplifier.  This  level  can  be  changed  by  the 
user  with  a  potentiometer.  See  Section  VIII  for  further  details. 

The  digital  circuitry  consists  of  a  microcontroller  and  memory.  There  are  digital- 
to-analog  output  converters/amplifiers  which  allow  the  user  to  monitor  the  EEG  and 
EOG  channels. 
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FIGURE  VII.B.l:  HARDWARE  BLOCK  DIAGRAM 
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FIGURE  VII.B.2:  ISOSWITCH  BOARD  SCHEMATIC 
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FIGURE  VII.B.3:  ISOSWITCH  BOARD 
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FIGURE  VII.B.4:  INSTRUMENTATION  AMP  &  HORIZONTAL  STIMULUS  SCHEMATIC 
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FIGURE  VII.B.5:  INSTRUMENTATION  AMP  &  HORIZONTAL  STIMULUS  BOARD 
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FIGURE  VII.B.6:  ANALOG  AMPLIFIERS  SCHEMATIC 
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FIGURE  VII.B.7:  AUTOCALIBRATION  &  VERTICAL  AUTOSTIMULATION 

SCHEMATIC 
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FIGURE  VII.B.8:  ANALOG  AMPLIFIER/ A UTOCA L/A UTOSTIM  BOARD  SILKSCREEN 


GMS  Engineering  Carp. 


38 


EEG  Artifact  Rejection  System 


FIGURE  VII.B.9:  DIGITAL  CIRCUIT  SCHEMATIC 
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FIGURE  VII.B.10:  ANALOG  INPUT  SIGNAL  LIMITER/BUFFER  SCHEMATIC 
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FIGURE  VII.B.l  1:  DIGITAL  BOARD  SILKSCREEN 
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GHS  ENGINEERING  CORP. 

DIGITRL  BOARD _ 

AUTOMATED  DESIGN  SERVICES 
SILKSCREEN  i  T 


FIGURE  VII.B.12:  RELAY  BOX  AUXILIARY  CIRCUIT  SCHEMATICS 
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s«4Kirjy" 


FIGURE  VII. B. 13:  RELAY  BOARDS 
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FIGURE  VII.B.14:  RELAY  BOX  CONNECTOR  SCHEMATIC 
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FIGURE  VII. B. 15:  37  PIN  ”D"  CONNECTOR  PINOUT  DESIGNATION 
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C3'  =  ACT  1  on  Spectrum  32  Headbox. 

Cz'  =  ACT  2  on  Spectrum  32  Headbox. 

C4'  =  ACT  3  on  Spectrum  32  Headbox. 

Gnd  =  Iso  Gnd  on  Spectrum  32  Headbox. 

NOTE:  In  order  to  use  the  Prime  electrodes 
(  C3',  Cz',  or  C4' )  you  must  specify  the  correct 
active  input  (ACT  1,  2,  or  3)  on  the  Spectrum  32. 


c. 


SOFTWARE  DESCRIPTION 


There  are  two  components  of  the  software  for  the  EEG  Artifact  Rejection  System.  These 
components  are  the  EPROM  code  which  is  contained  in  the  unit  and  the  PC  based  code 
which  is  contained  on  the  disk.  The  complete  software  are  appended  here. 

The  EPROM  code  ("EEG"  written  in  C,  and  "FFT_FOR"  written  in  assembler)  consists 
of  six  parts.  The  first  part  is  the  serial  communications  routine.  This  allows  the  user  to 
control  what  the  unit  does  as  well  as  selection  of  the  desired  channel  and  the  LED  level.  This 
also  is  the  vehicle  for  data  transfer.  The  second  part  is  the  service  of  the  hardware  and 
housekeeping  routines.  This  part  is  necessary  for  all  software.  The  processor,  memory,  and 
peripheral  circuits  require  specific  signals  and  protocols.  This  is  accomplished  in  the 
housekeeping  and  hardware  service  routines.  The  batteries  are  checked  here  as  well.  The 
third  and  fourth  parts  are  similar  to  each  other.  They  perform  the  real-time  monitoring 
function.  The  fourth  part  adds  the  calibration  pulses  if  the  user  desires  them.  The  fifth  part 
performs  the  direct  interrogation.  The  vertical,  horizontal,  and  calibration  drives  are  output, 
and  the  electrode  signals  are  recorded  and  stored  for  later  processing.  The  last  part  is  the 
application  of  the  model  for  correction  of  the  EEG. 

The  PC-based  code  ("EARS”  and  "EEG"  which  consists  of  "SHELL,  "DCALC2",  "SGM", 
and  "FTUNE2"  all  written  in  Fortran)  consists  of  four  parts.  The  first  part  is  the  serial 
communications  routine.  This  is  the  PC  side  of  the  communications  described  above.  This 
program  is  "EARS".  The  second  (DCALC2),  third  (SGM),  and  fourth  (FTUNE2)  parts  are 
all  contained  in  "EEG".  The  second  part  allows  the  user  to  enter  the  geometrical 
measurements,  and  calculates  the  model  parameters  from  these  measurements.  The  third  part 
takes  the  interrogation  data  and  calculates  the  medium  transfer  function  for  the  specific  test 
subject.  The  last  part  statistically  fine  tunes  the  model  coefficients  obtained  in  the  second 
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part  using  the  transfer  function  obtained  in  the  third  part.  The  coefficients  and  the  transfer 


function  must  be  then  sent  back  to  the  main  processor  unit  before  the  correction  can 


commence.  This  is  accomplished  by  option  "C"  in  "EARS". 


FIGURE  VII.C.l:  SOFTWARE  BLOCK  DIAGRAM. 
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VIII. 


CALIBRATION  PROCEDURES 


By  choosing  the  "P"  (Calibration  Pulses)  option  on  the  main  menu,  the  user  can  tune  the 
gains  and  offsets  while  viewing  the  input  signal  with  a  common  calibration  pulse  riding  on 
all  the  channels.  This  amplitude  of  this  calibration  pulse  can  be  adjusted  by  the  user.  This 
may  be  required  when  the  channel  gains  are  changed  to  the  upper  or  lower  extremes.  This 
procedure  is  as  follows. 

The  EARS  main  processor  unit  must  be  opened  by  removing  the  seven  screws  on  the  face 
opposite  the  potentiometers.  Once  the  top  is  removed,  the  several  printed  circuit  boards  in 
the  unit  will  become  visible.  On  the  side  of  the  main  printed  circuit  board  closest  to  the  RS- 
232  connector  is  a  jumper  and  two  monitoring  pins.  The  jumper  can  be  moved  to  the  next 
position,  which  bypasses  the  relay  and  continuously  applies  the  calibration  battery  to  the 
circuit.  The  monitoring  pins  can  be  used  to  measure  the  exact  voltage  of  the  calibration 
pulse.  A  microvolt  meter  must  be  used  for  this  purpose.  The  potentiometer  control  on  the 
outside  opposite  face  of  the  unit  can  be  used  to  adjust  the  voltage.  The  jumper  MUST  be 
placed  back  in  the  factory  position  after  this  is  complete  and  before  the  unit  is  closed. 
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DOS  3.30  (038-M)  C96  COMPILER  VI. 1,  COMPILATION  OF  MODULE  EEG 
OBJECT  MOOULE  PLACED  IN  EEG.obj 

COMPILER  INVOKED  BY:  D:\196\C96\C96.EXE  EEG.C96  MOOEL(196)  OPTIMIZE(2)  REGISTERS(IOO) 
stmt  level  incl 


EEG.C96  -  EOG  Artifact  Removal  System 
Implemented  in  Intel  iC96  for  the  80C196KA 

Steven  M.  Falk 
Jeffrey  C.  Sigl 

Created:  August  22,  1989 

Version  No. :  1 .0 

Last  Update:  February  12,  1990 

GMS  Engineering  Corporation 
8940 -D  Route  108 
Columbia,  MD  21045 
(301)  995-0508 


0000 

0100 

1FFE 

2000 

2014 

2018 

2019 

2020 
2030 
2040 
2080 

6000 


MEMORY  MAP 
SP.SFRS, REGISTERS 


EXTERNAL  MEMORY 


I/O  PORTS  3  &  4 
LOW  INTRPT  VECTORS 


RESERVED  (INTEL) 
CCR 


RESERVED  (INTEL) 
EPROM  SECURITY  KEY 


HIGH  INTRPT  VECTORS 
RESERVED  (INTEL) 


CODE 


I  / 
I  / 
■1/ 

|\ 

I  \ 

I  \ 


8K  RAM 


16K  ROM 
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EOOO 


FFFF 


| 32K  RAM 

I 


Port  pin  assignments 


Pin 

Direction 

Function 

P1.x=1 

P1.x=0 

P1.0 

output 

Cal 

OH 

OFF 

P1.1 

output 

Stiml 

ON 

OFF 

PI. 2 

output 

Stim2 

ON 

OFF 

PI. 3 

output 

Relay  Address 

(LSB) 

Pi. 4 

output 

Relay  Address 

Pi. 5 

output 

Relay  Address 

Pi. 6 

output 

Relay  Address 

PI. 7 

output 

Relay  Address 

(MSB) 

P2.5 

output 

LEO 

P2.2 

input 

Event  Trigger 

In 

Input 

Channels 

#  Stored 

in  data_buffer[#]  1)  l) 

Output  Chan 

0 

EEG 

EEG  uEEG(t) 

1 

HEOG 

HEOG 

ETO(t) 

2 

VEOGu 

VEOGu 

cEEG(t-T) 

3 

VEOGl 

VEOGl 

ETO(t-T) 

4 

5 

6 

Event  Trigger  Out 

HEOG(t-T) 

VEOGu(t-T) 

VEOGl(t-T) 

This  routine  must  be  linked  with: 


e_int.obj  (macro  El) 
CCR.ABS  (ccb) 

user. lib  (patched  c96.lib) 
fft_for.obj  (the  FFT) 
cstart.obj  (main  module) 
plm96. I ib 


i.e., 

eeg.obj,cstart.obj,ccr.abs,e_int.obj,fft_for.obj,  & 
user. I ib,plm96. I ib  to  eeg  ixref  & 
ro(2000H-2013H,2018H-2018H,2030H-203FH,2080H-5FFFH)  & 
ra(1 AH- 1FFDH( STACK), 6000H- OOF FFH) 
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1 

1  rtine_flag  =0  =  real-time 

1  rtine_flag  =1  =  interrogation 

1  rtine_flag  =2  =  collection  of  epochs  for  fine  tuning 

1  rtine_flag  =3  =  correction 

1  rtine_flag  =7  =  real-time  with  cal  pulses 

1  rtine_flag  =99  =  stop 

1 

1  V 

1 
1 
1 

1  /*  Headers  V 

1 

1  (((include  <800196. h>  /*  80C196  I/O  registers  */ 

34  ^include  <ctype.h> 

35  ^include  <setjmp.h> 

36  #include  <stdio.h> 

44  ^include  <stdlib.h> 

52  ^include  <string.h> 

65 

65 

65 

65 


65 

/*  Definitions 

*/ 

65 

65 

#def i ne 

PRESCALE 

64 

65 

65 

#def i ne 

BUFF  START 

64 

65 

#def ine 

BUFF  DIFF 

384 

/*  0.75*SAMP  NUM  */ 

65 

#def ine 

CHAN  BASE 

0x08 

/*  ADC  Channel  base  (ADC0)  */ 

65 

#def i ne 

CHAN  HUM 

4 

/*  number  of  input  (ADC)  channels  */ 

65 

#def ine 

CHAN  HI 

3 

/*  CHAN  NUM- 1  */ 

65 

#def i ne 

CHAN  PI 

5 

/*  CHANNUM+1  */ 

65 

#def ine 

DELTA  T 

0x0E371 

65 

#def ine 

LEFT_LIM 

128 

/*  0.25*SAMP  NUM  */ 

65 

#def ine 

RIGHT  LIM 

383 

/*  (0.75*SAMP  NUMJ-1  */ 

65 

#def ine 

SAMP  NUM 

512 

/*  number  of  samples  */ 

65 

#def ine 

TRO  PULSE 

4 

/*  30  msec  */ 

65 

#def ine 

CAL  NUM 

60 

/*  cal/stim  samples  in  buffers  */ 

65 

#def i ne 

CAL  NUM  Ml 

59 

/*  CAL_NUM- 1  */ 

65 

#def ine 

CAL  ON 

39 

/*  Stim  on  */ 

65 

#def i ne 

CAL  OFF 

19 

/*  Stim  off  */ 

65 

#def ine 

CAL  REP 

60 

/*  cal/stim  repetitions  for  avg  */ 

65 

#def ine 

DELTA  T  FAST 

OxOF8DE 

65 

#def ine 

SAMP  LIM 

200 

/*  sample  time  out  limit  */ 

65 

#def ine 

SAMP  X2 

1024 

65 

#def ine 

LOW  BATT 

650 

65 

#def ine 

E  LOW  BATT 

600 

65 

#def ine 

H 

0 

65 

#def ine 

V 

1 

65 

#def ine 

R 

0 

65 

#def ine 

I 

1 

65 

65  /*  Interrupt  service  function  assignments  */ 

65 

65  (fpragma  interrupt  (nmi_int=31 )  /*  NMI  interrupt  */ 
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65  #pragma  interrupt  (extint=29)  /*  EXTINT  Pin  interrupt  */ 

65  #pragma  interrupt  (receive=25)  /*  Serial  Port  Receive  interrupt  */ 

65  (Kpragma  interrupt  (samp=1)  /*  A2D  CONVERSION  COMPLETE  interrupt  */ 

65  IHpragma  interrupt  (time1=0)  /*  TIMER1  OVERFLOW  interrupt  */ 

65 

65 

65  /*  Function  declarations  */ 

65 

65  void  main(void); 

66 

66  void  nmi_int(void);  /*  Interrupt  Service  Routines  */ 

67  void  extint(void); 

68  void  receive(void); 

69  void  timel(void); 

70  void  samp(void); 

71 

71  void  seriat(void); 

72  void  senddatal int); 

73  int  recvdatal void); 

74  void  dac( int, int); 

75  void  err(void); 

76  void  fft_for(void); 

77 
77 
77 

77  /*  External  functions  */ 

77 

77  extern  void  enabint(void); 

78 
78 

78  /*  SFP  Images  */ 

78 

78  unsigned  char  im  iocO,  im_ioc1,  im_ioc?  im_sp_stat,  importl,  pmvim; 

79  register  unsigned  char  status  temp;  /*  defined  in  USER. LIB  */ 

80 
80 
80 

80  /*  Global  variables  */ 

80 


80  register  long  int  tltempl,  tltemp2,  tltemp3,  tltemp4,  tdenom; 

81  register  int  gainmag,  newpt,  alphal,  alpha2,  sgv_mag,  sgh_mag,  tunept,  templreg, 

temp2_reg; 

82  register  char  chan,  buff  num,  buff_not,  chan_end; 

83 

83  register  int  dv[8J; 

84 

84  int  calcounter,  calmode,  repcount,  cal  moni tor,  calcnt,  correctcnt; 

85  char  timerlflag,  ri_flag,  correctf lag,  error,  loop_ftag; 

86  char  restartf lag,  tro_counter,  rtineflag; 

87  int  battlevel,  battcnt,  itemp; 

88  int  xreal[512],  ximag[512); 

89  long  int  lltempl,  Utemp2,  lltemp3,  lltemp4,  denom; 

90  int  outbuf fer (2) ISAMPNUM] ; 

91  int  alpfia [SAMP  NUM]  ; 

92  int  eogf  12)12]  ISAMP  NUMJ  ; 

93  int  databuf fer (CHANP1] 12] ISAMPNUMJ ,  bat t_vol t [3J ; 

94  int  datatuner CCHANNUM] [SAMP_X2] ; 

95  int  gssu"1  [CAL  NUM] ; 
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96 

97 

98 

98 

98 

98 

98 

99 
100 
101 
102 
102 
102 
102 
102 
102 
102 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 
103 


int  gssum2 [CALNUM] ; 

int  gs  [3]  [2]  tCAL_NUM]  ; 

/*  Messages  V 

const  char  messltl  =  "  Type  ENTER  to  cont inue. . .\n\r"; 

const  char  cist]  =  {'\033‘ , '  [• ,  *2' , '  J‘>;  /*  Esc,t,2,J  */ 

const  char  bell  =  '\007'; 

const  char  stimportt]  =  {0x02,  0x04,  0x01};  /*  i/o  lines  for  stim/cal  */ 

/*  Square  Root  Table  */ 

const  int  sqrtt]  =  {  100,  121,  144,  169,  196,  225,  256,  289, 


324,  361,  400,  441,  484,  529,  576,  625,  676,  729,  784,  841,  900, 
961,  1024,  1089,  1156,  1225,  1296,  1369,  1444,  1521}; 


/*  Super  Gaussian  Window,  power  =  9  */ 

const  int  supgauU  =  {16,  20,  26,  32,  40,  49,  61,  74, 

90,  108,  129,  154,  182,  214,  251,  293, 

339,  391,  449,  514,  585,  662,  747,  840, 

940,  1048,  1165,  1289,  1423,  1564,  1714,  1872, 

2039,  2213,  2396,  2587,  2785,  2991,  3203,  3422, 

3647,  3879,  4115,  4357,  4603,  4853,  5106,  5363, 

5622,  5883,  6146,  6409,  6674,  6938,  7202,  7465, 

7728,  7988,  8246,  8503,  8756,  9006,  9253,  9497, 

9736,  9972,  10203,  10429,  10651,  10869,  11081,  11288, 

11490,  11687,  11878,  12065,  12245,  12421,  12591,  12756, 

12916,  13070,  13220,  13364,  13503,  13637,  13767,  13891, 

14011,  14126,  14237,  14343,  14445,  14543,  14636,  14726, 

14812,  14894,  14973,  15048,  15120,  15188,  15253,  15316, 

15375,  15431,  15485,  15536,  15585,  15631,  15675,  15717, 

15756,  15794,  15829,  15863,  15895,  15925,  15953,  15980, 

16006,  16030,  16052,  16074,  16094,  16113,  16131,  16148, 

16164,  16179,  16193,  16206,  16218,  16230,  16241,  16251, 

16260,  16269,  16278,  16285,  16293,  16300,  16306,  16312, 

16317,  16323,  16327,  16332,  16336,  16340,  16343,  16347, 

16350,  16353,  16355,  16358,  16360,  16362,  163e>4,  16366, 

16367,  16369,  16370,  16371,  16373,  16374,  16375,  16376, 

16376,  16377,  16378,  16378,  16379,  16380,  16380,  16380, 

16381,  16381,  16381,  16382,  16382,  16382,  16382,  16383, 

16383,  16383,  16383,  16383,  16383,  16383,  16383,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  *>384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 

16384,  16384,  16384,  16384,  16384,  16384,  16384,  16384, 
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103 

16384, 

16384, 

16384,  16384, 

16384, 

16384, 

16384,  16384, 

103 

16384, 

16384, 

16384,  16384, 

16384, 

16384, 

16384,  16384, 

103 

16384, 

16384, 

16384,  16384, 

16384, 

16384, 

16384,  16384, 

103 

16384, 

16384, 

16384,  16384, 

16384, 

16384, 

16384,  16384, 

103 

16384, 

16383, 

16383,  16383, 

16383, 

16383, 

16383,  16383, 

103 

16383, 

16382, 

16382,  16382, 

16382, 

16381, 

16381,  16381, 

103 

16380, 

16380, 

16380,  16379, 

16378, 

16378, 

16377,  16376, 

103 

16376, 

16375, 

16374,  16373, 

16371, 

16370, 

16369,  16367, 

103 

16366, 

16364, 

16362,  16360, 

16358, 

16355, 

16353,  16350, 

103 

16347, 

16343, 

16340,  16336, 

16332, 

16327, 

16323,  16317, 

103 

16312, 

16306, 

16300,  16293, 

16285, 

16278, 

16269,  16260, 

103 

16251, 

16241, 

16230,  16218, 

16206, 

16193, 

16179,  16164, 

103 

16148, 

16131, 

16113,  16094, 

16074, 

16052, 

16030,  16006, 

103 

15980, 

15953, 

15925,  15895, 

15863, 

15829, 

15794,  15756, 

103 

15717, 

15675, 

15631,  15585, 

15536, 

15485, 

15431,  15375, 

103 

15316, 

15253, 

15188,  15120, 

15048, 

14973, 

14894,  14812, 

103 

14726, 

14636, 

14543,  14445, 

14343, 

14237, 

14126,  14011, 

103 

13891, 

13767, 

13637,  13503, 

13364, 

13220, 

13070,  12916, 

103 

12756, 

12591, 

12421,  12245, 

12065, 

11878, 

11687,  11490, 

103 

11288, 

11081, 

10869,  10651, 

10429, 

10203, 

9972,  9736, 

103 

9497, 

9253,  9006,  8756,  8503 

.  8246, 

7988, 

7728, 

103 

7465, 

7202,  6938,  6674,  6409 

,  6146, 

5883, 

5622, 

103 

5363, 

5106,  4853,  4603,  4357 

.  4115, 

3879, 

3647, 

103 

3422, 

3203,  2991,  2785,  2587 

.  2396, 

2213, 

2039, 

103 

1872, 

1714,  1564,  1423,  1289 

.  1165, 

1048, 

940, 

103 

840,  747,  662, 

585,  514,  449 

.  391, 

339, 

103 

293,  251,  214, 

182,  154,  129 

.  108, 

90, 

103 

74,  61 

,  49,  40 

,  32,  26,  20. 

16); 

104 

104 

104 

104 

104 

void  main(void) 

104  { 

105  1  int  1,  j,  k; 

106  1 

106  1  /*  System  Configuration  */ 

106  1 

106  1  intmask  =  0x00;  /*  mask  all  interrupts  */ 

107  1  imaskl  =  0x00; 

108  1  enab_int();  /*  enable  globsl  interrupts  */ 

109  1 

109  1  im_ioc0  =  0x00;  /*  set  I/O  control  registers  */ 

110  1  iocO  =  im_ioc0; 

111  1 

111  1  im_ioc1  =  0x00; 

112  1  iocl  =  im_ioc1; 

113  1 

113  1  im_ioc2  =  0x00; 

114  1  ioc2  =  im_ioc2; 

115  1 

115  1  importl  =  0x00;  /*  initial  Portl  config.  */ 

116  1  ioportl  =  importl; 

117  1  pmwm  =  127; 

118  1  pwm_control  =  pmwm; 

119  1 
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158 

2 

for  (  j=0;  j  <8;  j++  )  < 

159 

3 

dac{j,i  temp- 128); 

160 

3 

for  (  k=0;  k<10;  k++  ) 

161 

3 

i 

161 

3 

> 

161 

2 

> 

161 

1 

161 

1 

for  (  i=0;  1  <  8;  i++) 

162 

1 

dac<  i,  64  ); 

163 

1 

163 

1 

163 

1 

163 

1 

/*  Go  to  main  menu  */ 

163 

1 

163 

1 

serialO; 

164 

1 

loop  flag  =  0; 

165 

1 

ri  flag  =  0; 

166 

1 

ipendl  8=  -0x02; 

167 

1 

167 

1 

167 

1 

/*  Initialize  Timerl  */ 

167 

1 

167 

1 

im  iocl  |=  0x04;  /*  enable  TIMER1  ovrflow  intrpt  */ 

168 

1 

iocl  =  im  iocl; 

169 

1 

169 

1 

timerl  flag  =  0;  /*  timerl  flag  */ 

170 

1 

loop  flag  =  0; 

171 

1 

wsr  =  OxOF;  /*  alternate  window  */ 

172 

1 

timerl  =  DELTA  T;  /*  load  7.813  msec  timer  */ 

173 

1 

wsr  *  0; 

174 

1 

174 

1 

int  mask  |=  0x03;  /*  unmask  TIMER1,  A2D  DOME,  */ 

175 

1 

imaskl  |=  0x22;  /*  EXTINT,  and  RI  intrpt  V 

176 

1 

176 

1 

176 

1 

/*  Loop  endlessly . */ 

176 

1 

176 

1 

/*  Wait  for  Timerl  Overflow  (every  7.813  msec)  to  start  another  sampling 

176 

1 

sequence.  The  TIMER1  interrupt  handling  routine  starts  the  CHAN  NUM  channel 

176 

1 

sweep:  First  the  last  data  point  is  written  out  to  the  DAC,  then  each  ADC  is 

176 

1 

sampled  on  an  interrupt  driven  basis.  */ 

176 

1 

176 

1 

176 

1 

next: 

176 

1 

176 

1 

/*  Has  the  serial  port  received  a  character?  */ 

176 

1 

176 

1 

if  (  (ri  flag  ==  1)  |  (rtine  flag  ==  99)  )  ( 

177 

2 

int  mask  &=  -0x03;  /*  mask  TIMER1  &  A2D  DONE  intrpts  */ 

178 

2 

imaskl  8=  -0x22;  /*  mask  EXTINT  8  RI  intrpts  */ 

179 

2 

serialO; 

180 

2 

loop  flag  =  0; 

181 

2 

ri  flag  =  0; 

182 

2 

correct  flag  =  0; 

183 

2 

correct  cnt  =  2; 

184 

2 

timerl  flag  =  0; 

185 

2 

ipendl  8=  -0x02; 
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186  2 

187  2 

188  2 

189  2 

190  2 

191  2 

192  2 

192  1 

192  1 

193  1 

194  2 

195  2 

196  2 

197  2 

198  2 

198  1 

198  1 

198  1 

199  1 

200  1 

201  1 

202  1 

202  1 

202  1 

202  1 

202  1 

203  2 

204  2 

205  2 

205  1 

205  1 

206  1 

207  2 

208  2 

209  2 

209  1 

209  1 

209  1 

209  1 

209  1 

210  2 

211  2 

212  2 

212  2 

212  2 

212  2 

212  2 

213  2 

214  2 

214  1 

214  1 

215  1 

216  2 

217  2 

218  2 

218  1 

218  1 


wsr  =  OxOF; 
timerl  =  DELTA  T; 
wsr  =  0; 

intmask  |=  0x03; 
imnskl  |=  0x22; 
goto  wait; 


/*  alternate  window  */ 


/*  unmask  TIMER1  S  A2D  DONE  intrpts  */ 
/*  unmask  EXT INT  &  R!  intrpts  */ 


battlevel  =  0; 

for  (  i=0;  i<3;  i++  )  { 

if  (  battvoltti]  <  LOUBATT  ) 
battlevel  =  300; 

if  <  battvoltCi]  <  ELOWBATT  ) 
battlevel  =  70; 

> 


if  (  restart_flag  ==  1  ) 
goto  restart; 
if  (  correct_flag  1=  1  ) 
goto  wait; 


/*  Load  EOG-VU,  window  the  data,  multiply  by  PRE_SCALE,  &  transform  */ 
for  (  i_temp=0;  i_temp  <  512;  i_temp++  )  { 

xreal Ci_temp]  =  ( int )(((( long)  data_buf fer [2] [buf f_not) [i_temp] )*supgauli 
ximag[i_temp]  =  0; 

> 

fft_for(); 
if  (  error  !=  0  )  ( 
errO; 

goto  restart; 

> 

/*  Save  VU-EOG(w)  */ 


for  (  i_temp=0;  itemp  <  512;  i_temp++  )  { 

eogf [0] [0] [i_temp3  =  (int)  (((long)  xreal [i_temp]  *  10)  /  PRE_SCALE); 
eogf [0] [1] [i_temp]  =  (int)  (((long)  ximag[i_temp]  *  10)  /  PRE_SCALE); 


/*  Load  EOG-VL,  window  the  data,  multiply  by  PRE_SCALE,  8  transform  */ 

xreal Ci_temp]  =  (int)((((long)  data_buffer [3][buff_not] [i_tempj )*supgau[i 
ximag[i_temp]  =  0; 

) 

fftforO; 
if  (  error  !=  0  )  ( 
err(); 

goto  restart; 


temp] )/256); 


temp] )/256); 
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218  1 

218  1 

21S  1 

219  2 

219  2 

220  2 

221  2 

221  2 

221  2 

222  2 

222  2 

223  2 

223  2 

224  2 

225  2 

225  2 

226  2 

227  2 

228  3 

229  3 

230  3 

231  3 

232  3 

233  3 

234  3 

234  2 

235  2 

235  1 

235  1 

235  1 

235  1 

235  1 

235  1 

235  1 

236  2 

237  2 

238  2 

238  1 

238  1 

239  1 

240  2 

241  2 

242  2 

242  1 

242  1 

243  2 

244  2 

245  2 

245  1 

245  1 

245  1 

245  1 

245  1 

245  1 

245  1 

245  1 

246  2 


/*  alpha(w)  =  EOG-VU(w)  /  EOG-VL(w)  ;  alpha(w)  is  scaled  by  10  */ 

for  (  i_temp=0;  i_tenp  <  512;  i_temp++  )  < 

xreal [i_temp]  =  (int)  (((long)  xreal [i_temp]  *  10)  /  PRE_SCALE); 
ximag[i_temp]  =  (int)  (((long)  ximag[i_temp]  *  10)  /  PRE_SCALE); 

lltempl  =  gain_mag  *  ((eogf  [0]  [0]  ti_  temp]  *  (long)  xreal [i_temp] ) 
+  (eogf  (0]  [1]  [i_temp]  *  (long)  ximag[i_temp]  )); 
denom  =  (  (((long)  xreal [i_temp] )  *  xreal  [i_temp] )  + 

(((long)  ximag[i_temp] )  *  ximag[i_temp] )  )  *  10; 

if  (  U tempi  <  0  ) 

l l tempi  =  -lltempl; 

if  (  denom  ==  0  ) 

xrealti_temp]  =  40; 
else  { 

Utemp2  =  1 1  tempi  /  denom; 
if  (  lltemp2  >  40  ) 
lltemp2  =  40; 
if  (  lltemp2  <  10  ) 
lltemp2  =  10; 

xreal Ii_temp]  =  (int)  lltemp2; 

> 

alphati  temp]  =  xreal  [itemp]; 


/*  Load  EOG-H(L-R),  window  the  data,  multiply  by  PRE_SCALE,  &  transform  */ 
for  (  i_temp=0;  i_temp  <  512;  i_temp*+  )  ( 

xreal (i_temp]  =  ( int )(((( long)  data_buf fer [1] tbuf f_not] Ii_temp] )*supgau(i_temp) )/256); 
ximag (itemp]  =  0; 

> 

fft_for(); 
if  (  error  1=  0  )  f 
err(); 

goto  restart; 


for  (  i_temp=0;  i_temp  <  512;  i_temp++  )  { 

eogf  [1]  (OHi  temp]  =  (int)  (((long)  xreal  [i  temp]  *  10)  /  PRE  SCALE); 
eogf [1] [1] [itemp]  =  (int)  (((long)  ximag[i_temp]  *  10)  /  PRE_SCALE); 


/* . V 

/*  Correct  the  EEG  */ 

for  (  i_temp=0;  i_temp  <  SAMP_NUM;  i_temp++  )  ( 
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246  2 

247  2 

248  2 

248  2 

248  2 

248  2 

248  2 

248  2 

248  2 

249  2 

250  2 

251  3 

252  3 

253  3 

254  3 

254  2 

254  2 

254  2 

254  2 

255  2 

256  2 

256  2 

257  2 

258  2 

258  2 

258  2 

258  2 

258  2 

258  2 

258  2 

258  2 

258  2 

259  2 

260  2 

261  3 

262  3 

263  3 

264  3 

264  2 

264  2 

264  2 

264  2 

265  2 

266  2 

266  2 

266  2 

266  2 

267  2 

268  2 

268  2 

268  2 

268  2 

268  2 

269  2 

269  1 

269  1 

269  1 


alphal  =  (alpha ti_temp]  +  10)/2;  /*  alpha  +  1  */ 

alpha2  ®  (alpha Ci_temp]  •  10)/2;  /*  alpha  -  1  */ 

/*  (alpha+1)*D1v  +  (alpha-1)*D2v 

DSGv(w)  =  . .  *  SGv(w) 

(alpha+1 )*D3v  +  (alpha- 1)*D4v 

V 

lltempl  =  (((long)  dv[0] )*alpha1)  +  (((long)  dv[1]  )*alpha2); 
lltemp2  =  (((long)  dv[2] )*alpha1 )  +  (((long)  dv[3]  )*alpha2); 
if  (  lltemp2  ==  0  )  { 
error  =  100; 
err( ); 

goto  restart; 

> 

/*  DSGv(w)  *  EOGvu(w)  */ 

lltemp3  =  1 1  tempi  *  eogf  [0]  [0]  [i_temp] ; 
lltemp4  =  l l tempi  *  eogf  t0H1Hi_temp]  ; 

xreal li_temp]  =  (int)  (  (lltemp3  /  lltemp2)  &  65535  ); 
ximagti  temp]  =  (int)  (  (Utemp4  /  lltemp2)  &  65535  >; 


/*  (alpha+1 )*D1h  +  (alpha- 1)*D2h 

DSGh(w)  =  .  *  SGh(w) 

(alpha+1 )*D3h  +  (alpha- 1)*D4h 

*/ 

lltempl  =  (((long)  dv(41 )*alpha1 )  ♦  (((long)  dv[5) )*alpha2); 
lltemp2  =  (((long)  dv[6] )*alpha1 )  +  (((long)  dv[7] )*alpha2); 
if  (  lltemp2  ==  0  )  { 
error  =  101; 
err(); 

goto  restart; 

> 

/*  DSGh(w)  *  EOGh(w)  */ 

lltemp3  =  lltempl  *  eogf [1] [0] ti_temp] ; 
lltemp4  =  ((tempi  *  eogf Cl] C1J (i_temp] ; 

/*  (OSGv(W)  *  EOuyj(w))  ♦  (DSGh(w)  *  EOGh(tt))  */ 

xreal  [i_temp]  +=  (int)  (  (Utemp3  /  lltemp2)  &  65535  ); 
ximag(i_temp)  +=  (int)  (  (Utemp4  /  lltemp2)  &  65535  ); 


/*  Complex  Conjugate  */ 

ximagti_temp]  =  -ximagli_temp] ; 

) 


/*  Inverse  Transform  */ 
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269  1 

269  1 

270  1 

271  2 

272  2 

273  2 

273  1 

273  1 

273  1 

273  1 

274  2 

275  2 

275  1 

275  1 

275  1 

275  1 

275  1 

276  2 

276  2 

277  2 

277  1 

277  1 

277  1 

277  1 

277  1 

278  2 

279  2 

279  1 

279  1 

280  1 

280  1 

280  1 

281  1 

282  1 

282  1 

283  1 

283  1 

283  1 

283  1 

284  1 
286 

286 

286 

286 

286 

286 

286 

286 

287  1 

288  1 

289  1 

290  1 

291  1 
293 

293 

293 

293 


fft_for();  /*  inverse  FFT  */ 

if  <  error  1=  0  )  C 
err(); 

goto  restart; 


/*  Remove  Window  */ 

for  (  i_temp=BUFF_START;  itemp  <  (SAMP_NUM-BUFF_START);  i_temp++  )  C 

xreal ti_temp]  =  (int)(<  ((long)  xreal [i_temp] )  *  83886)  /  supgauti_temp] ); 

) 


/*  Moving  Average  Filter  (7  point)  of  corrector  */ 


for  (  i_temp=BUFF_START;  itemp  <  (SAMPNUM-BUFFSTART );  i_temp++  )  ( 

alphali_temp]  =  (xreal [i_temp-33  +  xreal [i_temp-2]  +  xreal [i_temp-1J  +  xreal [i_temp] 


) 


xreal [i_temp+1]  +  xreal ti_temp+2)  +  xreal [i_temp+3] )  /  7; 


/*  EEGcorr(t)  =  EEGobs(t)  -  IFFTC  (SGv  *  EOGvu)  +  (SGh  *  EOGh)  >  */ 

for  (  i_temp=BUFF_START;  i_temp  <  (SAMP_NUM-BUFF_START);  i_temp++  )  C 

out_buffer[buff_not] ti_temp]  =  data_buffer[0] [buff_not] ti_temp]  -  alphalitemp) ; 

> 

correct_f log  =  0; 
wait: 

loop_f lag- - ; 

while  (  timer1_flag  ==  0) 

t 

timer1_flag  =  0; 

/*  Loop  again  */ 

goto  next;  /*  wait  for  another  TIMER1  INTERRUPT  */ 

> 


/* . */ 

/*  nmi_int  Function  -  NMI  interrupt  handler  */ 

void  nmi_int(void) 

< 

imaskl  =  0x00;  /*  disable  interrupts  */ 

intjnask  =  0x00; 
error  =  3; 
err(); 

> 


/*  error  */ 
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293 

293 

293 

293 

293 

293 

293 

294  1 

295  1 
297 

297 

297 

297 

297 

297 

297 

297 

297 

297 

298  1 

299  1 

300  1 
302 

302 

302 

302 

302 

302 

302 

302 

302 

302 

302 

303  1 

304  1 

304  1 

304  1 

304  1 

305  2 

306  2 

307  2 

308  2 

308  1 

309  2 

310  2 

311  2 

312  2 

312  1 

312  1 

313  2 

314  2 

315  2 

315  1 

315  1 

315  1 

315  1 

316  1 


. . . 

/*  extint  Function  -  EXTINT  interrupt  handler  */ 

/*  -  senses  an  external  trigger  on  TRI  */ 

void  extint(void) 

< 


trocounter  =  TRO_PULSE;  /*  load  TRO  counter  */ 

> 


/*  . V 

/*  Receive  Function  -  Serial  Port  Receive  interrupt  handler  */ 

void  receive(void) 
t 


ri_flag  =  1; 
getcharO; 

> 


/*  TIME1  •  TIMER1  overflow  interrupt  handler 
•  Operates  in  conjunction  with  SAMP  */ 

void  timel(void) 

t 

fnt  i,  j,  k; 

/*  Reset  TIMER1  */ 

if  (  rtineflag  1=1  )  € 
wsr  =  OxOF; 
timerl  =  DELTA_T; 
wsr  =  0; 

) 

else  t 

wsr  =  OxOF; 

timerl  =  DEITA_T_FAST; 
wsr  =  0; 

) 

if  <  (loopflag  >  0)  &  <rtine_flag  !=  3>  )  < 
error  =  88; 
err(); 

> 

/*  Is  TRO  high  ?  */ 

if  (  trocounter  >  0  ) 
tro_counter--; 
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317  1 

317  1 

313  1 

319  1 

320  1 

321  2 

322  3 

323  3 

324  3 

325  3 

326  3 

326  2 

326  1 

326  1 

326  1 

326  1 

327  2 

328  2 

329  2 

330  2 

330  2 

331  3 

332  3 

333  3 

333  2 

335  3 

336  3 

337  3 

337  2 

339  3 

340  3 

341  3 

341  2 

342  2 

34?  2 

343  3 

344  3 

345  3 

345  3 

346  4 

347  4 

348  4 

349  5 

350  5 

351  5 

352  5 

353  5 

354  5 

354  4 

355  4 

355  4 

356  5 

357  5 

358  5 

358  4 

359  4 

360  4 


battcnt++; 

if  (  battcnt  >=  32000  ) 

battcnt  =  0;  , 

if  (  bctt level  1  =  0  )  ( 

if  (  (battcnt  X  battlevei)  ==  0  >  ( 

if  (  ((battcnt  /  battlevei)  X  2)  ==  0  ) 
pwm_control  =  pmwm; 

else 

pwm  control  =  0; 

} 

) 

/*  Interrogate  ?  */ 

if  (  rtineflag  ==  1)  { 

i  =  CALNUMM1  -  cal_counter; 

gssunlti]  +=  data  buf fertKcal  modeX2)/c8l_mode)*2)  [buffnum]  [new_pt-1]; 
gssum2[ij  +=  data_buffert(cal_mode%2)+2-(cal_mode/2)] [buffnun]  (new  pt-IJ; 

if  (  calcounter  >  CALON  )  { 

importl  &=  ~(stim_port(3-cal_mode] ); 
ioportl  =  importl; 

> 

else  if  (calcounter  >  CAL_OFF  )  { 

importl  |=  (stim_port(3-cal_mode] ); 
ioportl  =  inportl; 

> 

else  if  (calcounter  >=  0  )  { 

importl  &=  ~(stim_portt3-caljnode] ); 
ioportl  =  importl; 

> 

cal_counter- 

if  (calcounter  <  0)  ( 

cal  counter  =  CAL_NUM_M1; 
rep_count -  - ; 

if  (  repcount  ==  0  )  { 

intjnask  *=  -0x01;  /*  mask  TIMER1  intrpt  */ 

cal_mode- - ; 

if  (  calmode  ==  0  )  ( 
rtine_flag  =  2; 
tune_pt  =  0; 
new_pt  =  BUFFSTART; 
buff_num  =  0; 
buff_not  =  1; 

> 

reP_count  -  CAL_REP;  /*  reset  average  counter  */ 

for  (  i=0;  i  <  CALNUM;  i++)  <  /*  average  */ 

gs[2-cal_mode] CO] ti]  =  (gssumlCil  /  CALREP); 
gs  12  cal  mode]  [1]  ti]  =  (gssum2[i]  /  CAL  REP); 

) 

templ  reg  =  (gs [2 -cal  mode] CO] [CAL  NUM  M1]  +  gs [2-csl_mode] 10] 10] )  /  2; 
temp2_reg  =  (gs(2-cal_modeJ [1] ICAL_NUM_M1]  +  gs 12 -cal  mode] [1] [0] )  /  2; 
for  (  i=0;  i  <  CALNUM;  i++)  ( 
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361 

362 

363 

364 

365 
365 

365 

366 

367 

368 
368 
368 
368 

368 

369 

370 

371 

372 
372 

374 

375 

376 
376 

378 

379 

380 

380 

381 

381 

382 

383 

384 

385 
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gsC2-caljnode] [0] til  *=  temp1_reg; 
gst2-cal  jnode]  113  li]  -=  temp2_reg; 
gssumlli]  =  0; 
gssum2[i]  =  0; 

> 

> 

loop_flag  =  0; 
intjnask  |=  0x01; 
goto  skiploop; 

> 

> 


if  (  rtine_flag  ==  7  )  < 

if  (  cal_cnt  >  100  )  < 

importl  &=  ~(stim_port[2J ); 
ioportl  =  importl; 

> 

else  if  (cal_cnt  >  90  )  { 

importl  |=  (stim_port[23 ); 
ioportl  =  importl; 

) 

else  if  (cal_cnt  >=  0  )  { 

importl  8=  ~(stim_port[2] ); 
ioportl  =  importl; 

> 

calcnt--; 

if  (cal_cnt  <  0}  { 
calcnt  =  150; 
caljnonitor--; 
if  (  cal  monitor  ==  0  )  ( 
intjnask  &=  -0x01; 
rtineflag  =  0; 

pr i nt f ( "\n\r\n\r\n\r  REAL  TIME  MONITORING. . .\n\r\n\r\n\r"); 

loopflag  =  0; 
intjnask  |=  0x01; 
goto  skiploop; 

> 

) 

> 


1  /*  Initiate  the  first  A/0  Conversion  */ 

1 

1  chan  =  0;  /*  load  ADC  chan  offset  */ 

1  intjnask  |=  0x02; 

1  adcommand  =  CHAN_8ASE;  /*  start  1st  ADC  conversion  */ 

1  loop_flag  =  2; 

1  skiploop: 

1  timer1_flag  =  1; 

1  ) 
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/*  SAMP  •  A2D  CONVERSION  COMPLETE  interrupt  handler. 

•  Operates  in  conjunction  with  T1ME1. 

-  The  last  sample  in  the  data_buffer  is  written  out  to  the 
corresponding  DAC,  converted  to  8  bits. 

-  reads  the  sample  from  the  ADC. 

•  returns  a  10  bit  value. 

-  start  a  conversion  on  the  next  channel  of  the  sweep.  */ 
void  samp(void) 


register  int  temp_reg; 

/*  Write  out  the  data  to  the  DAC 

-  convert  raw  data  from  10  bits  to  8  bits;  */ 

if  (  chan  ==  0  )  { 

dac(  0,  (data_buffer[0] [buff_num] [new_pt-1]  »  2)  ); 
if  (  rtine_flag  ==  3  )  C 

dac(  1,  <data_buffer[4] [buff_numj [new_pt-1J  >>  2)  ); 
dacC  2,  (out_buffer[buf f_num] [new_pt]  >>  2)  ); 
dac(  3,  (data_bufferC4][buf f_num] [newpt]  »  2)  ); 
dac(  4,  (data_buffer[0] [buf f_num] [new_pt]  »  2)  ); 
dac(  5,  (data_buffer[1] [buff_num] [newpt]  »  2)  ); 
dac(  6,  (data_buffer[2] [buff_num] [newpt)  »  2)  ); 
dac(  7,  (databuffer[3] Cbuffnum] [newpt]  »  2)  ); 

> 

else  { 

dacC  1,  (data_buffer[1] [buff_numj [newpt- 1]  »  2)  ); 
dac(  2,  (data_buffer(2] [buff_num]  [new_pt-1]  »  2)  ); 
dac(  3,  (data_buffer(3] [buff  nun]  [new  pt - 13  »  2)  ); 

> 

> 


/*  Read  a  new  sample  from  the  ADC;  convert  from  0  -  1023  to  512  */ 

if  (  chan  <  4  )  { 

tempreg  =  ad_result_hi ; 

tempreg  =  (((  tempreg  «  8)  +  ad_result_lo)  »  6)  -  512; 
data_buffer[chan]  [buf  f_nun]  [new_pt]  =  temp_reg; 

> 

else  if  (  chan  <  7  )  { 

tempreg  =  ad_result_hi; 

temp_reg  =  ((  temp_reg  «  8)  +  ad_result_lo)  »  6; 
batt_volt [chan-4]  =  temp_reg; 

} 

if  (  chan  ==  3  )  { 

data_buffer[4] [buff_num] [new_pt]  =  -200; 
if  (  tro_counter  >  0  ) 

data_buffer[4] [buff_num] [new_pt]  =  200; 


/*  If  necessary,  store  it  in  the  other  buffer  as  well  */ 

if  (  (rtineflag  ==  2)  &  (chan  <  CHAN_NUK)  ) 

data_tuner[chan] [tune_pt]  =  databuffer [chan] [buf f_nun] [new_pt] ; 
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if  (  newpt  <  LEFTJ.1M  ) 

data_buf fer tchanltbuf f_not] [new_pt+BUFF_DI FF]  =  temp_reg; 
else  if  (  new  pt  >  RIGHT_UM  ) 

datajxiffertchen]  [buff_not]  tnew_pt-BUFF_DIFF]  =  temp_reg; 

/*  Either  start  a  new  conversion,  or  update  the  data  pointers  */ 


chan_end  =  8; 
if  (  rtineflag  ==  1  ) 
chan_end  =4; 
chan++; 

if  (  chan  <  chan_end  ) 

adcoronand  =  CHAN_BASE  +  chan;  /* 

else  C  /*  Update  data 

if  (  newpt  <  ( S AMP_NUM -BUFF_START)  )  { 
new_pt++; 
tmept++; 

if  (  (tune_pt  ==  SAMP_X2)  &  (rtine_flag  == 
rtine_flag  =  99; 
if  (  rtine_flag  ==  3  )  < 
if  (  newpt  ==129  )  ( 

if  (  correct_flag  ==  1  )  { 
error  =  80; 
err(); 

> 

else  ( 

if  (  correct_cnt  >  0  ) 
correct_cnt- -; 
if  (  correct_cnt  ==  0  ) 
correct_flag  =  1; 

> 


} 

> 

else  { 

newpt  =  BUFFSTART; 
tempreg  =  buffnot; 
buffnot  =  buffnvm; 
buff_num  =  temp_reg; 

> 


loopflag--; 
intmask  &=  -0x02; 


start  next  conversion  */ 
buffer  pointers  */ 


2)  ) 


) 


/* . */ 

/*  Serial  Port  Conmuni cat  ions  Function  */ 


void  serial(void) 
t 

int 


resp,  i,  j,  k,  resptmp,  jk,  tempwait; 
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467 

467 

467 

1 

1 

1 

int_mask  =  0; 

468 

1 

pwm_control  =  pmwm; 

469 

1 

if  (  rtine_flag  ==  98  ) 

470 

1 

goto  menuu; 

471 

1 

if  (  rtine_flag  ~  99  )  { 

472 

2 

printf ("Press  'PgDn',  select  '7',  and 

' .\n\r" 
473 

>; 

2 

printf ("Interrogation  data  will  take  fi 

474 

2 

sndgs: 

474 

2 

resp  =  getcharO; 

475 

2 

printf("Xc",resp); 

476 

2 

resp  -=  65; 

477 

2 

if  (  resp  !=  0  ) 

478 

2 

goto  sndgs; 

479 

2 

printf("\n"); 

480 

2 

for  (  i=0;  i<3;  i++  )  { 

481 

3 

for  (  j=0;  j<2;  j++  )  { 

482 

4 

for  (  k=0;  k<CAL_NUM;  k++  )  { 

483 

5 

senddata(gsti)  [j]  [k] ); 

484 

5 

for  (  jk=0;  jk<10000;  jk++ 

tempwait  =  0; 


for  (  i=0;  i <CHAN_NUM;  i++  )  < 
for  (  j=0;  j<SAMP_X2;  j++  )  ( 
senddata(data_tunerti]  tj] ); 
for  (  jk=0;  jk<10000;  jk++  ) 
tempwait  =  0; 

> 

> 

printf ("Push  'ESC',  ALT-X,  'Y',  and  then  run  the  program  EEG.\n\r"); 
goto  menuu; 


menu: 

pr i nt f ( "\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r 


GMS  Engineering  Corporation\n\r"); 


printfC  EEG  Artifact  Rejection  System\n\r\n\r\n\r\n\r\n\r"); 

printfC  N  -  Channel  Number  Selection\n\r"); 

printfC  L  •  LEO  Light  Level\n\r"); 

printfC  R  -  Real  Time  Monitoring\n\r"); 

printfC  P  -  Calibration  Pulses\n\r"); 

printfC  I  -  Interrog8tion\n\r"); 

printfC  C  -  Correction\n\r\n\r\n\r\n\r\n\r“); 

printfC  Enter  RESPONSE  >  "); 

menuu: 

resp  =  getcharO; 
printf("Xc",resp); 

if  <  (resp  ==  * L • )  |  (resp  ==  'l')  )  C 
lvl4: 

resptmp  =  pmwm  +  1; 

printf("\n\rCurrent  light  level  is:  %d\n\rEnter  new  level  (1-255)  >  ".resptmp); 
resptmp  =  0; 
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508 

2 

resp  =  getcharO; 

509 

2 

printf ("%cn, resp); 

510 

2 

if  (  resp  ==  13  ) 

511 

2 

goto  Ivll; 

512 

2 

if  (  (resp<48)  |  (resp>=58)  ) 

513 

2 

goto  Ivltt; 

514 

2 

resptmp  =  (resptmp*10)  +  (resp-48); 

515 

2 

goto  Ivlll; 

516 

2 

Ivll: 

516 

2 

if  (  resptmp  >  255  ) 

517 

2 

goto  lvl4; 

518 

2 

if  (  resptmp  ==  0  ) 

519 

2 

resptmp  =  pmwm  +  1; 

520 

2 

resptmp- - ; 

521 

2 

pmwm  =  resptmp; 

522 

2 

pwm_control  =  pmwm; 

523 

2 

> 

523 

1 

523 

1 

else  if  (  (resp  ==  'R')  |  (resp  ==  'r')  >  * 

525 

2 

pr int f ( "\n\r\n\r\n\r  REAL  TIME  MONITORING. .  An\r\n\r\n\r" 

526 

2 

rtine_flag  =  0; 

527 

2 

goto  menuend; 

528 

2 

) 

528 

1 

528 

1 

else  if  (  (resp  ==  'P')  1  (resp  ==  'p')  )  { 

530 

2 

pr i nt f ( "\n\r\n\r\n\r  Calibra*'on  Pulses\n\r\n\r\n\r"); 

531 

2 

rtine_flag  =  7; 

532 

2 

cal_monitor  =  300; 

533 

2 

calcnt  =  150; 

534 

2 

goto  menu_end; 

535 

2 

> 

535 

1 

535 

1 

else  if  (  (resp  ==  '!•)  1  (resp  ==  • i • >  )  { 

537 

2 

cal  counter  =  CAL_NUM  Ml; 

538 

2 

calmode  =  3; 

539 

2 

rep_count  =  CAL_REP; 

540 

2 

rtineflag  =  1; 

541 

2 

for  (  i=0;  i<CAL_NUM;  i++  )  C 

542 

3 

gssumlli]  =  0; 

543 

3 

gssun2[ij  =  0; 

544 

3 

> 

544 

2 

pr i nt f ( "\n\r\n\r\n\r  I NTERROGAT I NG . . . \n\r\n\r\n\r" ) ; 

545 

2 

goto  menuend; 

546 

2 

> 

546 

1 

else  if  (  (resp  ==  'C')  1  (resp  ==  'c')  )  { 

548 

2 

printf ("\n\r\n\r"); 

549 

2 

printf("Press  'PgUp',  select  '7',  and  then  type  'correct' .\n\r") 

550 

2 

sgvmag  =  recvdataO; 

551 

2 

sgh_mag  =  recvdataO; 

552 

2 

gainjnag  =  recvdataO; 

553 

2 

for  (  1=0;  i<8;  i++  ) 

554 

2 

dv[i]  =  recvdataO; 

555 

2 

dvtO]  =  ((int)  (((long)  dv[0]  *  sgv  mag)  /  10))  &  65535; 

556 

2 

dv[4]  =  ((int)  (((long)  dv(4)  *  sghmeg)  /  10))  S  65535; 

557 

2 

rtine_flag  =  3; 

558 

2 

pr i nt?( "\n\r\n\r\n\r  CORRECTING. . .\n\r\n\r\n\r"); 

559 

2 

resp  =  getcharO; 

J 
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", resptmp); 
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569  2 

570  2 

571  2 

572  2 

573  2 

574  2 

574  2 

575  2 

576  2 

577  2 

578  2 
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583  1 
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585 
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587  1 

588  1 
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589  1 

590  1 

591  2 

592  2 

593  2 

593  1 

594  1 

595  1 

596  2 

597  2 

598  2 

599  2 


goto  menu_end; 

> 

else  If  <  (resp  ==  • M • }  |  (resp  ==  'n')  )  { 
cher>4: 

resptmp  =  (importl  »  3)  +  1; 

printf ("\n\rCurrent  channel  number  is:  Xd\n\rEnter  new  charnel  nunber  (1-24)  > 
resptmp  =  0; 

channn: 

resp  =  getcharO; 
pr i nt  f ( "%c" , resp) ; 
if  (  resp  ==  13  ) 
goto  chann; 

if  (  (resp<48)  |  (resp>=58)  ) 
goto  channn; 

resptmp  =  <resptmp*10)  +  (resp-48); 
goto  channn; 

chann: 

if  <  resptmp  -  24  > 
goto  chan4; 
if  (  resptmp  ==  0  ) 

resptmp  =  (importl  »  3)  +  1; 
resptmp- - ; 
importl  &=  0x07; 
importl  |=  (resptmp  «  3); 
ioportl  =  importl; 

> 

goto  menu; 
menu  end: 


/* . */ 

/*  SENDOUT  DATA  Function  */ 

void  senddata(data_out) 

int  data_out; 

< 

int  i,  j,  temp_data,  temp_sign,  data2_out,  data_first; 

char  temp_chr; 

temp_sign  =  0; 
if  (  data_out  <  0  )  ( 
temp_sign  =  1; 
data  out  =  -date_out; 

) 

data2_out  =  0; 
data_first  =  data_out; 
for  (  i=4;  i>=0;  i--  )  ( 
data^first  -=  data2_out; 
temp_data  =  data_first; 
for  (  j  =  i ;  j>0;  j -  -  ) 
temp_data  /=  10; 
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631 

631 


data2_out  =  temp_data; 
for  (  j=i;  j>0;  J--  ) 
data2_out  *=  10; 
printfC'Xd",  temp  data); 

> 

printf ("%d\n",  temp_sign); 


/* . 

/*  RECEIVE  DATA  Function  */ 

int  recvdata(void) 

int  i,  j,  temp_data,  data_in; 

char  temp_chr; 

temp_data  =  0; 
for  (  i=4;  i>=0;  i--  )  { 
retryl : 

temp_chr  =  getcharO; 
if  (  tempchr  ==  32  ) 
temp_chr  =  48; 

if  (  (tempchr  <  48)  |  (temp  chr  >=  58)  ) 
goto  retryl; 

data_in  =  (int)  tempchr  -  48; 
for  (  j=0;  j<i ;  j*+  ) 
data_in  10; 
temp  data  +=  data  in; 

> 

retry2: 

tempchr  =  getcharO; 

if  (  (temp  chr  <  48)  |  (temp  chr  >  49)  ) 
goto  retry2; 
if  (  temp  chr  ==  49  ) 

temp_data  =  -temp_data; 
return  (  tempdata  ); 


/*-•-- 

/* 

DAC  Function 

*/ 

void 

dac(channe(. 

level) 

int 

channel. 

level; 

( 


*/ 
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648  1 

649  1 

650  1 

650  1 

651  2 

652  2 

652  1 

653  1 

654  1 
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655  1 

656  1 

657  1 

658  1 

659  1 

660  1 

661  1 

662  1 

663  1 

664  1 


int  dace; 

/*  Decode  DAC  address  */ 

if  (  channel  <=  3  ) 

dace  =  0xE013  +  channel; 

else 

dace  =  0xE007  +  channel; 

if  (  (dace  %  2)  ==  0) 
dace  -=  2; 

/*  Write  out  the  sample  to  the  DAC  */ 

if  (  level  >  127  ) 
level  =  127; 

else  if  (  level  <  -128  ) 
level  =  -128; 

memset(  dace,  (level+128),  1  );  /*  +-  128  into  0  -  256  */ 

> 


/* . */ 

/*  Error  Function  */ 

void  err(void) 

( 

int  i; 

int  mask  =  0x00;  /*  mask  all  interrupts  */ 

imaskl  =  0x02; 

for  (  i=0;  i  <  CHANNUM;  i++)  ( 

dac<  i,  0  );  /*  zero  DACs  */ 

printf(cls);  /*  clear  the  screen  */ 

printf("\n\n\r\t\t\t\t%cERROR\n\n\r",bel  l); 
printf("\t\t\t\tNo.  Xi\n\r", error); 

intjnask  &=  -0x03;  /*  mask  TIMER1  &  A2D  DONE  intrpts  */ 

imaskl  &=  -0x22;  /*  mask  EXTINT  8  RI  intrpts  */ 

rtine_flag  =  98; 

serialO; 

ri_flag  =  0; 

ipendl  &=  -0x02; 

int_mask  |=  0x03;  /*  unmask  TIMER1  &  A2D  DONE  intrpts  */ 

imaskl  |=  0x22;  /*  unmask  EXTINT  S  RI  intrpts  */ 

restart_flag  =  1; 
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MODULE  INFORMATION: 


CODE  AREA  SIZE 

=  16DFH 

5855D 

CONSTANT  AREA  SIZE 

=  080EH 

20620 

DATA  AREA  SIZE 

=  6FDBH 

286350 

STATIC  REGS  AREA  SIZE 

=  005AH 

900 

OVERLAY ABLE  REGS  AREA  SIZE 

*  OOOAH 

10D 

MAXIMUM  STACK  SIZE 

=  0094 H 

1480 
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0  WARNINGS 


0  ERRORS 
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1 


DOS  3.30  (038-N)  MCS-96  MACRO  ASSEMBLER,  VI. 2 

SOURCE  FILE:  FFT_FOR.A96 
OBJECT  FILE:  FFT_F0R.08J 

CONTROLS  SPECIFIED  IN  INVOCATION  COMMAND:  <none> 

ERR  LOC  OBJECT  LINE  SOURCE  STATEMENT 

1  FFT_FO  MOOULE  STACKS! ZE(6) 

2 

3  ;FFT  ALGORITHM  FROM  INTEL  APPLICATIONS  NOTE,  AP-275,  BY  IRA  HORDON 

4  ;  "EMBEDDED  CONTROL  APPLICATIONS",  INTEL  CORP,  1988. 

5 

0000  6  RSEG 

7  EXTRN  error 

8 

0024  9  OSEG  at  24H 


0024 

10 

TMPR: 

dsl 

1 

0028 

11 

TMPI : 

dst 

1 

002C 

12 

TMPR1: 

dsl 

1 

0030 

13 

TMPI 1 : 

dsl 

1 

0034 

14 

XRTMP: 

dsl 

1 

0038 

15 

XITMP: 

dsl 

1 

003C 

16 

WRP: 

dsw 

1 

003E 

17 

WIP: 

dsw 

1 

0040 

18 

PUR: 

dsw 

1 

0042 

19 

IN  CNT : 

dsw 

1 

0044 

20 

NDIV2: 

dsw 

1 

21 

0046 

22 

KPTR: 

dsw 

1 

0048 

23 

KN2: 

dsw 

1 

004A 

24 

N  SUB  K: 

dsw 

1 

004C 

25 

RK: 

dsw 

004 E 

26 

RNK: 

dsw 

1 

0050 

27 

SHFT  CNT: 

dsb 

1 

0051 

28 

LOOP  CNT: 

dsb 

1 

29 

30 


0000  31  DSEG 

32 

33  EXTRN  XREAL,  XI MAG 

34 

35  ;  XREAL,  XIMAG:  Base  addresses  for  512  16-bit  signed 

36  ;  entries  for  real  and  imaginary  numbers,  respectively. 

37 

38  SEJECT 
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LOC  OBJECT 

LINE 

SOURCE 

STATEMENT 

39 

0000 

40 

CSEG 

41 

42 

PUBLIC  fft 

for  ;  Starting  point  for  FFT  algorithm 

43 

44 

1 

;;;;  START  FOURIER  CALCULATIONS 

0000 

45 

FFT_for: 

;;;;  400  '  INITIALIZATION  OF  LOOP 

0000  1100 

E 

46 

clrb 

error 

47 

0002  FC 

48 

clrvt 

0003  810151 

49 

Idb 

loop_cnt,  #1 

0006  B 10850 

50 

tdb 

shft  cnt,  #8 

0009  A1 000244 

51 

Id 

ndiv2,  #512 

52 

1 

;;;;  410  K=0 

000D 

53 

OUT_LOOP: 

0000  0146 

54 

clr 

kptr 

55 

1 

;;;;  420  IF  LOOP  >  EXP  THEN  700 

OOOF  990951 

56 

cmpb 

loop  cnt,  #9  ;  512  =  2*9 

0012  DA0220A3 

57 

! 

bgt 

UNWEAVE 

58 

59 

0016 

60 

MID_LOOP: 

;;;;  430  INCRNT=0 

0016  0142 

61 

clr 

i n_cnt 

62 

63 

;;;;  440  'CALCULATIONS  BEGIN  HERE 

0018 

64 

I NLOOP : 

0018  65020042 

65 

add 

in  cnt,  #2  ;  450  INCNT=INCNT+1 

66 

;;;;  460  P=BR(INT(K/(2ASHIFT)>) 

001C  A04640 

67 

Id 

pwr,  kptr 

001 F  085040 

68 

shr 

pwr,shft  cnt  ;  Calculate  mult  factors 

0022  71FE40 

69 

andb 

pwr,  #1 11111 10b 

0025  A341FC0040 

R 

70 

(d 

pwr,  brevfpwrj 

71 

? 

;;;;  470  WRP=WR(P)  :  UIP=WI(P)  :  KN2=X+N2 

002A  A341FC043C 

R 

72 

gw: 

Id 

wrp,  wrtpwr] 

002 F  A341FE083E 

R 

73 

Id 

wip,  wi [pwr] 

0034  44444648 

74 

add 

kn2,  kptr,  ndiv2 

75 

76 

» # 

Complex  multiplication  follows 

77 

78 

# 

;;;;  480  TMPR=(WRP*XR(KN2)-UIP*XI(KN2))/2 

0038  FE4F4900003C24 

E 

79 

gm: 

mul 

tmpr,  wrp,  xreal Ikn2] 

003 F  FE4F4900003E28 

E 

80 

mul 

tmpi,  wip,  ximagfkn2] 

0046  682A26 

81 

sub 

tmpr+2,  tmpi +2 

82 

r 

;;;;  490  TMPI=  (URP*XI (KN2)+WIP*XR(KN2) )/2 

0049  FE4F4900003 C2C 

E 

83 

mul 

tmprl,  wrp,  ximag[kn2J 

0050  FE4F4900003E28 

E 

84 

mul 

tmpi,  wip,  xreal Ikn2] 

0057  642E2A 

85 

add 

tmpi+2,  tmprl+2 

MCS-96  MACRO  ASSEMBLER  FFT  FO 


04/06/90  16:00:52  PAGE  3 
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LIME 

SOURCE 
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86 

87 

;;  high  byte  only  of  a  signed  multiply 

88 

;;  provides  an  effective  divide  by  two 

89 

005A  DC55 

90 

BVT 

ERR1  ;  branch  on  error 

91 

005C  A34700002C 

E 

92 

Id 

tmprl,  xreal  tkptr]  ;  500  TMPR1=XR(K)/2 

0061  0A012C 

93 

shra 

tmprl,  #1  ;  TMPI  1=XI 00/2 

0064  A34 7000030 

E 

94 

Id 

tmpil,  ximag tkptr] 

0069  0A0130 

95 

shra 

tmpil,  #1 

96 

97 

t 

;;;;  510  XROCN2)  =  TMPR1  =  TMPR1  -TMPR 

006C  48262C34 

98 

gr2: 

sub 

xrtmp,  tmprl,  tmpr+2 

0070  C349000034 

E 

99 

St 

xrtmp,  xreal [kn2] 

100 

« 

;;;;  520  XI (KN2)  =  TMPR1  =  TMPI 1  -TMPI 

0075  482A3038 

101 

gx2: 

sub 

xitmp,  tmpil,  tmpi+2 

0079  C349000038 

E 

102 

St 

xitmp,  ximag[kn2] 

103 

! 

;;;;  530  XR(K)  =  TMPR1  -  TMPR 

007E  44262C34 

104 

add 

xrtmp,  tmprl,  tmpr+2 

0082  C347000034 

E 

105 

St 

xrtmp,  xreal tkptr] 

106 

r 

;;;;  540  XI OO  =  TMPI1  +  TMPI 

0087  442A3038 

107 

gx: 

add 

xitmp,  tmpil,  tmpi+2 

008B  C347000038 

E 

108 

St 

xitmp,  ximag tkptr] 

109 

0090  DC23 

110 

BVT 

ERR2  ;  Branch  on  error 

111 

112 

# 

;;;;  560  k  =  x  +  i 

0092  65020046 

113 

ik: 

add 

kptr,  #2 

114 

115 

• 

;;;;  570  IF  incnt  <  n2  then  goto  450 

0096  884442 

116 

emp 

in  ent,  ndiv2 

0099  D6022778 

117 

} 

bit 

IN^LOOP 

118 

119 

• 

;;;;  580  K  =  K  +  N2 

0090  644446 

120 

add 

kptr,  ndiv2 

121 

» 

;;;;  590  IF  K  <  N1  THEN  GOTO  430 

OOAO  89"E0346 

122 

emp 

kptr,  #1022  ;;  N1  =  2  *(N-1) 

00A4  D602276E 

123 

! 

bit 

HID_LOOP 

124 

125 

# 

;;;;  600  loop  =  loop  +  1  :  N2  =  N2  /  2 

00A8  1751 

126 

ineb 

loop  ent  ;  605  SHIFT  =  SHIFT  +  1 

OOAA  0A0144 

127 

shra 

ndiv2,  #1 

OOAD  1550 

128 

deeb 

shft_cnt 

129 

* 

;;;;  610  GOTO  400 

OOAF  275C 

130 

br 

OUT_LOOP 

131 

132 

00B1  B10100 

E 

133 

ERR1: 

tdb 

error,  #01  ;  overflow  error 

0084  FO 

134 

ret 

00B5  B10200 

E 

135 

ERR2: 

Idb 

error,  #02  ;  overflow  error 

0088  FO 

136 

ret 

137 

138 

SEJECT 
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139 

140 

1 

00B9 

141 

UNWEAVE: 

142 

143 

f 

00B9  0146 

144 

clr 

00BB  A 100044 A 

145 

Id 

146 

OOBF 

147 

UNLOOP: 

148 

OOBF  A347FC004C 

R 

149 

Id 

150 

00C4  884C46 

151 

cmp 

00C7  D628 

152 

bge 

153 

00C9  A34 7000024 

E 

154 

Id 

OOCE  A347000028 

E 

155 

Id 

0003  A34D00002C 

E 

156 

Id 

00D8  A34D000030 

E 

157 

Id 

158 

OODD  C34D000024 

E 

159 

St 

OOE2  C34D000028 

E 

160 

St 

00E7  C34700002C 

E 

161 

St 

OOEC  C347000030 

E 

162 

St 

163 

164 

• 

1 

OOFI  65020046 

165 

ENDL:  add 

00F5  6902004A 

166 

sub 

00F9  D7C4 

167 

bne 

168 

OOFB  FO 

169 

RET 

170 

171 

;$nol ist 

OOFC 

172 

CSEG 

173 

OOFC 

174 

BREV: 

175 

OOFC  0000000200010003 

176 

DCU 

2*0,  ; 

0108  8001800340004002 

177 

DCW 

2*192 

0114  C0OOCOO2COO 1 C003 

178 

OCW 

2*96, 

0120  2001 2003A000A002 

179 

DCU 

2*144 

012C  6000600260016003 

180 

DCU 

2*48, 

0138  E001E003 10001 002 

181 

DCU 

2*240 

0144  9000900290019003 

182 

DCU 

2*72, 

0150  5001 5003D 0000002 

183 

DCW 

2*168 

015C  3000300230013003 

184 

DCU 

2*24, 

0168  BOO 1 B00370007002 

185 

DCW 

2*216 

0174  F000F002F001 F003 

186 

DCU 

2*120 

0180  0801080388008802 

187 

DCU 

2*132 

018C  4800480248014803 

188 

DCU 

2*36, 

0198  C801 C80328002802 

189 

OCW 

2*228 

01A4  A800A802A801 A803 

190 

DCU 

2*84, 

01B0  680 1 6803E800E802 

191 

DCW 

2*180 

01BC  1800180218011803 

192 

DCW 

2*12, 

01C8  9801980358005802 

193 

DCU 

2*204 

01D4  D800D802D801D803 

194 

DCU 

2*108 

01E0  380 1 3803B800B802 

195 

DCU 

2*156 

;;;;  700  •  REORDERING  STARTS  HERE 

;;;;  720  FOR  K  =  0  TO 

kptr 

n_sub_k,  #1024 

;;  Bit  reversal  of  the  transformed  array 

rk,  brevlkptr] 

kptr,  rk 
ENDL 

tmpr,  xreaUkptr] 
tmpi,  xi mag [kptr] 
tmprl,  xreallrk] 
tmpil,  ximaglrk] 

tmpr,  xrealtrk) 
tmpi,  ximagtrk] 
tmprl,  xrealtkptr] 
tmpil,  xi mag [kptr] 


kptr,  #2 


;;  950  NEXT  K 


n_sub_k,  #2 
UN  LOOP 
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01EC  7800780278017803  196  DCW 

01F8  F801 F80304000402  197  DCW 

0204  8400840284018403  198  DCW 

0210  44014403C400C402  199  DCW 

021C  2400240224012403  200  DCW 

0228  A40 1 A40364006402  201  DCW 

0234  E400E402E401E403  202  DCW 

0240  1401140394009402  203  DCW 

024C  5400540254015403  204  DCW 

0258  D401D40334003402  205  DCW 

0264  B4008402B401B403  206  DCW 

0270  74017403F400F402  207  DCW 

027C  0COOOCO2OCO1OC03  208  DCW 

0288  8C018C034C004C02  209  DCW 

0294  CC00CC02CC0 1 CC03  210  DCW 

02A0  2C0 1 2C03AC0OACO2  211  DCW 

02AC  6C006C026C016C03  212  DCW 

02B8  EC01EC031C001C02  213  DCW 

02C4  9C009C029C0 1 9C03  214  DCW 

02D0  5CO15CO3DC00DCO2  215  DCW 

02DC  3C003C023C0 1 3C03  216  DCW 

02E8  BC01BC037C007C02  217  DCW 

02F4  FC00FC02FC01FC03  218  DCW 

0300  0201020382008202  219  DCW 

030C  4200420242014203  220  DCW 

0318  C201C20322002202  221  DCW 

0324  A200A202A201A203  222  DCW 

0330  6201 6203E200E  202  223  DCW 

033C  1200120212011203  224  DCW 

0348  9201920352005202  225  DCW 

0354  D200D202D201D203  226  DCW 

0360  32013203B200B202  227  DCW 

036C  7200720272017203  228  DCW 

0378  F2O1F2030AOOOAO2  229  DCW 

0384  8A008A028AO 1 8A03  230  DCW 

0390  4 AO 1 4 A03CA00CA02  231  DCW 

039C  2A002A022A01 2A03  232  DCW 

03A8  AAO 1 AA036A006A02  233  DCW 

03B4  EA00EA02EA01EA03  234  DCW 

03C0  1 AO 1 1 A039A009A02  235  DCW 

03CC  5A005A025A015A03  236  DCW 

03D8  DA01CA033A003A02  237  DCW 

03E4  BA0OBA02BAO1BA03  238  DCW 

03FO  7A017A03FAOOFA02  239  DCW 

03FC  0600060206010603  240  DCW 

0408  8601860346004602  241  DCW 

0414  C600C602C60 1 C603  242  DCW 

0420  2601 2603A600A602  243  DCW 

042C  6600660266016603  244  DCW 

0438  E601E603 1600 1602  245  DCW 

0444  9600960296019603  246  DCW 

0450  5601560306000602  247  DCW 

045C  3600360236013603  248  DCW 

0468  B601B603 76007602  249  DCW 

0474  F600F602F601F603  250  DCW 

0480  0E010E038E008E02  251  DCW 

048C  4E004E024E014E03  252  DCW 


STATEMENT 

2*60,  2*316,  2*188,  2*444,  2*124,  2*380 
2*252,  2*508,  2*2,  2*258,  2*130,  2*386 
2*66,  2*322,  2*194,  2*450,  2*34,  2*290 
2*162,  2*418,  2*98,  2*354,  2*226,  2*482 
2*18,  2*274,  2*146,  2*402,  2*82,  2*338 
2*210,  2*466,  2*50,  2*306,  2*178,  2*434 
2*114,  2*370,  2*242,  2*498,  2*10,  2*266 
2*138,  2*394,  2*74,  2*330,  2*202,  2*458 
2*42,  2*298,  2*170,  2*426,  2*106,  2*362 
2*234,  2*490,  2*26,  2*282,  2*154,  2*410 
2*90,  2*346,  2*218,  2*474,  2*58,  2*314 
2*186,  2*442,  2*122,  2*378,  2*250,  2*506 
2*6,  2*262,  2*134,  2*390,  2*70,  2*326 
2*198,  2*454,  2*38,  2*294,  2*166,  2*422 
2*102,  2*358,  2*230,  2*486,  2*22,  2*278 
2*150,  2*406,  2*86,  2*342,  2*214,  2*470 
2*54,  2*310,  2*182,  2*438,  2*118,  2*374 
2*246,  2*502,  2*14,  2*270,  2*142,  2*398 
2*78,  2*334,  2*206,  2*462,  2*46,  2*302 
2*174,  2*430,  2*110,  2*366,  2*238,  2*494 
2*30.  2*286,  2*158,  2*414,  2*94,  2*350 
2*2.2,  2*478,  2*62,  2*318,  2*190,  2*446 
2*126,  2*382,  2*254,  2*510,  2*1,  2*257 
2*129,  2*385,  2*65,  2*321,  2*193,  2*449 
2*33,  2*289,  2*161,  2*417,  2*97,  2*353 
2*225,  2*481,  2*17,  2*273,  2*145,  2*401 
2*81,  2*337,  2*209,  2*465,  2*49,  2*305 
2*177,  2*433,  2*113,  2*369,  2*241,  2*497 
2*9,  2*265,  2*137,  2*393,  2*73,  2*329 
2*201,  2*457,  2*41,  2*297,  2*169,  2*425 
2*105,  2*361,  2*233,  2*489,  2*25,  2*281 
2*153,  2*409,  2*89,  2*345,  2*217,  2*473 
2*57,  2*313,  2*185,  2*441,  2*121,  2*377 
2*249,  2*505,  2*5,  2*261,  2*133,  2*389 
2*69,  2*325,  2*197,  2*453,  2*37,  2*293 
2*165,  2*421,  2*101,  2*357,  2*229,  2*485 
2*21,  2*277,  2*149,  2*405,  2*85,  2*341 
2*213,  2*469,  2*53,  2*309,  2*181,  2*437 
2*117,  2*373,  2*245,  2*501,  2*13,  2*269 
2*141,  2*397,  2*77,  2*333,  2*205,  2*461 
2*45,  2*301,  2*173,  2*429,  2*109,  2*365 
2*237,  2*493,  2*29,  2*285,  2*157,  2*413 
2*93,  2*349,  2*221,  2*477,  2*61,  2*317 
2*189,  2*445,  2*125,  2*381,  2*253,  2*509 
2*3,  2*259,  2*131,  2*387,  2*67,  2*323 
2*195,  2*451,  2*35,  2*291,  2*163,  2*419 
2*99,  2*355,  2*227,  2*483,  2*19,  2*275 
2*147,  2*403,  2*83,  2*339,  2*211,  2*467 
2*51,  2*307,  2*179,  2*435,  2*115,  2*371 
2*243,  2*499,  2*11,  2*267,  2*139,  2*395 
2*75,  2*331,  2*203,  2*459,  2*43,  2*299 
2*171,  2*427,  2*107,  2*363,  2*235,  2*491 
2*27,  2*283,  2*155,  2*411,  2*91,  2*347 
2*219,  2*475,  2*59,  2*315,  2*187,  2*443 
2*123,  2*379,  2*251,  2*507,  2*7,  2*263 
2*135,  2*391,  2*71,  2*327,  2*199,  2*455 
2*39,  2*295,  2*167,  2*423,  2*103,  2*359 
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0498 

CE01CE032E002E02 

253 

DCW 

2*231,  2*487,  2*23,  2*279,  2*151,  2*407 

04A4 

AEOOAEO2AEO1AE03 

254 

DCW 

2*87,  2*343,  2*215,  2*471,  2*55,  2*311 

04B0  6E016E03EE00EE02 

255 

OCW 

2*183,  2*439,  2*119,  2*375,  2*247,  2*503 

04BC 

1E001E021E01 1E03 

256 

OCW 

2*15,  2*271,  2*143,  2*399,  2*79,  2*335 

04C8 

9E019E035E005E02 

257 

OCW 

2*207,  2*463,  2*47,  2*303,  2*175,  2*431 

0404 

DE00OE02DE01DE03 

258 

DCW 

2*111,  2*367,  2*239,  2*495,  2*31,  2*287 

04E0 

3E013E03BE00BE02 

259 

OCW 

2*159,  2*415,  2*95,  2*351,  2*223,  2*479 

04EC 

7E007E027E017E03 

260 

DCW 

2*63,  2*319,  2*191,  2*447,  2*127,  2*383 

04F8 

FE01FE03 

261 

262 

DCW 

2*255,  2*511 

04FC 

263 

264 

WR: 

265 

04FC 

FF7FF07FF57FE97F 

266 

OCW 

32767,  32765  ,  32757,  32745  ,  32728  ,  32705 

0508 

A67F867F61 7F377F 

267 

OCW 

32678,  32646,  32609,  32567,  32521,  32469 

0514 

9C7E5F7E107E0570 

268 

OCW 

32412,  32351,  32285,  32213,  32137,  32057 

0520 

E37C887C297CC578 

269 

OCW 

31971,  31880,  31785,  31685,  31580,  31470 

052C 

7C7A057A89790979 

270 

OCW 

31356,  31237,  31113,  30985,  30852,  30714 

0538 

6B 7708764 176A5 75 

271 

DCW 

30571,  30424,  30273,  30117,  29956,  29791 

0544 

8573077354729071 

272 

DCW 

29621,  29447,  29268,  29085,  28898,  28706 

0550 

5E6F966EC960F86C 

273 

OCW 

28510,  28310,  28105,  27896,  27683,  27466 

055C 

606A8B69A668BC67 

274 

OCW 

27245,  27019,  26790,  26556,  26319,  26077 

0568 

E864EE63F162F061 

275 

DCW 

25832,  25582,  25329,  25072,  24811,  24547 

0574 

D75EC75DB35C9C5B 

276 

DCW 

24279,  24007,  23731,  23452,  23170,  22884 

0580 

42581D57F555C954 

277 

DCW 

22594,  22301,  22005,  21705,  21403,  21096 

058C 

3351 FB4F8F4E814D 

278 

DCW 

20787,  20475,  20159,  19841,  19519,  19195 

0598 

B44969481C47CD45 

279 

DCW 

18868,  18537,  18204,  17869,  17530,  17189 

05A4  CE41 73401 73FB83D  280  DCW  16846,  16499,  16151,  15800,  15446,  15090 

05B0  8C392438BA364D35  281  OCW  14732,  14372,  14010,  13645,  13279,  12910 


05BC 

FB30872F 1 1 2E992C 

282 

DCW 

12539, 

12167,  11793,  11417,  11039,  10659 

05C8 

2628A8262825A623 

283 

DCW 

10278, 

9896,  9512,  9126 

8739, 

8351 

0504 

1A1 F931D0B1C821A 

284 

DCW 

7962,  7571,  7179,  6786, 

6393,  5998 

05E0 

E2155514C8123A1 1 

285 

OCW 

5602,  5205,  4808,  4410, 

4011,  3612 

05EC 

8C0CF80A6A09D907 

286 

DCW 

3212,  2811,  2410,  2009, 

1608,  1206 

05  F  8 

2403920100006EFE 

287 

DCW 

804,  402,  0,  -402,  -804 

-1206 

0604 

B8F927F896F605F5 

288 

OCW 

•1608, 

-2009,  - 

2410,  -2811,  -3212,  -3612 

0610 

55F0C6EE38EDABEB 

289 

OCW 

-4011, 

-4410,  - 

4808,  -5205,  -5602,  -5998 

061C 

07E77EE5F5E360E2 

290 

DCW 

•6393, 

-6786,  - 

7179,  -7571,  -7962,  -8351 

0628 

DDDD5ADCD8DA5809 

291 

DCW 

-8739, 

-9126,  - 

9512,  -9896,  -10278,  -10659 

0634 

E1D467D3EFD 17900 

292 

DCW 

-11039, 

-11417, 

-11793, 

-12167, 

-12539, 

-12910 

0640 

21CCB3CA46C9DCC7 

293 

DCW 

-13279, 

-13645, 

-14010, 

-14372, 

•14732, 

15090 

064C 

AAC348C2E9C080BF 

294 

OCW 

-15446, 

-15800, 

-16151, 

-16499, 

-16846, 

17189 

0658 

86B033BAE4B897B7 

295 

DCW 

-17530, 

-17869, 

-18204, 

-18537, 

-18868, 

19195 

0664 

C1B37FB241B105B0 

296 

DCW 

- ’9519, 

-19841, 

-20159, 

-20475, 

-20787, 

21096 

0670 

65AC37AB0BAAE3A8 

297 

DCW 

-21403, 

-21705, 

-22005, 

-22301, 

-22594, 

22884 

067C 

7EA564A440A339A2 

298 

DCW 

-23170, 

-23452, 

-23731, 

-24007, 

-24279, 

24547 

0688 

159F109E0F9D129C 

299 

DCW 

•24811, 

-25072, 

-25329, 

-25582, 

-25832, 

26077 

0694 

3 1 9944985A977596 

300 

DCW 

-26319, 

-26556, 

-26790, 

-27019, 

-27245, 

27466 

06AO  DD93089337926A91 

301 

DCW 

-27683, 

-27896, 

-28105, 

-28310, 

-28510, 

28706 

06  AC 

1E8F638EAC80F98C 

302 

DCW 

-28898, 

-29085, 

-29268, 

-29447, 

-29621, 

29791 

0688 

FC8A5B8ABF892389 

303 

DCW 

-29956, 

-30117, 

-30273, 

-30424, 

-30571, 

30714 

06C4 

7C87F7867786FB85 

304 

DCW 

-30852, 

-30985, 

-31113, 

-31237, 

-31356, 

31470 

0600 

A4843B8407837883 

305 

DCW 

-31580, 

-31685, 

-31785, 

-31880, 

-31971, 

32057 

060  C 

77822B82E381A181 

306 

DCW 

-32137, 

-32213, 

-32285, 

-32351, 

-32412, 

32469 

06E8 

F  780C9809F807A80 

307 

DCW 

-32521, 

-32567, 

-32609, 

-32646, 

-32678, 

32705 

06F4 

2880 1 7800B800380 

308 

DCW 

-32728, 

-32745, 

-32757, 

-32765, 

•32767, 

32765 

0700 

0B801 78028803 F 80 

309 

DCW 

-32757, 

-32745, 

-32728, 

-32705, 

-32678, 

32646 
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070C 

9F80C980F7802B81 

310 

DCW 

-32607,  -32567, 

-32521, 

-32469, 

-32412,  -32351 

0718 

E381 2B827782C782 

311 

DCW 

-32285,  -32213, 

-32137, 

-32057, 

-31971,  -31880 

0724 

D7833B84A484 1 285 

312 

DCW 

-31785,  -31685, 

-31580, 

-31470, 

-31356,  -31237 

0730 

7786F7867C870688 

313 

DCW 

-31113,  -30985, 

-30852, 

-30714, 

-30571,  -30424 

073C 

BF895B8AFC8AA18B 

314 

DCW 

-30273,  -30117, 

-29956, 

-29791, 

•29621,  -29447 

0748 

AC80638E 1 E8FDE8F 

315 

DCW 

-29268,  -29085, 

-28898, 

-28706, 

-28510,  -28310 

0754 

379208930D93B694 

316 

DCW 

-28105,  -27896, 

-27683, 

-27466, 

-27245,  -27019 

0760 

5 A9744983 1 99239A 

317 

DCW 

-26790,  -26556, 

-26319, 

-26077, 

-25832,  -25582 

076C 

0F9D109E159F1DA0 

318 

DCW 

-25329,  -25072, 

-24811, 

-24547, 

-24279,  -24007 

0778 

4DA364A47EA59CA6 

319 

DCW 

-23731,  -23452, 

-23170, 

-22884, 

-22594,  -22301 

0784 

0BAA37AB65AC98AD 

320 

DC.i 

-22005,  -21705, 

-21403, 

-21096, 

-20787,  -20475 

0790 

41B17FB2C1B305B5 

321 

DCW 

-20159,  -19841, 

-19519, 

-19195, 

-18868,  -18537 

079C 

E4B833BA868BDBBC 

322 

DCW 

-18204,  -17869, 

-17530, 

-17189, 

-16846,  -16499 

07A8 

E9C048C2AAC30EC5 

323 

DCW 

-16151,  -15800, 

-15446, 

-15090, 

-14732,  -14372 

0784 

46C9B3CA2 1 CC92C0 

324 

DCW 

-14010,  -13645, 

-13279, 

-12910, 

-12539,  -12167 

07C0 

EFD167D3E1D45DD6 

325 

DCW 

-11793,  -11417, 

-11039, 

-10659, 

-10278,  -9896 

07CC  D8DA5ADCDDDD61DF 

326 

DCW 

-9512,  -9126,  - 

8739,  -8351,  -7962,  -7571 

07D8 

F5E37EE507E792E8 

327 

DCW 

-7179,  -6786,  - 

6393,  -5998,  -5602,  -5205 

07E4 

38E0C6EE55F0E4F1 

328 

DCW 

-4808,  -4410,  - 

4011,  -3612,  -3212,  -2811 

07F0 

96F627F8B8F94AFB 

329 

DCW 

-2410,  -2009,  - 

1608,  -1206,  -804 

,  -402 

07FC 

00009201 2403B604 

330 

DCW 

0,  402,  804,  1206,  1608, 

2009 

0808  6A09F  BO A8C0C 1 COE 

331 

DCW 

2410,  2811,  3212,  3612, 

4011,  4410 

0814 

C8125514E2156E17 

332 

DCW 

4808,  5205,  5602,  5998, 

6393,  6786 

0820 

0B1C931D1A1F9F20 

333 

DCW 

7179,  7571,  7962,  8351, 

8739,  9126 

082C 

2825A8262628A329 

334 

OCW 

9512,  9896,  10278,  10659,  11039, 

11417 

0838  1 1 2E872F  FB306E32  335  DCW  11793,  12167,  12539,  12910,  13279,  13645 

0844  BA3624388C39F23A  336  DCW  14010,  14372,  14732,  15090,  15446,  15800 

0850  173F7340CE4 12543  337  DCW  16151,  16499,  16846,  17189,  17530,  17869 

085C  1C476948B449FB4A  338  DCW  18204,  18537,  18868,  19195,  19519,  19841 

0868  BF4EFB4F335 16852  339  DCW  20159,  20475,  20787,  21096,  21403,  21705 

0874  F5551D5742536459  340  DCW  22005,  22301,  22594,  22884,  23170,  23452 

0880  B35CC750D75EE35F  341  DCW  23731,  24007,  24279,  24547,  24811,  25072 

088C  F162EE63E864DD65  342  DCW  25329,  25582,  25832,  26077,  26319,  26556 

0898  A6688B696D6A4A6B  343  DCW  26790,  27019,  27245,  27466,  27683,  27896 

08A4  C960966E5E6F2270  344  DCW  28105,  28310,  28510,  28706,  28898,  29085 

08B0  54720773B5735F74  345  DCW  29268,  29447,  29621,  29791,  29956,  30117 

08BC  4176D8766B77FA77  346  DCW  30273,  30424,  30571,  30714,  30852,  30985 

08C8  8979057A7C7AEE7A  347  DCW  31113,  31237,  31356,  31470,  31580,  31685 

0804  297C887CE37C397D  348  DCW  31785,  31880,  31971,  32057,  32137,  32213 

08E0  1D7E5F7E9C7ED57E  349  DCW  32285,  32351,  32412,  32469,  32521,  32567 

08EC  617F867FA67FC17F  350  DCW  32609,  32646,  32678,  32705,  32728,  32745 

08F8  F57FFD7FFF7F  351  DCW  32757,  32765,  32767 

352 

353 

08FE  354  WI : 

355 

08FE  00006EFEDCFC4AFB  356  DCW  0,  -402,  -804,  -1206,  -1608,  -2009 

090A  96F605F574F3E4F1  357  DCW  -2410,  -2811,  -3212,  -3612,  -4011,  -4410 

0916  38EDABEB1EEA92E8  358  DCW  -4808,  -5205,  -5602,  -5998,  -6393,  -6786 

0922  F5E360E2E6E061DF  359  DCW  -7179,  -7571,  -7962,  -8351,  -8739,  -9126 


092E  D80A5809DA075DD6 

360 

DCW 

-9512,  -9896,  - 

10278,  - 

10659,  - 

11039, 

11417 

093A  EFD1790005CF92CD 

361 

DCW 

-11793,  -12167, 

-12539, 

-12910, 

•13279, 

-13645 

0946  46C9DCC774C60EC5 

362 

DCW 

-14010,  -14372, 

-14732, 

-15090, 

-15446, 

-15800 

0952  E9C080BF32BEDBBC 

363 

DCW 

-16151,  -16499, 

-16846, 

-17189, 

•17530, 

-17869 

095E  E4B897B74CB605B5 

364 

DCW 

-18204,  -18537, 

-18868, 

-19195, 

-19519, 

-19841 

096A  4 1 B 1 05B0CDAE98AD 

365 

DCW 

-20159,  -20475, 

-20787, 

-21096, 

-21403, 

-21705 

0976  0BAAE3A8BEA79CA6 

366 

DCW 

-22005,  -22301, 

-22594, 

-22884, 

-23170, 

-23452 
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0982  4DA339A229A1 1DA0  367  PCW 

098E  0F9D129C189B239A  368  DCW 

099A  5A9775969395B694  369  OCW 

09A6  37926A91A290DE8F  370  OCW 

09B2  AC80F98C4B8CA18B  371  OCW 

09BE  BF 89288995 880688  372  DCW 

09CA  7786FB858485 1 285  373  DCW 

09D6  078378831D83C782  374  DCW 

09E2  E 381 A181 6481 2B81  375  DCW 

09EE  9F 807A805 A803 F80  376  DCW 

09FA  0B8003800 1800380  377  DCW 

OA06  28803 F805A807A80  378  DCW 

0A12  F7802B816481A181  379  DCW 

0A1E  7782C782 1 D837883  380  DCW 

0A2A  A48412858485FB85  381  DCW 

0A36  7C87068895 882889  382  DCW 

OA42  FC8AA18B4B8CF98C  383  DCW 

OA4E  1E8FDE8FA2906A91  384  DCW 

OA5A  DD93B69493957596  385  OCW 

0A66  3 1 99239A 1 89B 1 29C  386  DCW 

0A72  159F1DA029A139A2  387  DCW 

0A7E  7EA59CA68EA7E3A8  388  DCW 

0A8A  65AC98ADCDAE05BO  389  DCW 

0A96  C1B305B54CB697B7  390  DCW 

0AA2  86BBDBBC32BE8DBF  391  DCW 

OAAE  AAC30EC574C60CC7  392  OCW 

OABA  21CC92CDO5CF79D0  393  DCW 

0AC6  E 1D45DD6DAD758D9  394  DCW 

0AD2  DDDD61DFE6E06DE2  395  DCW 

OADE  07E792E81EEAABEB  396  DCW 

OAEA  55F0E4F174F305F5  397  OCW 

0AF6  B8F94AFBDCFC6EFE  398  DCW 

0B02  2403B6044806D907  399  DCW 

OBOE  8C0C1C0EAB0F3A11  400  DCW 

OBI A  E2156E17F918821A  401  DCW 

0B26  1 A1 F9F202322A623  402  DCW 

0B32  2628A3291 F2B992C  403  DCW 

0B3E  FB306E320  F334035  404  OCW 

0B4A  8C39F 23A563CB83D  405  DCW 

0B56  CE4 1 25437A44CD45  406  DCW 

0B62  B449FB4A3F4C814D  407  OCW 

0B6E  335 1 68529B53C954  408  DCW 

0B7A  A25B6459825A9C5B  409  DCW 

0B86  075EE35FEB60F061  410  OCW 

0B92  E864DD65CF66BC67  411  DCW 

0B9E  606A4A6B236CF86C  412  OCW 

OBAA  5E6F2270E2709D71  413  DCW 

0BB6  B5735F740475A575  414  DCW 

0BC2  6B77FA7784 780979  415  DCW 

OBCE  7C7AEE7A5C7BC57B  416  DCW 

OBDA  E37C397D897DD57D  417  DCW 

0BE6  9C7ED57E097F377F  418  DCW 

0BF2  A67FC17FD87FE97F  419  DCW 
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-23731,  -24007, 

-24279, 

-24547, 

-24811, 

-25072 

-25329,  -25582, 

-25832, 

-26077, 

-26319, 

-26556 

-26790,  -27019, 

-27245, 

-27456, 

-27683, 

-27896 

-28105,  -28310, 

-28510, 

-28706, 

-28898, 

-29085 

-29268,  -29447, 

-29621, 

-29791, 

•29956, 

-30117 

-30273,  -30424, 

-30571, 

-30714, 

-30852, 

-30985 

-31113,  -31237, 

-31356, 

-31470, 

-31580, 

-31685 

-31785,  -31880, 

-31971, 

-32057, 

-32137, 

-32213 

•32285,  -32351, 

-32412, 

-32469, 

-32521, 

-32567 

-32609,  -32646, 

-32678, 

-32705, 

-32728, 

-32745 

■32757,  -32765, 

-32767, 

-32765, 

-32757, 

-32745 

-32728,  -32705, 

-32678, 

-32646, 

-32609, 

-32567 

-32521,  -32469, 

-32412, 

-32351, 

-32285, 

-32213 

-32137,  -32057, 

-31971, 

-31880, 

-31785, 

-31685 

-31580,  -31470, 

-31356, 

-31237, 

-31113, 

-30985 

-30852,  -30714, 

-30571, 

-30424, 

-30273, 

-30117 

-29956,  -29791, 

-29621, 

-29447, 

-29268, 

-29085 

-28898,  -28706, 

-28510, 

23310, 

-28105, 

-27896 

-27683,  -27466, 

-27245, 

•27019, 

-26790, 

-26556 

-26319,  -26077, 

•25832, 

-25582, 

-25329, 

-25072 

-24811,  -24547, 

-24279, 

-24007, 

-23731, 

-23452 

-23170,  -22884, 

-22594, 

-22301, 

-22005, 

-21705 

-21403,  -21096, 

-20787, 

-20475, 

-20159, 

-19841 

-19519,  -19195, 

-18868, 

-18537, 

-18204, 

-17869 

-17530,  -17189, 

-16846, 

-16499, 

-16151, 

-15800 

-15446,  -15090, 

-14732, 

-14372, 

-14010, 

-13645 

-13279,  -12910, 

-12539, 

-12167, 

-11793, 

-11417 

-11039,  -10659, 

-10278, 

-9896, 

9512,  - 

9126 

-8739,  -8351,  -7962,  -7571,  -7179,  -6786 
-6393,  -5998,  -5602,  -5205,  -4808,  -4410 
-4011,  -3612,  -3212,  2811,  -2410,  -2009 

-1608,  -1206,  -804,  -402,  0,  402 
804,  1206,  1608,  2009,  2410,  2811 
3212,  3612,  4011,  4410,  4808,  5205 
5602,  5998,  6393,  6786,  7179,  7571 
7962,  8351,  8739,  9126,  9512,  9896 
10278,  10659,  11039,  11417,  11793,  12167 
12539,  12910,  13279,  13645,  14010,  14372 
14732,  15090,  15446,  15800,  16151,  16499 
16846,  17189,  17530,  17869,  18204,  18537 
18868,  19195,  19519,  19841,  20159,  20475 
20787,  21096,  21403,  21705,  22005,  22301 
22594,  22884,  23170,  23452,  23731,  24007 
24279,  24547,  24811,  25072,  25329,  25582 
25832,  26077,  26319,  26556,  26790,  27019 
27245,  27466,  27683,  27896,  28105,  28310 
28510,  28706,  28898,  29085,  29268,  29447 
29621,  29791,  29956,  30117,  30273,  30424 
30571,  30714,  30852,  30985,  31113,  31237 
31356,  31470,  31580,  31685,  31785,  31880 
31971,  32057,  32137,  32213,  32285,  32351 
32412,  32469,  32521,  32567,  32609,  32646 
32678,  32705,  32728,  32745,  32757,  32765 


OBFE  FF7FF07FF57FE97F  420  DCW  32767,  32765,  32757,  32745,  32728,  32705 


OCOA  A67F867F617F377F  421  DCW  32678,  32646,  32609,  32567,  32521,  32469 


0C16  9C7E5F7E1D7ED57D  422  DCW  32412,  32351,  32285,  32213,  32137,  32057 


0C22  E37C887C297CC578  423  OCW  31971,  31880,  31785,  31685,  31580,  31470 
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0C2E  7C 7A05 7A89790979  424  OCW 

0C3A  6B77D87641 76A575  425  DCW 

0C46  B5 73077354729071  426  DCW 

0C52  5E6F966EC960F86C  427  DCW 

0C5E  6D6A8B69A6688C67  428  DCW 

0C6A  E864EE63F162F061  429  DCW 

0C76  D75EC75DB35C9C5B  430  DCW 

0C82  42581057F555C954  431  OCW 

0C8E  3351FB4FBF4E814D  432  DCW 

0C9A  B34969481C47CD45  433  DCW 

0CA6  CE4 173401 73FB83D  434  DCW 

0CB2  8C392438BA364D35  435  DCW 

OCBE  FB30872F 1 1 2E992C  436  DCW 

OCCA  2628A8262825A623  437  DCW 

0C06  1A1F931D0B1C821A  438  DCW 

0CE2  E2155514C8123A11  439  DCW 

OCEE  8C0C F B0A6A09D907  440  DCW 

OCFA  240392010000  441  DCW 

442 

ODOO  443  END 
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31356,  31237,  31113,  30985,  30852,  30714 
30571,  30424,  30273,  30117,  29956,  29791 
29621,  29447,  29268,  29085,  28898,  28706 
28510,  28310,  28105,  27896,  27683,  27466 
27245,  27C19,  26790,  26556,  26319,  26077 
25832,  25582,  25329,  25072,  24811,  24547 
24279,  24007,  23731,  23452,  23170,  22884 
22594,  22301,  22005,  21705,  21403,  21096 
20787,  20475,  20159,  19841,  19519,  19195 
18867,  18537,  18204,  17869,  17530,  17189 
16846,  16499,  16151,  15800,  15446,  15090 
14732,  14372,  14010,  13645,  13279,  12910 
12539,  12167,  11793,  11417,  11039,  10659 
10278,  9896,  9512,  9126,  8739,  8351 
7962  ,  7571,  7179  ,  6786  ,  6393  ,  5998 
5602,  5205,  4808,  4410,  4011,  3612 
3212,  2811,  2410,  2009,  1608,  1206 
804,  402,  0 
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SYMBOL  TABLE  LISTING 


NAME  VALUE  ATTRIBUTES 

BREV .  OOFCH  CODE  REL  WORD 

ENDL .  00F1H  COOE  REL  ENTRY 

ERR1 .  00B1H  COOE  REL  ENTRY 

ERR2 .  00B5H  COOE  REL  ENTRY 

ERROR .  .  REG  EXTERNAL 

FFT  FO .  .  MODULE  STACKS1ZE<6) 

FFT  FOR .  OOOOH  COOE  REL  PUBLIC  ENTRY 

GM.~ .  0038H  COOE  REL  ENTRY 

GR2  .  006CH  COOE  REL  ENTRY 

GW .  002AH  COOE  REL  ENTRY 

GX .  0087H  COOE  REL  ENTRY 

GX2  .  0075H  COOE  REL  ENTRY 

IK .  0092H  COOE  REL  ENTRY 

IN  CNT .  0042H  OVERLAY  ABS  WORD 

IN  LOOP .  0018H  COOE  REL  ENTRY 

KN2  .  0048H  OVERLAY  ABS  WORD 

KPTR .  0046H  OVERLAY  ABS  WORD 

LOOP_CNT .  0051H  OVERLAY  ABS  BYTE 

MIO  LOOP .  0016H  COOE  REL  ENTRY 

N_SU0_K  .  004AH  OVERLAY  ABS  WORD 

ND1V2  .  0044H  OVERLAY  ABS  WORD 

OUTLOOP .  OOODH  CODE  REL  ENTRY 

PWR .  0040H  OVERLAY  ABS  WORD 

RK .  004CH  OVERLAY  ABS  WORD 

RNK .  004EH  OVERLAY  ABS  WORD 

SHFT  CNT .  0050H  OVERLAY  ABS  BYTE 

TMPI .  0028H  OVERLAY  ABS  LONG 

TMPI1  .  0030H  OVERLAY  ABS  LONG 

TMPR .  0024H  OVERLAY  ABS  LONG 

TMPR1  .  002CH  OVERLAY  ABS  LONG 

UN  LOOP .  OOBFH  COOE  REL  ENTRY 

UNWEAVE .  00B9H  COOE  REL  ENTRY 

WI .  08FEH  COOE  REL  WORD 

WIP .  003EH  OVERLAY  ABS  WORD 

WR .  04FCH  COOE  REL  WORD 

WRP .  003CH  OVERLAY  ABS  WORD 

XI MAG .  .  DATA  EXTERNAL 

XITMP .  0038H  OVERLAY  ABS  LONG 

XREAL .  .  DATA  EXTERNAL 

XRTMP .  0034 H  OVERLAY  ABS  LONG 


ASSEMBLY  COMPLETED,  NO  ERROR(S)  FOUNO. 
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DOS  3.30  (038-N)  MCS-96  MACRO  ASSEMBLER,  Vi. 2 

SOURCE  FILE:  EJNT.A96 
OBJECT  FILE:  EJNT.OBJ 

CONTROLS  SPECIFIED  IN  INVOCATION  COMMAND:  <none> 


ERR  LOC  OBJECT 


SOURCE  STATEMENT 
$PAGELENGTH<51) 

ST1TLEC  ENA8LE  80C196  GLOBAL  INTERRUPTS") 

E_INT  MOOULE  STACKS12EI12) 

PUBLIC  ENAB  I  NT 


LINE  SOURCE  STATEH 

1  $PAGELENGTH<51 ) 

2  STITLEC  ENA8LE  80 

3 

4  E_INT  MC 

5 

6 

7  ;  E_1NT.A96 

8  ; 

9  ;  Version  1. 

10  ; 

11  ;  Jeffrey  C. 

12  ; 

13  ;  GMS  Engine 

14  ;  8940-0  Rou 

15  ;  Columbia, 

16  ; 

17  ; 

18  ;  ================ 

19  ; 

20  ;  COMMON  DEF 

21 

22  SINCLUDE(8096. INC) 


Version  1.0 


Jeffrey  C.  Sigl 


July  26,  1989 


GMS  Engineering  Corporation 
8940-0  Route  108 
Columbia,  Maryland  21045 


COMMON  DEFINITIONS 


22  SINCLUDE(8096. INC)  ;80C196  REGISTER  DEFINITIONS 

=  1  2  3 

*1  24  ; 

=1  25  ;  8096. INC  -  DEFINITION  OF  SYMBOLIC  NAMES  FOR  THE  I/O  REGISTERS  OF  THE 

*1  26  ;  8096  AND  THE  80C196 

=1  27  ;  (C)  INTEL  CORPORATION  1983 

=  1  2  8 

.a****************************************************************************** 


=1 

29 

I 

=1 

30 

;/* 

=1 

31 

;  *  8096  SFR's 

=1 

32 

;  V 

0000 

=1 

33 

RO 

EQU 

00H:WORD 

R  ZERO  REGISTER 

0002 

=1 

34 

AOCOMMAND 

ECU 

02H:BYTE 

W 

0002 

=1 

35 

AD  RESULT  LO 

EQU 

02H:BYTE 

R 

0003 

=1 

36 

AD  RESULT  HI 

EOU 

03H.-BYTE 

R 

0003 

=1 

37 

HSI  MODE 

EQU 

03H:BYTE 

W 

0004 

=1 

38 

HSO  TIME 

EQU 

04H :WORD 

U 

0004 

=  1 

39 

HSI_TIME 

EQU 

04 H: WORD 

R 

0006 

=1 

40 

HSO  COMMAND 

EOU 

06H:BYTE 

W 

0006 

=1 

41 

HSI  STATUS 

EOU 

06H.-BYTE 

R 

0007 

=1 

42 

SBUF 

EQU 

07H:BYTE 

R/W 

MCS-96  MACRO  ASSEMBLER 


ENABLE  80C196  GLOBAL  INTERRUPTS 


04/06/90  16:17:04  PAGE  2 


ERR  LOC  OBJECT 

LINE 

SOURCE 

STATEMENT 

0008 

=1 

43 

I NT  MASK 

EOU 

08H:BYTE 

R/W 

0009 

=1 

44 

I NT  PENDING 

EQU 

09H:BYTE 

R/W 

000A 

=1 

45 

WATCHDOG 

EQU 

OAH:BYTE 

W 

WATCHDOG  TIMER 

OOOA 

=  1 

46 

TIMER1 

EQU 

OAH:UORD 

R 

OOOC 

=1 

47 

TIMER2 

EQU 

OCH:WORO 

R 

OOOE 

=  1 

48 

BAUD  RATE 

EQU 

OEH:BYTE 

W 

OOOE 

=1 

49 

IOPORT0 

EQU 

OEH:BYTE 

R 

000  F 

=1 

50 

I0P0RT1 

EQU 

OFH:BYTE 

R/W 

0010 

=  1 

51 

IOPORT2 

EQU 

10H:BYTE 

R/W 

0011 

=1 

52 

SP  CON 

EQU 

11H:BYTE 

W 

0011 

=1 

53 

SP  STAT 

EQU 

11H:BYTE 

R 

0015 

=1 

54 

IOCO 

EQU 

15H:BYTE 

W 

0015 

=1 

55 

IOSO 

EQU 

15H.-BYTE 

R 

0016 

=1 

56 

I0C1 

EQU 

16H:BYTE 

W 

0016 

=1 

57 

I0S1 

EQU 

16H:BYTE 

R 

0017 

=1 

58 

PWM  CONTROL 

EQU 

17H:BYTE 

W 

0018 

=1 

59 

SP 

EQU 

18H.-W0RD 

R/W 

=1 

60 

=1 

61 

;/* 

=1 

62 

;  *  80C196  SFR's 

=1 

63 

;  */ 

OOOB 

=1 

64 

I0C2 

EOU 

0BH:8YTE 

W 

=1 

65 

;TIMER2 

EQU 

0CH:WORD 

R/W 

0012 

=1 

66 

I  PEND 1 

EQU 

12H:BYTE 

R/W 

0013 

=1 

67 

IMASK1 

EQU 

13H:BYTE 

R/W 

0014 

=1 

68 

WSR 

EQU 

14H:BYTE 

R/W 

0017 

=  1 

69 

I0S2 

EQU 

17H:BYTE 

R 

70 

71 

72 

•  ~~ 

t 

73 

;  CODE 

SEGMENT 

74 

0000 

75 

CSEG 

76 

0000 

77 

ENAB  I NT: 

78 

0000  FB 

79 

El 

; ENABLE 

INTRPTS 

0001  FO 

80 

RET 
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SYM80L  TABLE  LISTING 


NAME  VALUE  ATTRIBUTES 

AD_COMMAND .  0002H  NULL  ABS  BYTE 

AD_RESULT_HI .  0003H  NULL  ABS  BYTE 

ADRESULTLO .  0002H  NULL  ABS  BYTE 

BAUORATE .  000EH  NULL  ABS  BYTE 

E_INT  .  .  MODULE  STACKSIZE(12) 

ENAB_INT .  0000 H  CODE  REL  PUBLIC  ENTRY 

MSI  MODE .  0003H  NULL  ABS  BYTE 

HS1STATUS .  0006H  NULL  ABS  BYTE 

HSI_TIME .  0004H  NULL  ABS  WORD 

HS0_C0MMAND .  0006H  NULL  ABS  BYTE 

HSO_T  I  ME .  0004 H  NULL  ABS  WORD 

I  MASK  1 .  0013H  NULL  ABS  BYTE 

INT  MASK .  0008H  NULL  ABS  BYTE 

INT  PENDING .  0009H  NULL  ABS  BYTE 

IOCO .  0015H  NULL  ABS  BYTE 

I0C1 .  0016H  NULL  ABS  BYTE 

I0C2 .  OOOBH  NULL  ABS  BYTE 

I OPORTO  . .  OOOEH  NULL  ABS  BYTE 

I0P0RT1 .  OOOFH  NULL  ABS  BYTE 

I0P0RT2 .  0010H  NULL  ABS  BYTE 

1050  .  001 5H  NULL  ABS  BYTE 

1051  .  0016H  NULL  ABS  BYTE 

1052  .  0017H  NULL  ABS  BYTE 

I  PEND  1 .  001 2H  NULL  ABS  BYTE 

PWMCONTROL .  001 7H  NULL  ABS  BYTE 

RO .  OOOOH  NULL  ABS  WORD 

SBUF .  0007H  NULL  ABS  BYTE 

SP .  0018H  NULL  ABS  WORD 

SP  CON .  001 1H  NULL  ABS  BYTE 

SP_STAT .  001 1H  NULL  ABS  BYTE 

TIMER1 .  OOOAH  NULL  ABS  WORD 

TIMER2 .  OOOCH  NULL  ABS  WORD 

WATCHDOG .  OOOAH  NULL  ABS  BYTE 

WSR .  0014H  NULL  ABS  BYTE 


ASSEMBLY  COMPLETED,  NO  ERROR(S)  FOUND. 
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DOS  3.30  (038-N)  MCS-96  MACRO  ASSEMBLER,  VI. 2 

SOURCE  FILE:  STATUS.A96 
OBJECT  FILE:  STATUS. OBJ 

CONTROLS  SPECIFIED  IN  INVOCATION  COMMAND:  <none> 


status_temp 

rseg 

1  ,-Global  status  register 

end 


SYMBOL  TABLE  LISTING 


NAME  VALUE  ATTRIBUTES 

STATUS  TEMP .  0000H  REG  REL  PUBLIC  BYTE 


ASSEMBLY  COMPLETED,  NO  ERROR(S)  FOUND. 


ERR  LOC  OBJECT 

LINE 

SOURCE  STATEMENT 

1 

public 

0000 

2 

0000 

3 

status_tenp:  OSB 

0001 

4 

HCS-96  MACRO  ASSEMBLER  GETCHAR 
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DOS  3.30  (038-N)  MCS-96  MACRO  ASSEMBLER,  VI. 2 

SOURCE  FILE:  GETCHAR. A96 
08JECT  FILE:  GETCHAR. OBJ 

CONTROLS  SPECIFIED  IN  INVOCATION  COMMAND:  <none> 


ERR  LOC  OBJECT 

LINE 

SOURCE  STATEMENT 

1 

Sdebug 

2 

Snolist 

Include  (8096. inc) 

51 

001C 

52 

tmpO 

equ 

1CH:word 

0006 

53 

RI_pos 

equ 

06H:byte 

OOBF 

54 

RI_mask 

equ 

0BFH:byte 

55 

56 

extrn 

status_temp 

57 

public 

getchar 

58 

0000 

59 

CSEG 

60 

0000  901100 

E 

61 

getchar: 

orb 

status_temp,  SP_STAT 

0003  3600FA 

E 

62 

jbc  status_temp,  RI_pos,  getchar 

0006  B0071C 

63 

Idb  tmpO,  sbuf 

0009  71BF00 

E 

64 

andb  status_tenf>,  #RImask 

000C  FO 

65 

ret 
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SYMBOL  TABLE  LISTING 


NAME  VALUE  ATTRIBUTES 

AD_COMMAND .  0002H  NULL  ABS  BYTE 

AD_RESULT_HI .  0003H  NULL  ABS  BYTE 

AD_RESULTJ.O .  0002H  NULL  ABS  BYTE 

BAUD_RATE .  000EH  NULL  ABS  BYTE 

GETCHAR .  OOOOH  COOE  REL  PUBLIC  ENTRY 

HSI  MOOE .  0003H  NULL  A8S  BYTE 

HSI  STATUS .  0006H  NULL  ABS  BYTE 

HSI_TIME .  0004 H  NULL  ABS  WORD 

HSO  COMHAND .  0006H  NULL  ABS  BYTE 

HSO  TIME .  0004 H  NULL  ABS  WORD 

I  MASK  1 .  0013H  NULL  ABS  BYTE 

INT  MASK .  0008H  NULL  ABS  BYTE 

INT  PENDING .  0009H  NULL  ABS  BYTE 

IOCO .  0015H  NULL  ABS  BYTE 

I0C1 .  0016H  NULL  ABS  BYTE 

I0C2 .  OOOBH  NULL  ABS  BYTE 

I  OPORTO .  OOOEH  NULL  ABS  BYTE 

I0P0RT1 .  OOOFH  NULL  ABS  BYTE 

I0P0RT2 .  001  OH  NULL  ABS  BYTE 

1050  .  001 5H  NULL  ABS  BYTE 

1051  .  0016H  NULL  ABS  BYTE 

1052  .  0017H  NULL  ABS  BYTE 

I  PEND  1 .  0012H  NULL  ABS  BYTE 

PWM  CONTROL .  001 7H  NULL  ABS  BYTE 

RO .  OOOOH  NULL  ABS  WORD 

RI_MASK .  OOBFH  NULL  ABS  BYTE 

RI_POS .  0006H  NULL  ABS  BYTE 

SBl'F .  0007H  NULL  ABS  BYTE 

.  0018H  NULL  ABS  WORD 

.  001 1H  NULL  ABS  BYTE 

SF_  .  001 1H  NULL  ABS  BYTE 

STATU  MP .  .  NULL  EXTERNAL 

TIMER1 .  OOOAH  NULL  ABS  WORD 

TIMER2 .  OOOCH  NULL  ABS  WORD 

TMPO .  001CH  NULL  ABS  WORD 

WATCHDOG .  OOOAH  NULL  ABS  BYTE 

WSR .  0014H  NULL  ABS  BYTE 

ASSEMBLY  COMPLETED,  NO  ERROR(S)  FOUND. 


MCS-96  MACRO  ASSEMBLER  PUTCHAR 
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SOURCE  FILE:  PUTCHAR. A96 
OBJECT  FILE:  PUTCHAR. OBJ 

CONTROLS  SPECIFIED  IN  INVOCATION  COMMAND:  <none> 


LOC  OBJECT 

LINE 

SOURCE  STATEMENT 

1 

(debug 

2 

$nol 1st 

include  (8096. inc) 

51 

0005 

52 

TI_pos 

equ 

05H:I 

00OF 

53 

T!_mask 

equ 

OOFH 

54 

55 

extrn 

status_temp 

56 

public 

putchar 

57 

0000 

58 

CSEG 

59 

0000  901100 

E 

60 

putchar: 

orb 

0003  3500FA 

E 

61 

j  be 

0006  B3 180207 

62 

(db 

000A  71DF00 

E 

63 

andb 

0000  FO 

64 

ret 

OOOE 

65 

end 

statustemp,  SPSTAT 

status_temp,  TI_pos,  putchar 
sbuf,  2fspJ 
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SYMBOL  TABLE  LISTING 


NAME  VALUE  ATTRIBUTES 

AD_COMMAND .  0002H  NULL  ABS  BYTE 

AD_RESULT_H! .  0003H  NULL  ABS  BYTE 

ADRESULTLO .  0002H  NULL  ABS  BYTE 

BAUD_RATE .  000EH  NULL  ABS  BYTE 

HSI  MODE .  0003H  NULL  ABS  BYTE 

HSI  STATUS .  0006H  NULL  ABS  BYTE 

HSI  TIME .  0004 H  NULL  ABS  WORD 

HSOCOMMAND .  0006H  NULL  ABS  BYTE 

HSO_TIME .  0004 H  NULL  ABS  WORD 

IMASK1 .  0013H  NULL  ABS  BYTE 

1  NT  MASK .  0008H  NULL  ABS  BYTE 

1NTPEND1NG .  0009H  NULL  ABS  BYTE 

IOCO .  0015H  NULL  ABS  BYTE 

I0C1 .  0016H  NULL  ABS  BYTE 

I0C2 .  OOOBH  NULL  ABS  BYTE 

l OPORTO .  OOOEH  NULL  ABS  BYTE 

I0P0RT1 .  OOOFH  NULL  ABS  BYTE 

I0P0RT2 .  001  OH  NULL  ABS  BYTE 

1050  .  0015H  NULL  ABS  BYTE 

1051  .  0016H  NULL  ABS  BYTE 

1052  .  0017H  NULL  ABS  BYTE 

I  PEND  1 .  001 2H  NULL  ABS  BYTE 

PUTCHAR .  OOOOH  CODE  REL  PUBLIC  ENTRY 

PWMCONTROL .  001 7H  NULL  ABS  BYTE 

RO .  OOOOH  NULL  ABS  WORD 

SBUF .  0007H  NULL  ABS  BYTE 

SP .  0018H  NULL  ABS  WORD 

SP  CON .  001 1H  NULL  ABS  BYTE 

SP  STAT .  001 1H  NULL  ABS  BYTE 

STATUS_TEMP  .  .  NULL  EXTERNAL 

T I  MASK .  CODFH  NULL  ABS  BYTE 

TI_POS .  0005 H  NULL  ABS  BYTE 

TIMER1 .  OOOAH  NULL  ABS  WORD 

TIMER2 .  OOOCH  NULL  ABS  WORD 

WATCHDOG .  OOOAH  NULL  ABS  8YTE 

WSR .  0014H  NULL  ABS  BYTE 


ASSEMBLY  COMPLETED,  NO  ERROR(S)  FOUND. 


APPENDIX  B:  PC  SOFTWARE  LISTING 


GMS  Engineering  Corp. 


EEC  Artifact  Rejection  System 


SHELL 

PAGE  1 
03-22-90 
17:43:47 

Line#  Source  Line  Hierosoft  FORTRAN  Optimizing  Compiler  Version  4.00 

1  tdeclare 

2  c  shell. for 

3 

4  c  Driver  for  the  EEG  artifact  correction  system. 

5 

6 

7  c  Created: December  19,  1989 

8  c  Last  Update:  February  14,  1989 

9 

10  c  Steven  M.  Falk 

11  c  Jeffrey  C.  Sigl 

12  c  GMS  Engineering  Corporation 

13  c  8940 -D  Route  108 

14  c  Columbia,  HD  21045 

15  c  (301)  995-0508 

16 

17  program  shel l 

18 

19  c  Data  Structures 

20 

21  character*6  word 

22  character*1  chr(6),  cls(4),  capps,  cappi,  cr,  cappr,  cappz 

23  character^  bell,  resp 

24  integer  i,  j 

25  integer*4  tmpvb,  tmpvbb,  tmpva 

26  integer  step,  iresp,  ichanl,  ichann 

27  real*4  dlv,  d2v,  dlh,  d2h 

28  real *4  d3v,  d4v,  d3h,  d4h 

29  real*4  vmag,  hmag,  gain,  ccmax 

30  real*4  tmpvbl(ll) 

31 

32  c  Functions 

33 

34  integer  ichar 

35  integer*4  int4 

36 

37  c  Data  Relations 

38 

39  equivalence  (word.chr) 

40 

41  c  Data  Initialization 

42 

43  data  word  /6H  / 

44  data  capps,  cappi,  cappr,  cappz  /'S',  'I',  'R‘,  'A1/ 

45  data  ichanl,  ichann  /I,  0/ 

46  cls(1>  =  8#33 

47  cls(2)  =  8#133 

48  cls(3)  =  8#62 

49  cls(4)  =  8#112 

50  bell  =  8#7 

51 

52  step  =  0 

53 

54  dlv  =  0. 

55  d2v  =  0. 

56  dlh  =  0. 
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.  ine# 

Source 

Line  Microsoft  FORTRAN  Optimizing  Compiler  Version  4.00 

57 

d2h  =  0. 

58 

d3v  =  0. 

59 

d4v  =  0. 

60 

d3h  =  0. 

61 

d4h  =  0. 

62 

63 

64 

c  Clear  the  screen 

65 

66 

2 

write(*,3)(els(i), i=1,4> 

67 

3 

formate  *,4*1) 

68 

69 

open(11,f i le='drive' ,status='old' ) 

70 

66 

read(11,67)resp 

71 

67 

format(al) 

72 

if  (  resp  .ne.  cappz  )  then 

73 

goto  66 

74 

end  if 

75 

76 

call  dcalc(cls,d1v,d2v,d1h,d2h,d3v(d4v,d3h,d4h,step,bell,ccniax) 

77 

78 

call  sgm(vmag,0) 

79 

call  sgm(hmag,1) 

80 

call  sgm(gain,2) 

81 

82 

call  ftune(d1v,d2v,d1h,d2h,d3vfd4v,d3h,d4h,vmag,hmag,gain,ccmax) 

83 

84 

close(ll) 

85 

86 

if  (  gain  .eq.  -9999.  )  then 

87 

goto  1357 

88 

endif 

89 

90 

open( 1 2 , f i l e= ' correct ' , status? ■ new 1 > 

91 

92 

tmpvbl(l)  =  vrrtag 

93 

tmpvbl(2)  =  hmag 

94 

tmpvbl(3)  =  gain 

95 

tmpvbl(4)  =  dlv 

96 

tmpvbl(5)  =  d2v 

97 

tmpvbl(6)  =  d3v 

98 

tmpvbl(7)  =  d4v 

99 

tmpvbl(8)  =  dlh 

100 

tmpvbl(9)  =  d2h 

101 

tmpvbl(IO)  =  d3h 

102 

tmpvbl(ll)  =  d4h 

103 

do  152  i=1 , 11 

104 

tmpva  =  int4(tmpvbl(i )) 

105 

trpvb  =  tmpva  *  10 

106 

if  (  tmpva  .It.  0  )  then 

107 

tmpvb  =  Impvb 

108 

tmpvb  =  tmpvb  *  1 

109 

endif 

110 

writef 12, 154) tmpvb 

111 

154 

formatf i6) 

112 

152 

continue 

continue 
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Line#  Source  Line  Microsoft  FORTRAN  Optimizing  Compiler  Version  4.00 

113 

114  close(12) 

115 

116  write(*,1356) 

117  1356  formate  Correction  Matrix  now  computed  for  the  selected  channel.' 

118  +  /'  EARS  Correction  opotion  can  now  be  run  (main  menu  option  C).' 

119 

120  1357  continue 

121 

122  stop 

123  end 


main  Local  Symbols 

Name 

Class 

Type 

Size 

Offset 

D4H . 

.  local 

REALM 

4 

0002 

IRESP  . 

.  local 

INTEGER*4 

4 

0006 

TMPVA  . 

.  local 

INTEGERM 

4 

000a 

WORD . 

.  local 

CHAR*6 

6 

000c 

TMPVB  . 

.  local 

INTEGERS 

4 

OOOe 

CAPPS  . 

.  local 

CHAR*1 

1 

0012 

I . 

.  local 

INTEGERM 

4 

0012 

CAPPI  . 

.  local 

CHAR*1 

1 

0013 

CAPPR  . 

.  local 

CHAR*1 

1 

0014 

CAPPZ  . 

.  local 

CHARM 

1 

0015 

J . 

.  local 

INTEGERM 

4 

0016 

1CHANL . 

.  local 

INTEGERS 

4 

0016 

ICHANN . 

.  local 

INTEGERS 

4 

001a 

D1V . 

.  local 

REALM 

4 

001a 

TMPVBB . 

.  local 

INTEGERM 

4 

001  e 

D2V . 

.  local 

REALM 

4 

0022 

D3V . 

.  local 

REALM 

4 

0026 

04V . 

.  local 

REALM 

4 

002a 

TMPVBL . 

.  local 

REALM 

44 

002e 

CR . 

.  local 

CHARM 

1 

005a 

HMAG . 

.  local 

REALM 

4 

005c 

GAIN . 

.  local 

REALM 

4 

0060 

BELL . 

.  local 

CHARM 

1 

0064 

CLS . 

.  local 

CHARM 

4 

0066 

I CHAR  . 

.  local 

INTEGERM 

4 

006a 

VMAG . 

.  local 

REALM 

4 

006e 

CCMAX  . 

.  local 

REALM 

4 

0072 

RESP . 

.  local 

CHARM 

1 

0076 

STEP . 

.  local 

INTEGERM 

4 

0078 

D1H . 

.  local 

REALM 

4 

007c 

02H . 

.  local 

REALM 

4 

0080 

D3H . 

.  local 

REALM 

4 

0084 

CHR . 

.  equiv 

CHARM 

6 

000c 
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Global  Symbols 


Name  Class  Type  Size  Offset 

DCALC . extern  ***  ***  *** 

FTUNE . extern  ***  ***  *** 

SGM . extern  ***  ***  *** 

main . FSUBRT  ***  ***  0000 

Code  size  =  03a8  (936) 

Data  size  =  006b  (107) 

Bss  size  =  0088  (136) 


No  errors  detected 


0CALC2 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


Sstorage:2 

(declare 

c  dcalc.for 

c 

c  a  program  to  compute  the  D's  for  the  EEG  correction  algorithm 

c 

c  October  12,  1989 

c 

c  Jeffrey  C.  Sigl 

subroutine  dcalc(cls,d1v,d2v,d1h,d2h,d3v,d4v,d3h,d4h, step, bell, 
♦ccmax) 

c  Data  Structures 


characters 

characters 

integer 

integer 

real*4 

real*4 

double  precision 
double  precision 
double  precision 
double  precision 
double  precision 
double  precision 
double  precision 
double  precision 
double  precision 


cls(4) 
bell,  resp 

i,  j,  n,  zline,  ezline,  step,  iresp 
stepl 

dlv,  d2v,  dlh,  d2h 
d3v,  d4v,  d3h,  d4h,  ccmax 
z,  h,  phi,  psi,  pi2,  capl,  raddeg 
betel,  betar,  gamal,  gamar 
theta,  r,  sigma,  eta 
betalv,  betarv,  gamalv,  gamarv 
thetav,  rv,  sigmav,  etav 
betalh,  betarh,  gamalh,  gamarh 
thetah,  rh,  sigmah,  etah 
denomv,  denomh,  temp,  temp2,  terml,  term2 
delta,  jeff,  negflag,  max,  scale 


c  Functions 


integer  ichar 

doi*»le  precision  dacos,  dsin,  dcos,  dsqrt,  dabs 


c  Data  Initialization 

pi2  =  1.5707963268000 
raddeg  =  57.295779513100 

n  =  1 
z  =  111. 
h  =  70. 
capl  =  8. 
betalv  =  33. 
betarv  =  94. 
betalh  =  49. 
betarh  =  144. 
zline  =  0 
beta l  =  208. 
betar  =211. 
ccmax  =0.05 
ezline  =  0 


Line# 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 
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c  Clear  the  screen 

2  write(*,3)(cls(i),i=1,4) 

3  formate  ',4a1) 

c  Write  the  Main  Menu 

wri te(*,50)z,h,capl,betalv,betarv,betalh,betarh,zline,betal , 
+betar,ccmax 

50  format(//20x, 1  GMS  Engineering  EEG-EOG  Artifact  Removal'//, 
+20x,'  Parameter  Menu'////, 

♦  '  1 '  ,5x, 'Distance  between  the  eyes  (mm) . ', 

+f9.2,/, 

♦'  2' ,5x, 'Distance  from  the  origin  to  the  stim  electrode  (mm)’, 

+f9.2/, 

+  '  3'  ,5x,  'Corneo- retina l  distance  (urn) . ', 

+f9.2,/, 

+'  4' ,5x, 'Distance  from  Left  Eye  to  VU-EOG  electrode  (mm)....', 

+f9.2,/, 

+  '  5' ,5x, 'Distance  from  Right  Eye  to  VU-EOG  electrode  (nm)...', 

+f9.2,/ , 

+'  6' ,5x, 'Distance  from  Left  Eye  to  H-EOG  electrode  (nm) . ', 

+f9.2,/, 

+  '  7' ,5x, 'Distance  from  Right  Eye  to  H-EOG  electrode  (inn) _ ', 

+f9.2,/ , 

+'  8', 5x, 'H-EOG  electrode  above(0)/below(1)  the  "Z"  line?....  ', 

+i5./. 

*'  9', 5x, 'Distance  from  Left  Eye  to  EEC  electrode  (mm) . ', 

+f9.2./, 

+  '  10' ,5x, 'Distance  from  Right  Eye  to  EEG  electrode  (mm) . ’, 

♦f9.2,/ , 

♦'  11 ' ,5x, 'Haxiimjn  cross-correlation  function  for  correction..', 
♦f9.2,/ , 

+'  12' ,5x, 'Physical  data  entered;  compute  Correction  Matrix'/, 
*/12x,'  Enter  Response  >'\) 


c  Read  the  response  &  take  appropriate  action 

read(*,'(i5)')iresp 

if  (iresp  .eq.  1)  then 
write(*,110) 
read(*,'(f12.5)')z 
if  (  z  .It.  0.  )  then 
z  =  0. 

call  error(bell) 

endif 

elseif  (iresp  .eq.  2)  then 
write(*,110) 
read(*,'(f12.5)')h 
if  (  h  .It.  0.  )  then 
h  =  0. 

call  error(bell) 


endif 
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1 

|  113 

1 

I  114 

elseif  (iresp  .eq.  3)  then 

115 

write(*, 110) 

1 

.  116 

110  format(/13x, 'Enter  new  value  >'\) 

117 

read(*,'(f12.5)')capl 

1 

|  118 

if  (  capl  .It.  0.  )  then 

119 

capl  =  8. 

120 

call  error(bell) 

1 

121 

endif 

122 

123 

elseif  (iresp  .eq.  4)  then 

124 

write(*,110) 

1  125 

read(*,'(f12.5)')betalv 

!  126 

if  <  betalv  .It.  0.  )  then 

■ 

1  127 

betalv  =  0. 

128 

call  error(bell) 

■ 

»  129 

endif 

130 

1 

1  131 

elseif  (iresp  .eq.  5)  then 

132 

write(*,110) 

1  133 

read(*, '(f 12.5)' )betarv 

1 

134 

if  (  betarv  .It.  0.  )  then 

1 

135 

betarv  =  0. 

136 

call  error(bell) 

137 

endif 

■  138 

1  139 

elseif  (iresp  .eq.  6)  then 

■  140 

write(*,110) 

141 

read(*,'(f12.5)')betalh 

m  142 

if  (  betalh  .It.  0.  )  then 

■  143 

betalh  =  0. 

■  144 

call  error(bell) 

145 

endi  f 

_  146 

H  147 

elseif  (iresp  .eq.  7)  then 

H  148 

write(*,110) 

149 

read(*,*(f12.5)')betarh 

150 

if  (  betarh  .It.  0.  )  then 

■  151 

betarh  =  0. 

■  152 

call  error(bell) 

153 

endif 

154 

■  155 

elseif  (iresp  .eq.  8)  then 

»  156 

write(*,110) 

■  157 

read(*,'(i5)')zline 

158 

if  (  (zline  .ne.  0)  .and.  (zline  .ne.  1))  then 

m  159 

zline  =  0 

■  160 

call  error(bell) 

■  161 

endif 

162 

_  163 

elseif  (iresp  .eq.  9)  then 

■  164 

write(*,110) 

■  165 

read(*,'(f12.5)')betal 

166 

if  (  betal  .It.  0.  )  then 

167 

betel  =  0. 

H  168 

1 

1 

call  error(bell) 
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169 

end  if 

170 

171 

elseif  (iresp  .eq.  10)  then 

172 

write(*,110) 

173 

read(*, ' (f12.5)')betar 

174 

if  (  betar  .It.  0.  )  then 

175 

betar  =  0. 

176 

call  error(bell) 

177 

endif 

178 

179 

elseif  (iresp  .eq.  11)  then 

180 

write(*,110) 

181 

read(*, l(f12.5)')ccmax 

182 

if  (  ccmax  .It.  0.  )  then 

183 

ccmax  =  0. 

184 

call  error(bell) 

185 

endif 

186 

187 

elseif  (iresp  .eq.  12)  then 

188 

if  ((betal. eq.0.). or. (betar. eq.0.))  then 

189 

write(*f 155)bell 

190 

155 

format(//,12x, 'EEG  electrode  distances'. 

191 

+ 

1  must  be  entered! '//,a1,12x. 

192 

+ 

’Type  ENTER  to  continue. . . • ) 

193 

read(*, '(bn,a1)')resp 

194 

goto  2 

195 

endif 

196 

goto  98 

197 

198 

else 

199 

call  error(bell) 

200 

201 

endif 

202 

203 

goto  2 

204 

205 

206 

c  Check 

if  the  VEOG  electrode  is  below  the  stim  electrode; 

207 

c 

if  so,  negate  the  Ov's 

208 

209 

98 

delta  =  dacos(  (z**2  -  betalv**2  +  betarv**2) 

210 

# 

/  (2.*z*betarv)  ) 

211 

if  (  (betarv*dsin(delta))  .It.  h  )  then 

212 

negflag  =  -1.000 

213 

c 

write(*,170) 

214 

c170 

format(/'  The  EOGv  electrode  is  below  the  STIM  electrode 

215 

else 

216 

negflag  =  1.000 

217 

c 

write(*,172) 

218 

cl  72 

formate  The  EOGv  electrode  is  above  the  STIM  electrode. 

219 

endif 

220 

221 

c  Compute  Vertical  EOG  parameters 

222 

223 

c 

write(*,23) 

224 

c23 

format(/ 1  EOGV') 
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225 

226  rv  =  dsqrt<  0.5  *  (  betalv**2  +  betarv**2  -  ((z**2)/2.)  )  ) 

227  c  write(*,200)rv 

228  c200  formate  r  *ff16.8) 

229 

230  temp  =  (rv** 2  ♦  (z/2)**2  -  betarv**2)  /  (rv  *  z) 

231  if  (  temp  .gt.  1.  )  then 

232  temp  =  1.0 

233  elseif  (  temp  .It.  -1.  )  then 

234  temp  =  -1.0 

235  endif 

236  thetav  =  dacos(  temp  ) 

237  c  write(*,210)thetav*raddeg 

238  c210  formate  theta  ' ,f16.8, '  degrees') 

239 

240  gamalv  =  {  rv/betalv)  *  dsinf  thetav  ) 

241  c  write(*,230)gamalv 

242  c230  formate  gamal  ',f16.8) 

243 

244  gamarv  =  (  rv/betarv)  *  dsin(  thetav  ) 

245  c  write(*,240)gamarv 

246  c240  formate  gamar  ’,f16.8) 

247 

248  sigmav  =  dsqrt(  rv**2  +  h**2  -  2. *rv*h*dsin{ thetav)  ) 

249  c  write(*,250)sigmav 

250  c250  formate  sigma  ' ,f16.8) 

251 

252  etav  =  (-rv/sigmav)  *  dcos(  thetav  ) 

253  c  write(*,260)etav 

254  c260  formate  eta  ' ,f16.8) 

255 

256  c  Compute  Horizontal  EOG  parameters 

257 

258  c  write(*,24) 

259  c24  format(/ '  EOGH 1 ) 

260 

261  rh  =  dsqrt(  0.5  *  (  betalh**2  +  betarh**2  -  ((z**2)/2.)  )  ) 

262  c  vrite(*,200)rh 

263 

264  temp  =  (rh**2  +  (z/2)**2  -  betarh**2)  /  (rh  *  z) 

265  if  (  temp  .gt.  1.  )  then 

266  temp  =1.0 

267  elseif  (  temp  .It.  -1.  )  then 

268  temp  =  -1.0 

269  endif 

270  thetah  =  dacos(  temp  ) 

271  if  (  zline  .eq.  1  )  thetah  =  -thetah 

272  c  write(*,210)thetah*raddeg 

273 

274  gamalh  =  (  rh/betalh)  *  dsin(  thetah  ) 

275  c  write(*,230)gamalh 

276 

277  gamarh  =  (  rh/betarh)  *  dsin(  thetah  ) 

278  c  write(*f240)gamarh 

279 

280 


sigmah  =  dsqrt(  rh**2  +  h**2  -  2.*rh*h*dsin(thetah)  ) 
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281 

c 

writef* 

,250)sigmah 

282 

283 

etah  = 

(-rh/sigmah)  *  dcos(  thetah  ) 

284 

c 

write(* 

,260) etah 

285 

286 

287 

c  Compute  EEG  electrode  parameters 

288 

289 

998 

do  1000 

i  =  1,  n 

290 

291 

r  =  dsqrt(0.5*(betal**2  +  betar**2  -  <(z**2>/2.)>) 

292 

c 

write(*,200)r 

293 

294 

temp  =  (  r**2  +  (z/2)**2  •  betar**2  )  /  (r  *  z) 

295 

if  (  temp  .gt.  1.  )  then 

296 

temp  =1.0 

297 

elseif  (  temp  .It.  -1.  )  then 

298 

temp  =  -1.0 

299 

endif 

300 

theta  =  dacos(  temp  ) 

301 

if  (  ezline  .eq.  1  )  theta  =  -theta 

302 

c 

wri te(*,210)theta*raddeg 

303 

304 

gamal  =  (  r/betal)  *  dsin<  theta  ) 

305 

gamar  =  (  r/betar)  *  dsin(  theta  ) 

306 

c 

wr i te( * , 230 ) gama l 

307 

c 

write(*,240)gamar 

308 

309 

sigma  =  dsqrtf  r**2  +  h**2  •  2.*r*h*dsin(theta)  ) 

310 

eta  =  (-r/sigma)  *  dcos(  theta  > 

311 

c 

write(*,250)sigma 

312 

c 

write(*,260)eta 

313 

314 

315 

c  Compute  the  demoninators  (V  &  H) 

316 

317 

psi  =  0. 

318 

temp  =  (sigma**(-2>)  *  (  -(dsqrtl  1.-(eta**2)  )  )) 

319 

temp2  =  <sigmav**( -2))  *  <  -<dsqrt(  1.-(etav**2)  ) 

320 

denomv  =  temp  /  temp2 

321 

322 

c 

writef*, 1212)denomv 

323 

c1212 

formate  denomv  =  ',f16.8) 

324 

325 

psi  =  pi2 

326 

temp  =  (sigma**( -2))  *  eta 

327 

temp2  =  (sigmah**( -2))  *  etah 

328 

denomh  =  temp  /  temp2 

329 

330 

c 

write(*,1213)denomh 

331 

c1213 

formate  denomh  =  •  r f  16.8) 

332 

333 

c  DIVi 

334 

phi  =  0. 

335 

terml  =  (betal**( -2))  *  gamal 

336 

term2  =  (betar**(-2))  *  gamar 

Line#  Source  Line 
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337 

dlv  =  0.5  *  negflag  *  (  terml  +  term2  )  /  denomv 

338 

339 

c 

OlHi 

340 

phi  =  pi2 

341 

terml  =  (betal**(-2))  *  (-dsqrt(1.-(gamat**2))) 

342 

term2  =  (betar**( -2))  *  (dsqrt(1.-(gamar**2))) 

343 

dlh  =  0.5  *  (  terml  ♦  term2  )  /  denomh 

344 

345 

c 

D2Vi 

346 

d2v  =  (  (l./betal)  +  (l./betar)  )  /  capl 

347 

d2v  =  (  d2v  /  denomv  ) 

348 

349 

c 

D2Hi 

350 

d2h  =  (  (l./betal)  ♦  (l./betar)  )  /  capl 

351 

d2h  =  d2h  /  denomh 

352 

353 

1000 

continue 

354 

355 

356 

C 

03V 

357 

phi  =  0. 

358 

terml  =  (betalv**(-2))  *  gamalv 

359 

term2  =  (betarv**(-2))  *  gamarv 

360 

d3v  =  0.5  *  (  terml  +  term2  ) 

361 

362 

C 

03H 

363 

phi  =  pi  2 

364 

terml  =  (betalh**(-2))  *  ( -dsqrt(1 . -(gamalh**2))) 

365 

term2  =  (betarh**(-2))  *  (dsqrt(1.  (gamarh**2)>) 

366 

d3h  =  0.5  *  (  terml  +  term2  ) 

367 

368 

C 

D4V 

369 

d4v  =  (  (l./betalv)  +  (l./betarv)  )  /  capl 

370 

371 

C 

D4H 

372 

d4h  =  (  (l./betalh)  +  (l./betarh)  )  /  capl 

373 

374 

375 

C 

Print 

out  the  data,  as  well  as  writing  it  to  the  ASCII  file 

376 

377 

C 

write(*,2000) 

378 

c2000 

format(///'  The  0"s  are:',/) 

379 

380 

C 

write(*,2060)i,  dlv,  dlh,  d2v,  d2h 

381 

c2060 

format(/'  EEG  Electrode  #  • , i3, / •  D1V  =  *,f16.5,/ 

382 

c 

+ 

'  D1H  =  1  ,f  16.5,/'  D2V  =  ■ , f  16.5,/'  02H  =  \g16.5) 

383 

384 

c 

write(*,2010)  d3v,  d3h,  d4v,  d4h 

385 

c2010 

format(/'  03V  =  *,f16.5,/'  D3H  =  ',f16.5. 

386 

c 

♦ 

/'  D4V  =  ',f 16.5,/'  D4H  =  «,f16.5,//) 

387 

388 

c 

Scale  0's 

389 

390 

if  (  dabs(dlv)  .gt.  10.  )  dlv  =  (d1v/dabs(d1v))*10. 

391 

if  (  dabs(d2v)  .gt.  10.  )  d2v  =  (d2v/dabs(d2v))*10. 

392 

if  (  dabs(dlh)  .gt.  10.  )  dlh  =  (d1h/dabs(d1h))*10. 
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393  if  (  dabs(d2h)  .gt.  10.  )  d2h  =  (d2h/dabs(d2h))*10. 

394  if  (  dabs(d3v)  .gt.  10.  )  d3v  =  (d3v/dabs(d3v))*10. 

395  if  (  dabs(d4v)  .gt.  10.  )  d4v  =  (d4v/dabs(d4v))*10. 

396  if  (  dabs(d3h)  .gt.  10.  )  d3h  =  (d3h/dabs(d3h))*10. 

397  if  (  dabs(d4h>  .gt.  10.  )  d4h  =  (d4h/dabs(d4h))*10. 

398 

399  max  =  0.000 

400  if  (  dabs(dlv)  .gt.  max  )  max  =  dabs(dlv) 

401  if  (  dabs(d2v)  .gt.  max  )  max  =  dabs(d2v) 

402  if  (  dabs(dlh)  .gt.  max  )  max  =  dabs(dlh) 

403  if  (  dabs(d2h)  .gt.  max  )  max  =  dabs(d2h) 

404  if  (  dabs(d3v)  .gt.  max  )  max  =  dabs(d3v) 

405  if  (  dabs(d4v)  .gt.  max  )  max  =  dabs(d4v) 

406  if  (  dabs(d3h)  .gt.  max  )  max  =  dabs(d3h) 

407  if  (  dabs(d4h)  .gt.  max  )  max  =  dabs(d4h) 

408 

409  scale  =  8192.000  /  max 

410  dlv  =  dlv  *  (-scale) 

411  d2v  =  d2v  *  (-scale) 

412  d3v  =  d3v  *  scale 

413  d4v  =  d4v  *  scale 

414  dlh  =  dlh  *  scale 

415  d2h  =  d2h  *  scale  *  0.01 

416  d3h  =  d3h  *  scale 

417  d4h  =  d4h  *  scale 

418 

419  c  write(*,2100) 

420  c2100  format(///'  The  D"s  for  EEG.C96  are  (in  order):',/) 

421 

422  c  write(*,2070)  scale*d1v,  scale*d2v,  scale*d3v, 

423  c  +  scale*d4v,  scale*d1h,  scale*d2h,  scale*d3h, 

424  c  +  scate*d4h 

425  c2070  format(/8( '  ',f9.0)) 

426 

427  write(*,3333) 

428  3333  format(////////////'  Please  wait...' 

429  +  ////////////) 

430 

431  step  =  3 

432 

433  return 

434  end 

0CAIC  local  Symbols 

Name  Class  Type  Size  Offset 

CCMAX . param  0006 

BEIL . param  000a 

STEP . param  OOOe 

04  H . param  0012 

D3H . param  0016 

04V . param  001a 

03V . param  001  e 

D2H . param  0022 
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DCALC  Local  Symbols 


Name 

Class 

Type 

Size 

Offset 

DIM . 

0026 

D2V . 

002a 

D1V . 

002e 

CIS . 

0032 

ZLINE  . 

.  local 

INTEGER*2 

2 

0002 

IRESP  . 

INTEGER*2 

2 

0004 

BETARV . 

REAL*8 

8 

0006 

SIGMAV . 

REAL*8 

8 

OOOe 

EZLINE . 

INTEGER*2 

2 

0016 

H . 

REAL*8 

8 

0018 

OENOMV . 

.  local 

REAL *8 

8 

0020 

I . 

INTEGER*2 

2 

0028 

J . 

INTEGER*2 

2 

002a 

PI2 . 

REAL*8 

8 

002c 

THETAV . 

REAL*8 

8 

0034 

N . 

INTEGER*2 

2 

003c 

R . 

REAL*8 

6 

003e 

RH . 

REAL*8 

8 

0046 

ETA . 

.  local 

REAL*8 

8 

004 e 

Z . 

REAL*8 

8 

0056 

JEFF . 

REAL*8 

8 

005 e 

CAPL . 

REAL*8 

8 

0066 

PHI . 

REAL*8 

8 

006e 

ETAH . 

REAL*8 

8 

0076 

GAMAL  . 

REAL*8 

8 

007e 

MAX . 

REAL*8 

8 

0086 

I CHAR  . 

INTEGER*2 

2 

008e 

RADOEG . 

REAL*8 

8 

0090 

SCALE  . 

REAL*8 

8 

0098 

TEMP2  . 

REAL*8 

8 

OOaO 

RV . 

REAL*8 

8 

00a8 

GAMAR  . 

REAL*8 

8 

OObO 

8ETAL  . 

REAL*8 

8 

00b8 

TERM1  . 

REAL*8 

8 

OOcO 

DELTA  . 

REAL*8 

8 

00c8 

TERM2  . 

REAL*8 

8 

OOdO 

GAMALH . 

REAL*8 

8 

00d8 

PS1 . 

REAL*8 

8 

OOeO 

NEGFLA . 

REAL*8 

8 

00e8 

STEP1  . 

INTEGER*2 

2 

OOfO 

BETAR  . 

REAL*8 

8 

00f2 

GAMAR H . 

REAL*8 

8 

OOfa 

BETALH . 

REAL*8 

8 

0102 

ETAV . 

REAL*8 

8 

010a 

SIGMA  . 

REAL*8 

8 

0112 

TEMP . 

REAL*8 

8 

011a 

BETARH . 

REAL‘8 

8 

0122 

THETA  . 

REAL*8 

8 

012a 

GAMAL V . 

REAL*8 

8 

0132 

SIGMAH . 

REAL*8 

8 

013a 

RESP . 

CHAR*1 

1 

0142 
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DCALC  Local  Symbols 


Name 

Class 

Type 

Size 

Offset 

DENOMH . 

REAL*8 

8 

0144 

THETAH . 

.  local 

R£AL*8 

8 

014c 

GAMARV . 

REAL*8 

8 

0154 

BETALV . 

.  local 

REAL*8 

8 

015c 

435 

437 

438 

C  '  ' 

subroutine  error(bell) 

439 

440 

characters  resp,  bell 

441 

442 

write(*,100)bell 

443 

100 

format(///12x,a1, 'Invalid  Response  !'//,13x. 

444 

♦'Type  ENTER  to  continue. . . ' ) 

445 

read(*,20)resp 

446 

20 

format(bn,a1) 

447 

448 

return 

449 

end 

ERROR  Local  Symbols 

Name 

Class 

BELL . 

RESP . 

.  par am 
.  local 

Global  Symbols 

Name 

Class 

0CALC  . 

ERROR  . 

.  FSUBRT 
.  FSUBRT 

Type 

Size  Offset 

CHAR*! 

0006 

1  0164 

Type 

Size 

Offset 

»** 

*** 

0000 

*** 

*** 

1359 

Code  size  =  1359  (4953) 
Data  size  =  0206  (518) 
Bss  size  =  0165  (357) 

No  errors  detected 


SGN 


Line# 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


PAGE  1 
03-22-90 
17:43:57 

Source  Line  Microsoft  FORTRAN  Optimizing  Compiler  Version  4.00 

{DECLARE 

SLARGE 

subroutine  sgm(magsuni,  orient) 
c  Jeffrey  C.  Sigl  December  26,  1989 


c  Data  Structures 

character*6 

character*! 

character*! 

character*! 

integer*2 

real *4 

complex*8 

integer 

c  functions 


resp 

bite(2) 

cappa,  cappb,  cappc,  cappd,  cappe,  cappf 
chr(6),  one,  sp,  cr,  retran,  oktran 
in,  istart,  i,  j,  k,  irec,  orient 
eeg(4096),  magsum 
sg(4096),  sgtemp 
tmprsp 


integer  ichar 

real*4  float,  cabs 

complex*8  cnplx 

c  Data  Relations 


equivalence  (resp.chr) 


data  cappa , cappb, cappc , cappd , cappe, cappf / ' A ' ,,B,,,C,,,D,,,E,,,F,/ 
data  one  /49/ 

data  sp,  retran,  oktran  /8#140,  0,  1/ 


open  (9, fi le=' scrtch.dat ' ,status=lneHl ,access=' sequential ', 
+form=' formatted' ) 


c  Load  EEG  data  -  Read  120  words 
j  =  2018 


do  33  i  =  1,  60 

17  read(11,105)in 

105  format(i6) 

c  write(*,105)in 

orient  =  in  -  ((in/10)  *  10) 
in  =  in/10 

if  (  orient  .eq.  1  )  then 
in  =  -in 

endif 


eeg(j)  =  float (in) 
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57 

j=j+1 

58 

59 

33 

continue 

60 

61 

c  Extend  ends 

62 

63 

do  20  i  =  1,  2017 

64 

eeg(i)  =  eeg(2018) 

65 

20 

cont  i  nue 

66 

67 

do  30  i  =  2078,  4096 

68 

eeg(i)  =  eeg(2077) 

69 

30 

cont i nue 

70 

71 

c  Window  &  FFT  eeg 

72 

73 

call  sgwin(eeg,  eeg,  9,  0.001,  12) 

74 

75 

do  60  i  =  1,  4096 

76 

sg(i)  =  cmplx<  eeg(i),  0.  ) 

77 

60 

cont i nue 

78 

79 

call  cfft(  sg,  12,  0,  1.0  ) 

80 

81 

do  2323  i  =  1,  4096 

82 

write<9,*)sg(i) 

83 

2323 

cont i nue 

84 

85 

86 

87 

c  Load 

EOG  data  •  Read  120  words 

88 

89 

j  =  2018 

90 

do  133  i  =  1,  60 

91 

92 

117 

read(11,105)in 

93 

94 

orient  =  in  -  ((in/10)  *  10) 

95 

in  =  in/10 

96 

if  (  orient  .eq.  1  )  then 

97 

in  =  -  in 

98 

endif 

99 

100 

eeg(j)  =  float(in) 

101 

102 

j=M 

103 

133 

continue 

104 

105 

c  Extend  ends 

106 

107 

do  120  i  =  1,  2017 

108 

eeg(i)  =  eeg(2018) 

109 

120 

cont i nue 

110 

111 

112 


do  130  i  =  2078,  4096 

eeg(i)  =  eeg(2077) 


Line#  Source  Line 
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113 

130 

cont i nue 

114 

115 

c  Window  &  FFT  eog 

116 

117 

call  sgwin(eeg,  eeg,  9,  0.001,  12) 

118 

119 

do  170  i=1,  4096 

120 

sg(i)  =  cmplxC  eeg(i),  0.  ) 

121 

170 

continue 

122 

123 

call  cfft<  sg,  12,  0,  1.0  ) 

124 

125 

126 

c  Divide  &  compute  the  average  magnitude 

127 

128 

magsum  =  0. 

129 

rewind  9 

130 

do  1013  i  =  6,  45 

131 

read(9,*)sgtemp 

132 

eeg(i)  =  100.  *  cabs(  sgtemp/sg( 

133 

magsum  =  magsum  +  eeg(i) 

134 

135 

c 

write(*,2345)i ,eeg(i ), magsum 

136 

c2345 

formate  ',i6,2('  ',f15.3,)> 

137 

138 

1013 

continue 

139 

magsum  =  magsum  /  40. 

140 

141 

9999 

clos' C9,status='delete' ) 

142 

143 

return 

144 

end 

SGM  Local  Symbols 


Name  Class  Type  Size  Offset 

ORIENT . param  0006 

MAGSUM . param  000a 

EEG  .  local  REAL*4  16384  0000 

I . local  1NTEGER*2  2  0002 

J . local  INTEGER*2  2  0004 

K . local  1NTEGER*2  2  0006 

SGTEMP . local  C0MPLEX*8  8  0008 

CAPPA . local  CHAR*1  1  000a 

CAPPB . local  CHAR*1  1  000b 

CAPPC . local  CHAR*1  1  000c 

CAPPD . local  CHAR*1  1  OOOd 

CAPPE . local  CHAR*1  1  OOOe 

CAPPF . local  CHAR*1  1  OOOf 

ONE . local  CHAR*1  1  0010 

CR . local  CHAR*!  1  0010 

SP . local  CHAR*1  1  0011 

RETRAN . local  CHAR*1  1  0012 

ISTART . local  INTEGER*2  2  0012 
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SGM  Local  Symbols 


Name 

Class 

Type 

Size 

Offset 

OUTRAN . 

.  local 

CHAR*1 

1 

0013 

IN . 

.  local 

INTEGER*2 

2 

0014 

IREC . 

.  local 

INTEGER*2 

2 

0016 

TMPRSP . 

.  local 

INTEGERS 

4 

0018 

ICHAR  . 

.  local 

INTEGER*4 

4 

001c 

RESP . 

.  local 

CHAR*6 

6 

0020 

SG . 

.  local 

COMPLEX*8 

32768 

4000 

BITE . 

.  local 

CHAR*1 

2 

cOOO 

CHR . 

.  local 

CHAR*1 

6 

0020 

Global  Symbols 

Name 

Class 

Type 

Size 

Offset 

CFFT . 

.  .  .  extern 

*** 

*** 

*** 

SGM . 

.  .  .  FSUBRT 

*** 

*** 

0000 

SGWIN  . 

.  .  .  extern 

*** 

*** 

*** 

Code  size  =  04b0  (1200) 
Data  size  =  0060  (96) 
Bss  size  =  0026  (38) 


No  errors  detected 


FTUNE2 
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1 

(LARGE 

2 

(declare 

3 

c 

f tune. for 

4 

c 

5 

A 

c 

Steven  M.  Falk 

O 

7 

subroutine  ftune(d1v,d2v,d1h,d2h,d3v,d4v,d3h,d4h, vmag,hmag,gain. 

8 

9 

+ccmax) 

10 

c  Data 

Structures 

11 

12 

character^ 

cls<4),  chr(6),  cbell 

13 

characters 

resp,  sp,  retran,  oktran,  cr 

14 

integer 

i,  j,  n,  k,  1,  tmprsp,  in,  mod 

15 

real *4 

dlv,  d2v,  dlh,  d2h 

16 

real *4 

d3v,  d4v,  d3h,  d4h,  ccmax 

17 

real*4 

cc,  ccrawv,  ccrawh,  ccold,  deltal. 

delta2,  ccrti 

18 

integer 

ichar 

19 

real*4 

float,  cabs,  real,  abs 

20 

complex*8 

cmplx 

21 

character 

yes,  yess 

22 

integer*2 

num,  iorder,  icnst2,  ilim,  orent, orient 

23 

real *4 

eogvu(512),  eogvl(512),  eogh(512). 

raweeg(512) 

24 

real*4 

eohvu(512),  eohvl(512),  eohh(512). 

raweeh(512) 

25 

real*4 

alpha(512),  temp(512) 

26 

real *4 

time,  alphal,  alpha2,  const,  intrvl 

,  ceeg 

27 

real*4 

vmag,  hmag,  gain,  raddeg,  order,  constl,  pi2 

28 

complex*8 

comp1(512),  comp2(512) 

29 

double  precision  scale,  max,  dabs 

30 

31 

equivalence 

(resp.chr) 

32 

33 

c  Data 

Initialization 

34 

35 

pi2  =  1 . 57079632680D0 

36 

raddeg  =  57.2957795131D0 

37 

data  yes. 

yess  t'y' , 'V'/ 

38 

data  sp,  retran,  oktran  /8#140,  0,  1/ 

39 

data  cbell  /  8#7/ 

40 

41 

c 

deltal  =  0.5 

*  dlv 

42 

deltal  =  100. 

43 

delta2  =  0.5 

*  d2v 

44 

do  401  i  =  1, 

512 

45 

46 

407 

read(1 1,402) in 

47 

402 

formatl i6) 

48 

49 

orent 

=  in  -  <( in/10)  •  10) 

50 

in  = 

in/10 

51 

if  ( 

orent  .eq.  1  )  then 

52 

in  =  -in 

53 

endif 

54 

55 

56 

c 

if  (abs(floatdn))  .gt.  470.)  then 
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57 

c 

gain  =  -9999. 

58 

c 

write  (*,4433) 

59 

C4433 

formate  Data  is  saturated.  Interrogate  aga 

60 

c 

goto  2222 

61 

c 

end  if 

62 

63 

raweeg(i)  =  float(in) 

64 

65 

401 

cont i nue 

66 

67 

do  1401 

i  =  1,  512 

68 

69 

1407 

read(1 1,1402) in 

70 

1402 

format(i6) 

71 

72 

orent  =  in  -  ((in/10)  *  10) 

73 

in  =  in/10 

74 

if  (  orent  .eq.  1  )  then 

75 

in  =  -  in 

76 

endif 

77 

78 

c 

if  (abs(f loat(in))  .gt.  470.)  then 

79 

c 

gain  =  -9999. 

80 

c 

write  (*,4433) 

81 

c 

goto  2222 

82 

c 

endif 

83 

84 

raweeg(i)  =  float(in) 

85 

86 

1401 

cont i nue 

87 

88 

do  411 

i  =  1,  512 

89 

90 

417 

read(11,412)in 

91 

412 

format( i6) 

92 

93 

orent  =  in  -  ((in/10)  *  10) 

94 

in  =  in/10 

95 

if  (  orent  .eq.  1  )  then 

96 

in  =  -in 

97 

endif 

98 

99 

c 

if  (abs(f loat(in))  .gt.  470.)  then 

100 

c 

gain  =  -9999. 

101 

c 

write  (*,4433) 

102 

c 

goto  2222 

103 

c 

endif 

104 

105 

eogh(i)  =  float(in) 

106 

107 

411 

continue 

108 

109 

do  1411 

i  =  1,  512 

110 

111 

1417 

read(11,1412)in 

112 

1412 

format(i6) 
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113 

114 

115 

116 

117 

118 

119 

120  c 

121  c 

122  c 

123  c 

124  c 

125 

126 

127 

128  1411 

129 

130 

131 

132  427 

133  422 

134 

135 

136 

137 

138 

139 

140 

141  c 

142  c 

143  c 

144  c 

145  c 

146 

147 

148 

149  421 

150 

151 

152 

153  1427 

154  1422 

155 

156 

157 

158 

159 

160 
161 

162  c 

163  c 

164  c 

165  c 

166  c 

167 

168 


orent  =  in  -  ((in/10)  *  10) 

in  =  in/10 

if  (  orent  -eq.  1  )  then 
in  =  -  in 

endif 

if  (abs(float(in))  .gt.  470.)  then 
gain  =  -9999. 

Mrite  (*,4433) 
goto  2222 

endif 

eogh(i)  =  float(in) 


continue 

do  421  i  =  1,  512 

readd  1,422)  in 

format(i6) 

orent  =  in  -  ((in/10)  *  10) 

in  =  in/10 

if  (  orent  .eq.  1  )  then 
in  =  -in 

endif 

if  (abs(f loat(in))  .gt.  470.)  then 
gain  =  -9999. 
write  (*,4433) 
goto  2222 

endif 

eogvu(i)  =  float(in) 


cont i nue 

do  1421  i  =  1,  512 

read(1 1,1422) in 

format(i6) 

orent  =  in  -  ((in/10)  *  10) 

in  *  in/10 

if  (  orent  .eq.  1  )  then 
in  =  -  in 

endif 

if  (abs(f loat(in))  .gt.  470.)  then 
gain  =  -9999. 
write  (*,4433) 
goto  2222 

endif 


eogvu(i)  =  float(in) 


Line#  Source  Line 
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169 

170  1421 

171 

172 

173 

174  437 

175  432 

176 

177 

178 

179 

180 
181 
182 

183  c 

184  c 

185  c 

186  c 

187  c 

188 

189 

190 

191  431 

192 

193 

194 

195  1437 

196  1432 

197 

198 

199 

200 
201 
202 

203 

204  c 

205  c 

206  c 

207  c 

208  c 

209 

210 
211 

212  1431 

213 

214 

215  c 

216 

217  11 

218 

219 

220 
221 
222 

223 

224 


continue 

do  431  i  =  1,  512 

read(1 1,432) in 
format( i6) 

orent  =  in  -  ((in/10)  *  10) 
in  =  in/10 

if  (  orent  .eq.  1  )  then 
in  =  -in 

end  if 

if  (abs(f loat(in))  .gt.  470.)  then 
gain  =  -9999. 
write  (*,4433) 
goto  2222 

endif 

eogvl(i)  =  float(in) 

continue 

do  1431  i  =  1,  512 

read(11,1432)in 
format( i6) 

orent  =  in  -  ((in/10)  *  10) 
in  =  in/10 

if  (  orent  .eq.  1  )  then 
in  =  -in 

endif 

if  (abs(f loat( in))  .gt.  470.)  then 
gain  =  -9999. 
write  (*,4433) 
goto  2222 

endif 

eogvl(i)  =  float(in) 

continue 


fine  tuning  algorithm 

cont i nue 

call  ccf(eogvu,raweeg,ccrawv) 
call  ccf(eogh,raweeg,ccrawh) 

call  Crct(d1v,d2v,d1h,d2h,d3v,d4v,d3h,d4h, vmag,hmag,gain, 
*eo gvu,eogvl ,eogh, raweeg, temp) 


Source  Line 
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225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 


call  ccf (eogvu, temp, ccold) 

if  (abs(ccold)  .gt.  abs(ccrawv))  then 
write(*,1109)cbell 

1109  formate  Data  Entry  Error! *,a1) 

gain  =  -9999. 
goto  2222 

endif 

if  (abs(ccold)  .gt.  0.8)  then 
write(*,1108)cbell 

1108  formate  Measurement  Error!',a1) 

gain  =  -9999. 
goto  2222 

endif 

call  ccf(eogh,temp,ccold) 

if  (abs(ccold)  .gt.  abs(ccrawh))  then 
write(*,1109)cbell 
gain  =  -9999. 
goto  2222 

endif 

if  (abs(ccold)  .gt.  0.8)  then 
wri te(*, 1 108)cbel l 
gain  =  -9999. 
goto  2222 

endif 

2431  wri te(*, 2460) 

2460  formate  Fine  tuning  geometric  VERTICAL  parameters  of  model...*) 
wri te(*, 2462) 

2462  formate  Iteration  01  02  CCF') 

call  ccf (eogvu, temp, ccold) 
l  =  0 

write(*,*)l,d1v,d2v,ccold 

2500  1=1+1 


dlv  =  dlv  +  deltal 

call  Crct(d1v,d2v,d1h,d2h,d3v,d4v,d3h,d4h, vmag,hmag,gain, 
♦eogvu, eogvl ,eogh, raweeg, temp) 
call  ccf (eogvu, temp, ce) 

ccrtio  =  ce  /  ccold 
if  (ccrtio  .It.  0)  then 

deltal  =  deltal  *  (-.5) 
goto  2501 

endif 

if  (abs(cc)  .gt.  abs(ccotd))  then 
deltal  =  deltal  *  (-1.) 
goto  2501 


endif 


Line#  Source  Line 
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281 

282  2501 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299  2511 

300 

301 

302 

303  2512 

304 

305  2513 

306 

307 

308  2518 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333  2600 

334 

335 

336 


ccold  =  cc 

d2v  =  d2v  +  delta2 

call  crct (dl v,d2v, dlh, d2h,d3v,d4v,d3h,d4h,vmag,hmag, gain, 

+eogvu,eogvl ,eogh, raweeg, temp) 
call  ccf (eogvu, temp, cc) 

ccrtio  =  cc  /  ccold 
if  (ccrtio  .It.  0)  then 

delta2  =  delta2  *  (-.5) 
goto  2511 

endif 

if  (abs(cc)  .gt.  abs(ccold))  then 
delta2  =  delta2  *  (-1.) 
goto  2511 

endif 

write(*,*)l,d1v,d2v,cc 
j  =  mod(l,10) 
if  (j  -eq.  0)  then 

wri te(*,2512) 

formate  Continue  iterations  ?  (Y/N)  *,\) 

read(*,2513)resp 

format(al) 

if  (resp  .ne.  'Y')  then 
write<*,2518) 

formatf  '  Fine  tuning  geometric  HORIZONTAL  parameters', 

+  '  of  model . . . 1 ) 

urite(*,2462) 

call  ccf (eogh, temp, ccold) 

1=0 

write(*,*)l ,d1h,d2h, ccold 
deltal  =  0.5  *  dlh 
deltaZ  =  0.5  *  d2h 
goto  2600 

endif 

endif 

ccold  =  cc 

if  (abs(cc)  .It.  ccmax)  then 
writer*, 2518) 
wri te(*, 2462) 
call  ccf(eogh, temp, ccold) 

1=0 

write(*,*)l,d1h,d2h,ccold 
deltal  =  0.5  *  dlh 
delta2  =  0.5  *  d2h 
goto  2600 

endif 
goto  2500 

1  =  1  +  1 

dlh  =  dlh  +  deltal 

cal l  crctfdlv, d2v, dlh, d2h,d3v,d4v,d3h,d4h, vmeg.hmag, gain. 


tine#  Source  Line 
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337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 


♦eogvu.eogvl ,eogh,raweeg, temp) 
call  ccf(eogh, temp,cc) 

ccrtio  =  cc  /  ccold 
if  (ccrtio  .It.  0)  then 

deltal  =  deltal  *  (-.5) 
goto  2601 

end  if 

if  (abs(cc)  .gt.  abs(ccold))  then 
deltal  =  deltal  *  (-1.) 
goto  2601 

end  if 

2601  ccold  =  cc 

d2b  =  d2h  +  delta2 

cal l  Crct(d1v,d2v,d1h,d2h,d3v,d4v,d3h,d4h,vmag,hmag,gain, 
+eogvu,eogvl ,eogh,raweeg, temp) 
call  ccf(eogh,temp,cc) 

ccrtio  =  cc  /  ccold 
if  (ccrtio  .It.  0)  then 

delta2  =  delta2  *  (-.5) 
goto  2611 

endif 

if  (abs(cc)  .gt.  abs(ccold))  then 
delta2  =  delta2  *  (-1.) 
goto  2611 

endif 

2611  write(*,*)l,d1h,d2h.cc 
j  =  mod(l,10) 
if  (j  .eq.  0)  then 

write(*,2512) 
read(*,2513)resp 
if  (resp  .ne.  'Y 1 )  then 
goto  2800 

endif 

endif 

ccold  =  cc 

if  (abs(cc)  .It.  ccmax)  then 
goto  2800 

endif 
goto  2600 

2800  cont i nue 


open( 1 2 , f i l e= 1 correeg.dat ' ) 
do  2803  i  =64 ,384 

wri te( 12,*) i , temp( i ), raweeg( i ), eogvu( i ), eogh( i ) 

continue 

close(12) 


2803 

c  Scale  0's 


4.00 
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393  max  =  0.000 

394  if  (  dabs(dlv)  .gt.  max  )  max  =  dabs(dlv) 

395  if  (  dabs(d2v)  .gt.  max  )  max  =  dabs(d2v) 

396  if  (  dabs(dlh)  .gt.  max  )  max  =  dabs(dlh) 

397  if  (  dabs(d2h)  .gt.  max  )  max  =  dabs(d2h) 

398  if  (  dabs(d3v)  .gt.  max  )  max  =  dabs(d3v) 

399  if  (  dabs(d4v)  .gt.  max  )  max  =  dabs(d4v) 

400  if  (  dabs(d3h)  .gt.  max  )  max  =  dabs(d3h) 

401  if  (  dabs(d4h)  .gt.  max  )  max  =  dabs(d4h) 

402 

403  scale  =  8192.000  /  max 

404  dlv  =  dlv  *  scale 

405  d2v  =  d2v  *  scale 

406  d3v  =  d3v  *  scale 

407  d4v  =  d4v  *  scale 

408  dlh  =  dlh  *  scale 

409  d2h  =  d2h  *  scale 

410  d3h  =  d3h  *  scale 

411  d4h  =  d4h  *  scale 

412 

413  2222  return 

414  end 

FTUNE  Local  Symbols 

Name  Class  Type  Size  Offset 

CCMAX . param  0006 

GAIN . param  000a 

HMAG . param  OOOe 

VMAG . param  0012 

D4H . param  0016 

D3H . param  001a 

04V . param  OOle 

03V . param  0022 

D2H . param  0026 

01 H . param  002a 

02V . param  002e 

01V . param  0032 

EOGVU . local  REALM  2048  0000 

CMPLX . local  C0MPLEX*8  8  0002 

CCRTIO . local  REALM  4  000a 

YES . local  CHAR*1  1  OOOd 

YESS . local  CHAR*1  1  OOOe 

l  ORDER . local  INTEGER*2  2  OOOe 

SP . local  CHARM  1  OOOf 

RETRAN . local  CHARM  1  0010 

CCRAWV . local  REALM  4  0010 

OKTRAN . local  CHARM  1  0011 

CBELL . local  CHARM  1  0012 

CC . local  REALM  4  0014 

CONST . local  REALM  4  0018 

ORENT . local  !NTEGER*2  2  001c 

l . local  INTEGERM  4  OOle 

J . local  INTEGERM  4  0022 
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FTUNE  Local  Symbols 


Name  Class  Type  Size  Offset 

P 12 . local  REAL*4  4  0026 

K . local  INTEGER*4  4  002a 

L . local  INTEGER*4  4  002e 

N . local  INTEGER*4  4  0032 

ORIENT . local  INTEGER*2  2  0036 

CEEG . local  REAL*4  4  0038 

CR . local  CHAR*1  1  003c 

ALPHA1 . local  REAL  *4  4  003e 

IN . local  INTEGER*4  4  0042 

ALPHA2 . local  REAL*4  4  0046 

CABS . local  REAL*4  4  004a 

DELTA1 . local  REAL*4  4  004e 

0ELTA2 . local  REAL*4  4  0052 

INTRVL . local  REAL*4  4  0056 

REAL . local  REAL*4  4  005a 

CCOLD . local  REAL*4  4  005e 

MAX . local  REAL*8  8  0062 

TMPRSP . local  INTEGER*4  4  006a 

RAODEG . local  REAL*4  4  006e 

I  CHAR . local  INTEGERS  4  0072 

SCALE . local  REAL*8  8  0076 

ILIM . local  INTEGER*2  2  007e 

TIME . local  REAL*4  4  0080 

NUM . local  INTEGER*2  2  0084 

ICNST2 . local  INTEGERS  2  0086 

CONST1 . local  REAL*4  4  0088 

CCRAUH . local  REAL*4  4  008c 

ORDER . local  REAL*4  4  0090 

EOHVU . local  REAL*4  2048  0800 

CHR . local  CHAR*1  6  1000 

C0MP1 . local  COMPLEX*8  4096  1006 

C0MP2 . local  C0MPLEX*8  4096  2006 

CLS . local  CHAR*1  4  3006 

EOGH . local  REAL*4  2048  300a 

EOHH . local  REAL*4  2048  380a 

ALPHA . local  REAL*4  2048  400a 

TEMP . local  REAL*4  2048  480a 

RAWEEG . local  REAL*4  2048  500a 

RAWEEH . local  REAL*4  204B  580a 

EOGVL . local  REAL*4  2048  600a 

EOHVL . local  REAL*4  2048  680a 

RESP . equiv  CHAR*1  1  1000 


415 

417 

418 

subroutine  ccf (ueeg,ceeg,cc) 

419 

420 

real *4 

ccfcov,sho 

421 

real *4 

ueeg(512),eeeg(512) 

4  22 

real *8 

a(512,6),b(6) 
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423  integer*2  i  i ,  j  j 

424 

425 

426  do  10  i i =1 ,512 

427  a(ii,1)  =  ueeg( i i ) 

428  a(ii,2)  =  ceeg(ii) 

429  a(ii,3)  =  a( i i , 1 )*a< i i ,2) 

430  a(ii,4)  =  a( i i ,  1 )+a( ii ,2) 

431  8(11,5)  =  a( i i , 1 )*a( i i , 1 ) 

432  a(ii,6)  =  a(ii,2)*a(ii,2) 

433  10  continue 

434 

435  do  20  ii=1,6 

436  b( i i)  =  0.0 

437  20  continue 

438 

439  do  30  ii =64, 448 

440  do  31  jj=1,6 

441  b(jj)  =  bt  j  j )  +  a(ii,jj) 

442  31  continue 

443  30  continue 

444  do  40  i i=1 ,6 

445  b(ii)  =  b(ii)  /  384.0 

446  40  continue 

447 

448  cov  =  b(3)  -  (b(1)  *  b(2)) 

449  sho  =  (b(5)-(b(1)*b(1)))  *  <b<6)-<b(2)*b<2)» 

450  if  (  sho  .le.  0.0  )  then 

451  sho  =  0.0 

452  endif 

453  sho  =  sho**. 5 

454  if  (  sho  -eq.  0.0  )  then 

455  cc  =  1.0 

456  goto  75 

457  endif 

458  cc  =  cov  /  sho 

459 

460  75  continue 

461 

462  return 

463  end 


CCF  local  Symbols 

Name  Class  Type  Size  Offset 

CC . param  0006 

CEEG . param  000a 

UEEG . param  000e 

•1 . local  INTEGER*2  2  0094 

JJ . local  INTEGER*2  2  0096 

COV . local  REAL*4  4  0098 

SHO . local  REAl*4  4  009c 

*  .  local  REAl*8  24576  7142 

*  . local  REAL*8  48  d142 


Microsoft  FORTRAN  Optimizing  Compiler  Version  4.00 


464 

465 

466  . . 

467 

468  subroutine  Crct(d1v,d2v,d1h,d2h,d3v,d4v,d3h,d4h, vmag,hmag,gain, 

469  +eogvu,eogvl,eogh,raweeg,temp) 

470 


471 

integer 

i,  j,  n,  k,  tmprsp,  in 

472 

real *4 

dlv,  d2v,  dlh,  d2h 

473 

real *4 

d3v,  d4v,  d3h,  d4h 

474 

real*4 

float,  cabs,  real 

475 

compl ex*8 

cmplx 

476 

real *4 

eogvu(512),  eogvl(512),  eogh(512),  raweeg(512) 

477 

real*4 

teogvu(52),  teogvl(512>,  teogh(512) 

478 

real *4 

alpha(512),  temp(512) 

479 

real *4 

time,  alphal,  alpha2,  const,  intrvl,  ceeg 

480 

real*4 

vmag,  hmag,  gain,  raddeg,  order,  constl,  pi2 

481 

complex*8 

comp1(512),  coirp2(512) 

482 

483 

484  c  Window  &  FFT  EOG-VU 

485 

486  call  sgwin(eogvu,  teogvu,  9,  0.001,  9) 

487  do  60  i  =  1,  512 

488  compl(i)  =  cmplx<  teogvu(i),  0.  ) 

489  60  cont i nue 

490  call  cfft(  compl,  9,  0,  1.0  ) 

491 

492 

493  c  Window  &  FFT  EOG-VL 

494 

495  call  sgwinfeogvl,  teogvl,  9,  0.001,  9) 

496  do  70  i  =  1,  512 

497  conp2(i)  =  cmplx(  teogvl(i),  0.  ) 

498  70  continue 

499  call  cfft<  comp2,  9,  0,  1.0  ) 

500 

501 

502  c  Compote  Alpha 

503 

504  do  80  i  =  1,  512 

505  tempi i)  =  (gain  *  cabs(  compl(i)  /  comp2(i)  ))  /  100. 

506 

507  if  (temp(i)  .It.  1.)  then 

508  temp( i )  =  1 . 

509  elseif  (temp(i)  .gt.  4.)  then 

510  tempi i)  =  4. 

511  endif 

512 

513  alpha(i)  =  temp(i) 

514  80  continue 

515 

516  goto  8989 

517 

518  c  Moving  Average  Filter  of  Alpha  (7  point) 


Line# 


Source  Line 
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519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 


alpha<1)  =  ((4.*temp( 1 ))  +  temp(2)  +  temp(3)  +  temp(4>)  /  7. 
alpha{2)  =  ((3.*temp(D)  +  temp(2)  +  temp(3)  +  temp(4) 

+  +  temp<5))  /  7. 

alpha<3)  =  ((2.*temp(1))  +  temp(2)  +  temp<3)  +  temp<4) 

+  +  temp(5)  +  temp(6))  /  7. 

do  90  i  =  4,  509 

alpha(i)  =  (temp(i-3)  +  temp(i-2)  +  temp(i-l)  +  temp(i)  ♦ 
+  temp(i+1)  +  tempCi+2)  +  temp(i+3)>  /  7. 

90  cont i nue 

alpha(510)  =  (temp<507)  +  temp(508)  +  temp(509)  + 

+  temp(510)  +  temp(511)  +  <2.*temp<512>))  /  7. 

alpha(511)  =  (temp(508)  +  temp(509)  +  temp(510>  + 

+  temp(511)  +  (3.*temp(512)))  /  7. 

alpha(512)  =  <temp(509)  +  temp(510)  ♦  temp(511)  + 

+  (4.*temp(512)>)  /  7. 


c  Window  &  FFT  EOG-H 

8989  call  sgwin(eogh,  teogh,  9,  0.001,  9) 

do  100  i  =  1,  512 

comp2(i)  =  cmplx.  teogh(i),  0.  ) 

100  continue 

call  cfft(  comp2,  9,  0,  1.0  ) 


c  = 
c 

c  = 


Correction 


do  1000  i  =  1,  512 

alphal  =  alpha( i )  +  1 . 
alpha2  =  alpha( i )  -  1 . 

c  Vertical  Component 

const  =  (d1v*alpha1  +  d2v*alpha2)  /  (d3v*atpha1  +  d4v*alpha2) 
compl(i)  =  const  *  vmag  *  compl(i)  /  100. 

c  Horizontal  Component 

const  *  (d1h*alpha1  +  d2h*alpha2)  /  (d3h*alpha1  +  d4h*alpha2) 
compl(i)  =  compl(i)  +  (const  *  hmag  *  comp2(i>)  /  100. 

tempi i )  =  1 . 


1000  continue 


c  Inverse  Transform 


call  cift(  comp 1,  9,  0,  1.0  ) 


Line#  Source  Line 
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575  c  Dewindow  Corrector  &  Subtract  from  Raw  EEG 

576 

577  call  sgwin(temp,  temp,  9,  0.001,  9) 

578  do  1100  i  =  1,  512 

579  temp(i)  =  real(comp1(i)>  /  temp(i) 

580  tempi i)  =  raweeg(i)  -  tempi i) 

581 

582  1100  continue 

583 

584  return 

585  end 

CRCT  Local  Symbols 


Name  Class  Type  Size  Offset 

TEMP . param  0006 

RAWEEG . param  000a 

EOGH . param  OOOe 

EOGVL . param  0012 

EOGVU . param  0016 

GAIN . param  001a 

HMAG . param  001  e 

VMAG . param  0022 

D4H . param  0026 

D3H . param  002a 

OAV . param  002e 

03  V . param  0032 

D2H . param  0036 

01 H . param  003a 

02V . param  003e 

01V . param  0042 

ALPHA . local  REAL*4  2048  0000 

CONST . local  REAL*4  4  OOaO 

I . local  1NTEGER*4  4  00a4 

J . local  INTEGER*4  4  00a8 

PI2 . local  REAL*4  4  OOac 

K . local  INTEGER*4  4  OObO 

N . local  INTEGER*4  4  00b4 

CEEG . local  REAL*4  4  00b8 

ALPHA1 . local  REAL*4  4  OObc 

IN . local  INTEGER*4  4  OOcO 

ALPHA2 . local  REAL*4  4  00c4 

INTRVL . local  REAL*4  4  00c8 

TMPRSP . local  INTEGER*4  4  OOcc 

RAODEG . local  REAL*4  4  OOdO 

TIME . local  REAL*4  4  00d4 

FLOAT . local  REAL*4  4  00d8 

CONST  1 . local  REAL*4  4  OOdc 

ORDER . local  REAL*4  4  OOeO 

TEOGH . local  REAL  *4  2048  0800 

TEOGVL . local  REAL*4  2048  d172 

TEOGVU . local  REAL*4  208  d972 

COMP1 . local  COMPLEX*8  4096  da42 

C0MP2 . local  C0MPLEX*8  4096  ea42 
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Global  Symbols 


Name 

Class 

CCF . 

CFFT . 

Cl  FT . 

CRCT . 

FTUNE  . 

.  FSUBRT 

SGUIN  . 

Type 

Size 

Offset 

*** 

★  ★★ 

1 1  f  2 

*#* 

★  ** 

*** 

*#* 

★  Hr* 

★  ★★ 

★  *★ 

★*★ 

1590 

*** 

★  ★★ 

0000 

*** 

*** 

*** 

Code  size  =  2282  (8834) 
Data  size  =  01b3  (435) 
Bss  size  =  00e4  (228) 


No  errors  detected 


on  4.00 


